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comparisons  shoved  that.  In  general,  the  calculated  base  pressures  are  greater 
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about  the  applicability  of  a  single  stream  recompression  criterion  to  the 
solution  of  the  two  stream  problem. 


SUMMARY 


The  base  pressure  computer  program  of  Addy  [4]  has  been  modified  to 
Include  the  effects  of  the  upstream  boundary  layers  upon  the  free  shear  layers 
and  to  Include  the  ONERA  angular  recompresslon  criterion,  as  Implemented  by 
Wagner  and  White  [12]  for  the  wake  closure  condition.  After  the  program 
modifications  were  completed,  an  extensive  set  of  calculations  were  made  to 
make  comparisons  with  the  calculations  of  Wagner  and  White  and  to  make  com¬ 
parisons  with  a  broad  spectrum  of  experimental  data.  The  comparisons  with  the 
calculations  of  Wagner  and  White  showed  very  good  agreement  of  the  two 
calculations.  Comparisons  with  the  experimental  data  showed  that,  in  general, 
the  calculated  base  pressures  are  greater  than  experimental  base  pressures. 

The  Office  National  d 'Etudes  et  de  Recherches  Aerospatlales  (ONERA)  angu¬ 
lar  criterion  was  originally  developed  for  two-dimensional  flow  through  a 
correlation  of  a  comprehensive  set  of  experimental  data.  Later,  a  factor  to 
convert  the  two-dimensional  criterion  to  an  axisymmetrlc  criterion  was  found 
through  the  correlation  of  a  set  of  experimental  data  for  axls3rmmetric  flow. 
Both  of  these  correlations  were  developed  using  data  from  a  single-stream 
impinging  upon  a  wall.  The  analysis  of  this  data  to  define  the  critical  angle 
is  unique,  and  in  turn,  its  application  to  a  single-stream  problem  is 
unambiguous.  However,  when  this  criterion  is  applied  to  the  solution  of  the 
two-stream  problem,  an  ambiguity,  either  real  or  apparent,  arises.  In 
general,  the  criterion  yields  different  reattachment  pressures  for  the  two 
streams.  Since  application  of  this  criterion  does  not  locate  the  reattachment 
points  at  different  axial  locations,  different  reattachment  point  pressures 
appear  to  be  physically  inconsistent.  This  raises  a  philosophical  question  as 
to  whether  or  not  any  single-stream  reattachment  criterion  can  be  rigorously 
applied  to  the  two-stream  reattachment  problem. 

Base  pressures  calculated  for  wide  ranges  of  model  parameters  were  com¬ 
pared  to  experimental  data.  These  comparisons,  in  general,  showed  the  calcu¬ 
lated  values  to  be  greater  than  the  experimental  values.  The  size  of  this 
deviation  does  not  vary  greatly  for  the  majority  of  the  comparisons  which  were 
made.  It  seems  that  a  relatively  modest  increase  of  the  ONERA  critical  angle 
would  bring  the  calculated  data  into  good  agreement  with  a  majority  of  the 
experimental  data  which  has  been  used.  No  such  modification  to  the  criterion 
has  been  attempted  since  it  could  only  be  done  on  an  empirical,  data  fitting 
basis. 
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1.0  INTRODUCTION 


The  concept  and  development  of  the  so-called  component  model  power-on 
base  pressure  analysis  was  pioneered  by  Dr.  H.  H.  Korst  and  his  coworkers  at 
the  University  of  Illinois.  Their  early  work  culminated  in  Reference  1,  which 
documents  the  component  analysis  essentially  as  It  Is  applied  today.  One  of 
the  attractive  aspects  of  their  method  Is  that  each  component  or  element  of 
the  total  number  of  components  which  comprise  the  solution  may  be  of  any 
selected  degree  of  sophistication.  Thus,  any  one  of  the  components,  may.  In 
general,  be  Improved  or  replaced  by  a  more  sophisticated  component  without 
revising  or  Improving  the  entire  program.  In  the  formulations  [1],  Dr.  Korst 
assumed  that  the  free  shear  layers  originated  at  the  separation  point  on  the 
body  or  nozzle  and  that  a  similar  velocity  profile  existed  at  every  point 
along  the  length  of  the  shear  layer,  l.e.,  the  upstream  boundary  layer  thick¬ 
ness  was  neglected.  He  also  assumed  that  the  recompression  process  across 
the  shocks  at  the  wake  closure  was  Isentroplc  to  the  stagnation  pressure  of 
the  dividing  streamline. 

Addy  [2]  developed  a  computer  program  to  Implement  Korsts's  analysis  Into 
a  working  engineering  tool  for  the  prediction  of  missile  power-on  base  pres¬ 
sure  at  supersonic  speeds.  This  program  utilizes  the  method  of  characteris¬ 
tics  to  calculate  the  Invlscld  flow  fields  and  the  numerical  solution  tech¬ 
niques  that  are  coded  Into  the  program  are  still  state-of-the-art.  However, 
the  upstream  boundary  layers  were  not  accounted  for  and  the  Isentroplc  recom- 
presslon  criteria  of  Korst  was  used  for  the  wake  closure  condition.  A  very 
limited  comparison  of  calculated  and  experimental  base  pressures  presented 
In  Reference  2  showed  reasonable  agreement  but  the  calculated  base  pressures 
were  generally  larger  than  the  experimental  data.  In  Reference  3,  Addy  devel¬ 
oped  an  empirical  modifying  factor  for  the  recompression  pressure  ratio  which 
was  a  function  of  the  nozzle  radius  to  body  radius  ratio.  The  factor  was 
developed  by  correlating  calculated  and  experimental  base  pressure  data  with 
different  values  of  this  ratio.  It  was  shown  that  this  geometric  parameter 
did  a  reasonable  overall  Job  of  correlating  the  data.  However,  where  a  suf¬ 
ficiently  large  data  base  exists  for  a  constant  geometric  ratio,  significant 
factor  variation  was  required  to  correlate  all  of  the  calculations.  Addy 
further  extended  the  utility  of  the  computer  program  by  modifying  the  external 
Invlscld  flow  field  calculation  to  account  for  boattailed  and  flared  after¬ 
bodies.  This  work  was  published  In  Reference  4  and  greatly  Increased  the 
utility  of  the  program  from  the  standpoint  of  the  number  of  configurations  to 
which  It  could  be  applied.  This  final  version  of  Addy's  program  still  did  not 
consider  the  upstream  boundary  layers,  used  the  recompression  factor  [3]  for 
cylindrical  afterbodies,  and  used  a  recompression  factor  of  one  (original 
Korst  criterion)  for  flared  and  boattailed  afterbodies.  An  extremely  limited 
comparison  of  calculated  and  experimental  data  for  a  boattalled  afterbody 
showed  the  calculated  data  to  be  reasonable,  but  no  general  trend  was 
established. 

As  discussed  above,  Addy's  computational  model  was  very  good  except  that 
It  neglected  the  upstream  boundary  layers  and  It  needed  a  better  recompression 
criterion.  The  computer  program  [4]  was  modified  to  Include  the  effects  of 
both  the  Internal  and  external  upstream  boundary  layers.  This  work,  along 
with  very  lliolted  calculations  Including  the  boundary  layer  effects,  are 
reported  In  Reference  5.  The  boundary  layer  effects  were  Included  by  modify¬ 
ing  the  shear  layer  mixing  equations,  by  calculating  an  apparent  origin  shift 
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for  the  shear  layer,  and  by  retaining  the  assumption  of  a  similar  velocity 
profile.  This  modification  prepared  the  program  for  an  updated  recompression 
criterion.  Calculated  results  reported  In  this  reference  were  for  a  wind  tun¬ 
nel  test  model  with  a  nominal  turbulent  boundary  layer  thickness  of  the  exter¬ 
nal  stream.  Calculations  were  made  for  this  cylindrical  afterbody  configura¬ 
tion  with  and  trlthout  the  boundary  layer  terms  using  Addy's  recompression 
factor  and  It  was  shown  that  Including  the  boundary  layer  effect  significantly 
Increased  the  calculated  base  pressure.  This  Is  consistent  with  other  results 
reported  In  the  literature.  It  Indicates  that  empirical  correlation  recom- 
resslon  factors  derived  without  the  boundary  layer  terms  In  the  solution,  have 
an  Inherent  boundary  layer  effect  Included  In  the  recompression  factor.  Thus, 
one  should  apply  and  evaluate  recompression  criteria  with  computational  models 
which  contain  the  boundary  layer  terms  and  are  as  complete  as  possible  In 
other  respect.  The  modified  Addy  computer  program  meets  these  criteria. 

The  goals  of  the  present  work  are:  (1)  to  Include  one  of  the  currently 
available  recompression  models  In  the  computer  program,  and  (2)  to  assess  the 
adequacy  of  the  model  by  comparing  calculated  and  experimental  base  pressures 
for  a  large  data  base.  Since  Korst's  original  development,  a  number  of 
recompression  criteria  have  been  proposed  by  different  authors.  Several  of 
these  will  be  briefly  discussed  as  a  short  survey  of  the  overall  activities  In 
this  area. 

One  of  the  earliest  proposals  for  a  recompression  criterion  was  made  by 
Nash  [6],  who  simply  proposed  that  the  Korst  Isentroplc  pressure  rise  be 
multiplied  by  a  constant  to  obtain  the  recompression  pressure  rise.  On  the 
basis  of  a  very  limited  amount  of  experimental  data,  Nash  established  a  value 
of  0.33  for  this  factor  In  the  supersonic  flow  case.  Another  approach  to  the 
recompression  criterion  Is  the  Page  criterion  which  Is  presented  In  Reference 
7,  with  References  8  and  9  discussing  application  techniques.  The  Page  cri¬ 
terion  Is  based  on  a  correlation  of  experimental  data  which  model  the  ratio  of 
the  flow  turning  angle  at  the  reattachment  point  to  the  total  turning  angle  as 
a  function  of  the  discriminating  streamline  velocity  ratio.  Carrlere  and 
Slrlelx  [10,  11]  developed  an  angular  reattachment  criterion  which  correlated 
the  effective  turning  angle  of  the  reattaching  streamline  as  a  function  of  the 
Invlscld  flow  Mach  number.  In  Reference  12,  Wagner  and  White  studied  the 
effects  of  the  Initial  boundary  layers  and  some  of  the  recompression  models. 
They  found  that  over  a  considerable  range  of  free  streamline  Mach  numbers,  the 
Page  criterion  and  the  Office  National  d'Etudes  et  de  Recherches  Aerospatlales 
(ONERA)  angular  criterion  gave  very  similar  reattaching  streamline  turning 
angles. 

The  Page  model  and  the  ONERA  angular  criterion  are  of  special  Interest 
because  both  are  presented  in  the  form  of  simple  correlations  which  are  based 
upon  comprehensive  sets  of  experimental  data.  Both  techniques  relate  the 
recompression  process  to  the  angular  deflection  at  the  reattachment  point. 

Each  also  assumes  that  the  stagnation  pressure  of  the  discriminating 
streamline  Is  equal  to  the  reattachment  static  pressure.  However,  in  the 
application  of  the  Page  model,  an  Iterative  procedure  Is  necessary  to  calcu¬ 
late  the  length  of  the  constant  pressure  mixing  region  [8]  even  for  the  single 
stream  problem.  Wagner  and  White  [12]  reported  serious  convergence  problems 
when  trying  to  Implement  the  Page  model  for  the  two-stream  base  pressure 
problem.  Therefore,  It  was  decided  to  restrict  the  present  Investigation  to 
consideration  of  the  ONERA  angular  criterion  for  the  recompression  calculation. 


2.0  ONERA  RECOHPRESSION  CRITERION 
2.1  G«neral  Discussion 

The  ONERA  angular  criterion  Is  based  upon  the  correlation  of  a  body 
of  data  froB  a  set  of  systematic  experiments  [11].  These  experiments  were 
arranged  so  that  they  corresponded  essentially  to  a  vanishing  upstream  boun¬ 
dary  layer.  In  addition,  the  change  In  reattachment  angle  caused  by  bloving 
Into  the  base  is  correlated  with  Korst's  basic  theory.  Thus,  the  recompres- 
slon  turning  angle  for  the  no-bleed,  no-lnltlal  boundary  layer  Is  to  be  taken 
from  an  experimental  data  correlation;  however,  the  effect  of  bleed  or  bound¬ 
ary  layer  influence  may  be  predicted  analytically.  With  these  empirically 
established  facts  at  hand.  It  was  postulated  that  the  basic  angular  criterion 
could  be  formulated  as 

Cq)  (1) 

For  the  usual  application  of  the  criterion.  It  Is  assumed  that  the  boundary 
layers  are  thin  with  only  a  small  bleed.  In  this  case.  It  Is  permissible  to 
linearize  Equation  (1)  so  that 

T  -  ■?  (M.)  +  C,  1^  I 


where  Cq  la  a  general  bleed  coefficient  defined  by 

r) 


where  o  Is  Included  in  the  definition  of  Cq  to  make  It  compatible  with  the 
upstream  boundary  layer  effects  on  th^  mixing  layer  equations  (see  Equation  33 
of  Reference  5).  The  critical  angle  T  (Mg)  for  vanishing  Cq  Is  obtained  from 
the  ONERA  data.  This  data  was  obtained  for  an  Isoenergetlc,  two-dimensional 
flow  with  the  ratio  of  specific  heats,  y*  "  Thus,  rigorously,  the  appli¬ 

cation  of  this  criterion  should  be  restricted  to  these  flows;  however,  equa¬ 
tions  for  the  non-lsoenergetlc  case  will  be  derived.  For  application  of 
Equation  (2)  to  the  practical  calculation  of  base  pressure,  the  function 
J'  (Mg)  needs  to  be  given  In  analytical  form.  An  empirical  curve  fit  for 
?  (Mg)  was  first  given  by  Sollgnac  and  Delery  In  Reference  13  as: 

'T(Mg)  -  57.3  (0.569  -  0.5096/Ma)  (^) 

where  ? Is  In  degrees.  Later,  Vagner,  In  Reference  14,  developed  a  polynomial 
curve  fit  of  the  data  of  the  form 


F(M-)  -  57.3  (0.044  +  0.172  M*  -  0.018  Mg^) 


This  equation  also  yields  ?  in  degrees.  The  free  streamline  Mach  number  for 
the  ONERA  data  varies  from  approximately  2.0  to  4.0  and  both  of  the  empirical 
equations  fit  the  data  very  well  i£  this  Mach  range.  However,  when  the 
equations  are  used  to  extrapolate  T  beyond  the  data  range,  sizable  differences 
occur,  and  there  is  no  substantial  evidence  to  suggest  that  either  is  superior 
to  the  other.  Therefore,  solutions  for  base  pressure  tihich  utilize  free 
streamline  Mach  numbers  outside  of  the  data  rmge  should  be  treated  with  some 
caution.  A  comparison  of  the  critical  angle  ?  deduced  from  the  above 
equations  and  other  sources  is  presented  in  Figure  1  which  is  taken  from 
Figure  6  of  Reference  14.  It  is  seen  that  not  only  do  the  equations  agree 
well  with  the  ONERA  data,  but  that  the  Page  model  also  yields  critical  angles 
which  agree  well  with  the  ONERA  data. 


2.2  Isoenergetic  Equations 

With  the  relationship  for  the  critical  angle  V  established  in  a 
usable  format,  we  turn  to  the  task  of  evaluating  the  partial  derivative  in 
Equation  (2).  As  discussed  above,  Rorst's  basic  theory  accurately  predicts 
the  change  in  the  turning  angle.  Since  the  ONERA  investigations  do  not  alter 
the  concept  of  isentropic  recompression  up  to  the  reattachment  point,  the 
angular  change  of  the  dividing  streamline,  AT,  up  to  the  reattachment  point, 
must  be  equal  to  the  reattachment  angle  calculated  %d.th  the  basic  Korst 
theory.  So,  we  have 


AT-  u<M^)  -  oXMr)  -  (6) 

where,  uis  the  Prandtl-Meyer  function  and  is  the  Mach  number  of  the  adja¬ 
cent  invlscld  streamline  at  the  reattachment  point.  Since  the  effect 
changes  only  the  streamline  angle  up  to  the  point  of  reattachment  and  the 
angle  change  from  downstream  of  the  reattachment  points  is  Independent  of  C^, 
then  there  results: 

aT  .  a(AT)  _  _  aa<MR)  _  _  dto  ^  0Mr 

9Cq  aCq  ?5Cq  dMj 

where  will  be  defined  later.  Now,  the  gradient  of  the  Prandtl-Meyer  func- 
tlon  is  given  by  Reference  15,  pg  466,  Equation  (15-3). 


d  (I)  - 


Vm2-1  dM^ 

2m2  1+  m2 


(8) 


Then 


Now  the  recompression  on  the  reattaching  streamline  is  still  assumed  to  be 
isentropic  and  the  total  pressure  on  the  discriminating  streamline  is  equal  to 
the  static  pressure  at  reattachment.  Then,  using  the  usual  Isentropic  rela¬ 
tions,  there  results 


PR/Poa 


Y 

TT 


V 


(10) 


where  is  the  Mach  number  on  the  discriminating  streamline.  This  Mach 
number  can  now  be  expressed  as 


Md^ 


“a 


(11) 


and  for  the  present  assumption  of  isoenergetlc  flow,  the  relationship  between 
the  temperature  ratio  and  the  Crocco  number  is  given  by 


(12) 


(13) 


The  relationship  between  the  Crocco  number  and  the  Mach  number  is  given  by 
Equation  (51)  of  Reference  16 


(14) 


and  when  this  result  is  substituted  into  Equation  13,  there  results 


After  substituting  Equation  (IS)  Into  the  right-hand  side  of  Equation  10  and 
simplifying 

Pod  r  1  w 


1  +  M32(l-4y2) 


The  left-hand  side  of  Equation  10  Is  now  combined  with  Equation  16  to  give 


r-1  M-2  y-l 


1  +  Mr2 


,  .  Y-1  M.2 


Equation  17  can  be  solved  for  to  obtain 

Mh2  -  (1  -  ^2)m^2  (18) 

To  continue  toward  the  evaluation  of  the  derivative  In  Equation  (2),  Equation 
(7)  must  be  further  expanded  by  the  chain  rule  to  obtain 

_ (19) 

acq  dMR*  *  acq 

Differentiating  Equation  (18)  results  In 

’d  “a  (20) 


Combining  Equations  (18),  (19),  and  (20)  with  Equation  (9)  evaluated  at  the 
reattachment  point  gives 


(21) 

where  we  eelll  here  the  requlrenent  to  eveluete  Thie  can  be  done  by 

writing  the  Integral  equation  for  the  mass  flow  In  the  shear  layer  up  to  the 
discriminating  streamline  In  the  form  (Equation  C-Al  of  Reference  17) 


Now  constants,  and  therefore  Independent  of  C-,  so 

when  Equation  (22  )  is  differentiated  with  respect  to  C^,  there  Is  obtained 


stb  p(t^)  <iia  -  1 


Expanding  Equation  (23)  by  the  chain  rule  and  solving  the  resultant  equation 
for  yields 

ac„ 


p(n^  ) 


The  velocity  profile  In  the  shear  layer  la  defined  by  the  error  function,  so 


♦  -  +  erf(T,)]- 7[l  +  ^  ^  eVdn] 


which  3rlelds  upon  differentiation 


a"  VT 

and  using  1 sen tropic  relations 
Pa  ^  1  - 
P(vJ  "  I  -  cJ 


Now,  when  Equations  (26)  and  (27)  are  combined  with  Equation  24 
1  -  1  1  .-,2^ 


1  -  c. 


The  final  step  is  to  substitute  Equation  (28)  Into  Equation  (21)  use  Equation 
(18)  for  the  definition  of  the  and  make  the  approximation  Hj  ■  n  j,  the 
result  is 


_  Vm«^  -  1 _  1  -  (tj-Ca- 

aCq  (l-*j2)  ^  1  +  ^  M^2  (l-<>j2)j  1-0^2  VT 

This  completes  the  derivation  of  the  equations  required  for  the  application  of 
the  ONERA  recompression  criterion  to  the  restricted  case  for  Isoenergetic  flow. 
Values  of  nj  and  6j  are  determined  by  Korst's  theory  In  the  course  of  each  base 
pressure  iteration  as  a  function  of  Cg,  and  therefore,  can  be  assumed  to  be 
known  for  application  to  the  solution  of  the  recompression  criterion. 


1  -  *l2c«2 


2.3  Non-I scene rge tic  Equations 


The  extension  of  the  application  of  the  ONERA  recompression  cri¬ 
terion  to  non-lsoenergetlc  flows  is  discussed  by  Delery  and  Sirlelx  [18]  and 
also  by  Wagner  [14].  Since  missile  design  cases  are  non-lsoenergetlc  flow 
cases,  the  equations  for  this  situation  are  given  herein.  Delery  and  Wagner 
use  somewhat  different  approaches  for  determining  the  critical  angle,  7,  for 
this  case,  but  they  use  approximately  the  same  approach  for  determining  the 
derivative  term.  Using  the  same  approach  as  Delery,  it  is  assumed  that 

7-  7(M3,Aa,Cq)  (30) 

and  when  is  small 

7-  7(Ma.yid)  +  Cq  ^Ma,JVi)  (31) 

where 

Ai Ab+ *d(i  -  V  (32) 

'•oa  '^oa  \  '^oa  / 


b7 

The  analysis  to  determine  the  derivative,  proceeds  as  in  the  Isoenergetic 

case,  but  %d.th  greater  complication  due  to  the  temperature  effect.  Since  the 
process  is  identical,  this  development  will  be  presented  in  a  summary  manner 
with  critical  steps  related  to  the  equations  in  Section  2.1.  The  first  sign! 
flcant  non-isoenergetic  effect  occurs  in  the  temperature  equation  for  the 
discriminating  streamline  Equation  (12)  which  is 


^  M  -  <a^Ca^ 

Ta  '  1  -  Ca2 


(33) 


And  when  this  relation  is  inserted  into  Equation  (11),  the  following  defini¬ 
tion  for  the  discriminating  streamline  Mach  number 


Aid  +  Ma2) 


(34) 


When  Equation  (34)  for  substituted  into  the  right-hand  side  of  Equatio 

(10)  and  the  resultant  is  simplified  and  solved  for  Mb, 


and  when  all  the  derivatives  along  with  Equation  35  are  substituted  Into 


Thus,  It  Is  seen  that  the  effect  of  non-lsoenergetic  flow  on  the  angular  gra¬ 
dient  can  be  correctly  calculated  by  the  above  analysis.  However,  this  effect 
upon  the  critical  angle.  In  the  absence  of  experimental  diita.  Is  subject  to 
some  conjecture.  Delery  [18]  proposed  a  modification  of  ?  for  gases  «d.th  a 
ratio  of  specific  heats,  y,  different  from  1.4.  However,  he  did  not  propose  a 
temperature  ratio  modification.  Wagner  [14]  following  the  proposal  of  Sirieix 
et  al  [11],  used  a  blowing  correction  to  a  modified  version  of  Equation  (31). 
This  modified  equation  Is 

7  -  Icorst  (Ca.Tod/Toa)  +  ^^a’^od^^oa^ 

where  the  blowing  correction,  e,  proposed  by  Wagner  Is  defined  by 

E  -  0.2934  -  0.1773  +  0.04132  -  0.00311  (40) 

In  order  for  T  to  be  approximately  equal  to  the  ONERA  7  when  «  0.  Wagner 
suggests  that  these  results  are  applicable  to  the  non-lsoenergetlc  case  and 
for  values  of  y  that  are  different  from  1.4.  These  modifications  to  the  basic 
criteria  are  easily  Implemented  within  the  basic  MICOM  Base  Pressure  program 
and  therefore,  were  selected  to  be  used  in  the  program  modification  and  sub¬ 
sequent  calculations. 


2.4  Modified  ONERA  Criterion 


Wagner  and  White  [12]  and  Wagner  [14]  developed  a  modified  ONERA 
criterion  to  account  for  Interference  effects  in  the  two  stream  base  pressure 
problems  which  are  not  present  in  the  single  stream  problem.  From  a  physical 
point  of  view,  it  appears  that  the  pressures  in  both  recompression  zones 
should  become  equal  at  the  point  of  reattachment.  The  ONERA  criterion  uses 
the  slip  streamline  angle  to  establish  the  turning  angle,  ?;  whereas  the 
modified  criterion  uses  this  direction  only  as  a  first  approximation  and  then 
varies  the  direction  until  the  pressure  ratio  becomes  equal  in  the  two  shear 
layers.  This  technique  uses  an  iteration  procedure  which  varies  the  reattach¬ 
ment  angle,  7,  that  produces  a  change  in  the  mass  flow  coefficient  when 
Equation  (2)  is  solved  for  C^.  This  leads  to  a  change  in  4^,  and  thus,  produ¬ 
ces  a  corresponding  change  in  the  pressure  ratio  from  Equation  (16). 

For  the  application  of  this  modified  criterion,  it  is  necessary  to 
compute  the  pressure  ratio  at  each  step.  Therefore,  the  velocity  ratio  is 
required  and  is  computed  from  a  relationship  involving  and  Cq.  Wagner  [14] 
showed  that  while  a  linear  relationship  is  satisfactory  for  Equation  (2),  it 
would  be  necessary  to  use  a  second  order  Taylor  series  expansion  to  achieve 
comparable  accuracy  in  the  calculation  of  Therefore, 


where,  for  the  non-laoenergetic  case 
aCq  1-Ca2 


and,  the  second  derivative  can  be  expressed  as. 


- 

^  '''oa  _  1  _l_  -rf  .  2  ^  ~  1  (^3) 

c!Cq2  aCq  “  l-c^2  'J ^  1  *  ^ 

Utilizing  Equations  (42)  and  (43)  with  Equation  (41),  there  results 


h  -  »j^Ca^  _L  e~V 
1  -  0^2  ^ 


which  completes  the  equations  necessary  to  implement  Wagner's  modified  ONERA 
criterion. 
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2.5  Axlsymmecrlc  Reattachoenc 


The  original  development  of  the  ONERA  angular  criterion  was  based 
upon  data  taken  for  two-dimensional  flow  and  the  critical  reattachment  angle 
presented  In  Figure  1  Is  restricted  to  the  two-dimensional  case.  Sollgnac  and 
Delery  [13]  conducted  a  series  of  systematic  tests  to  establish  a  set  of  reat¬ 
tachment  data  for  axlsymmetrlc  flow.  Using  similarity  considerations  and  this 
data  base,  they  were  able  to  establish  a  correlation  parameter  relating  the 
axlsymmetrlc  reattachment  data  to  the  two-dimensional  critical  angle. 
Therefore,  In  order  to  apply  the  ONERA  criterion  to  the  solution  of  the  axl¬ 
symmetrlc  base  pressure  problem,  this  correction  must  be  applied  to  the  two- 
dimensional  critical  angle.  This  correlation  Is 

AY-  -  "^-D  -  f(F)  (46) 

and  ONERA  found  that 

f(F)  -  2.7  -  3.5  tan  (5.4F  -  4.4)  (47) 


did  an  excellent  job  of  correlating  the  differences 
dimensional  and  axls3nBmetrlc  experimental  data.  In 
ratio  of  the  two-dimensional  to  axlsymmetrlc  spread 
given  by 


*^2-0  .  1 
♦ax  'W 


rds 


between  the  two- 
Equatlon  47,  F  Is  the 
rate  parameters  and  Is 


(48) 


where  L  Is  the  length  of  the  shear  layer  and  r^  Is  the  radius  at  the  reattach¬ 
ment  point.  Correlation  of  the  two  sets  of  data  by  this  function  Is  presented 
In  Figure  2  which  has  been  taken  from  References  14  and  18.  The  axlsymmetrlc 
data  was  obtained  from  three  experiments  as  Illustrated  by  the  Insets  on  the 
figure.  In  the  first  case,  there  was  a  large  angle  cone  with  slender  forward 
cylinder  and  a  small  angle  cone  ahead  of  the  cylinder.  The  forward  cone  was 
larger  In  diameter  than  the  cylinder  so  the  surface  streamlines  separated  from 
the  model  at  the  base  of  the  forward  cone.  Reattachment  occurred  on  the  sur¬ 
face  of  the  large  angle  cone.  A  second  series  of  experiments  was  performed 
with  an  under  expanded  supersonic  nozzle  which  was  surrounded  by  an  external 
shroud.  In  this  case,  the  plume  boundary  expanded  out  of  the  nozzle  and  the 
reattachment  was  on  the  external  shroud.  The  third  test  to  obtain  data  used 
an  axlsymmetrlc  body  with  a  rearward  facing  step  followed  by  a  cylindrical 
afterbody.  Streamlines  separated  from  the  forward  surface  at  the  step  and 
reattachment  was  on  the  downstream  cylinder.  The  function  F  must  be  calcu¬ 
lated  for  both  flow  boundaries  In  the  axlsymmetrlc  base  pressure  problem  and 
then  the  critical  angles  are  obtained  by  applying  the  correction.  All  of 
these  calculations  are  easily  performed  within  the  framework  of  the  basic 
MICOM  Base  Pressure  program. 

3.0  COMPUTER  PROGRAM  MODIFICATIONS 


With  the  completion  of  the  present  work,  there  have  been  two  stages  of 
modifications  on  the  Addy  base  pressure  program  which  are  documented  In 
Reference  4.  First,  the  effects  of  the  upstream  boundary  layer  and  the  effec¬ 
tive  shear  layer  origin  shift  were  Included  In  the  program  and  were  documented 
in  Reference  5.  In  the  present  work,  the  ONERA  recompression  criterion  and  a 


modified  ONERA  criterion  have  been  included  to  replace  the  empirical 
recompresaion  coefficient  derived  by  Addy{3].  Most  of  the  calculations  for 
the  new  recompression  criterion  have  been  included  in  the  subroutine  TJMIX  of 
the  original  program.  The  solutions  for  these  criteria  implemented  in  this 
subroutine  have  been  adapted  from  the  Wagner  and  White  program  which  was  also 
based  upon  the  MICOM  Base  Pressure  program. 

The  program  has  been  run  for  a  large  number  of  isoenergetic  cases  and  has  pro¬ 
ven  to  be  very  reliable  when  operating  in  this  mode  for  either  the  standard  or 
modified  ONERA  criterion.  Results  from  these  calculations  will  be  discussed 
and  compared  with  some  available  experimental  data  in  Section  4.0  of  this 
report.  A  limited  number  of  non-isoenergetic  cases  were  successfully  run  with 
both  internal  and  external  gases  having  a  ratio  of  specific  heats  equal  to 
1.4.  However,  when  the  ratio  of  the  inner  to  outer  stream  temperatures  became 
large,  the  program  would  fall.  The  program  would  not  operate  for  the  non- 
isoenergetic  case  when  the  internal  gas  ratio  of  specific  heats  was  not  equal 
to  1.4.  Therefore,  it  is  apparent  that  irrespective  of  the  performance  of  the 
angular  criterion  in  the  basic  Isoenergetic  case  for  which  it  was  developed, 
additional  development  will  be  required  before  it  can  be  used  for  the  general 
non-isoenergetic  missile  design  case. 

This  program  is  available  as  a  private  user  (P)  file  in  the  Perkln-Elmer 
account  which  was  assigned  to  the  delivery  order  under  which  this  study  was 
performed.  It  is  also  available  on  a  backup  tape  for  safekeeping.  The  file 
name  for  this  problem  is  'BPRESOC  which  stands  for  Base  Pressure-ONERA 
Criterion.  Even  though  the  program  is  named  ONERA  criterion,  it  will  also 
calculate  with  the  modified  criterion  by  selecting  the  appropriate  value  of  a 
logical  variable.  The  progrm  is  set  up  to  run  in  an  interactive  mode  on  the 
Perkin-Elmer  3230  tdth  the  option  to  either  use  a  data  file  for  most  of  the 
data  or  to  input  all  of  the  data  interactively.  A  definition  of  the 
Perkin-Elmer  interactive  screen  cues  and  input  data  is  presented  in  Appendix 
A.  The  complete  FORTRAN  listing  of  the  program  is  printed  in  Appendix  B  to 
provide  a  permanent  record  of  the  current  state  of  the  program. 

4.0  PRESENTATION  OF  RESULTS 

Calculations  have  been  performed  for  a  large  number  of  isoenergetic  cases 
for  comparison  with  experimental  data  and  %rlth  Wagner  and  White's  calculations 
which  are  presented  in  Reference  12.  Some  reasonably  recent  data  such  as  that 
presented  in  Reference  19  were  found  and  most  of  the  data  sources  that  Addy 
used  for  his  empirical  correlation  [3]  will  be  presented.  The  various  com¬ 
parisons  which  are  given  represent  a  rather  wide  cross  section  of  geometric 
configurations,  nozzle  configurations,  and  free  stream  Mach  numbers. 

Generally,  the  data  sample  seems  sufficiently  large  to  support  any  general 
conclusions  that  might  be  reached. 

4.1  Comparison  with  Wagner  and  White's  Results 

Wagner  and  White  [12]  compared  results  from  their  calculations  with 
experimental  data  for  several  cases.  Calculations  were  performed  with  the 
present  progrm  to  compare  with  the  Wagner  and  White  calculations  and  the 
appropriate  experimental  data  is  included  to  complete  the  theory  to  data 
comparison.  Figure  3  presents  the  comparison  of  the  two  sets  of  calculations 
and  experimental  data  for  a  cylindrical  afterbody  configuration  of  Agrell  and 


White  [19].  The  Wagner  and  White  calculations  were  extracted  from  Figure  5  of 
that  paper.  Wagner  and  White  Investigated  the  effects  of  origin  shift,  spread 
rate  parameter,  and  recompression  criterion.  From  these  and  various  other 
results.  It  was  concluded  that  the  ONERA  origin  shift  [11],  and  the  Korst  and 
Trlpp^  [20]  spread  rate  parameter  definition  would  be  used  exclusively  In  the 
current  Investigation.  The  two  sets  of  calculations  for  the  standard  and 
modified  ONERA  criteria  agree  extremely  well,  as  Indeed  they  should.  The  dif¬ 
ference  between  the  modified  and  standard  criteria  Is  small  at  low  pressure 
ratios  and  It  steadily  Increases  as  Che  pressure  ratio  Increases.  Agreement 
between  Che  modified  criterion  calculations  and  Che  experimental  data  Is  very 
good  at  Che  higher  pressure  ratios  and  not  quite  so  good  at  Che  lower  pressure 
ratios.  Also  Included  In  this  figure  as  a  reference  for  the  recompression 
criterion  effect,  are  calculations  using  Che  original  Korst  Isentroplc  recom- 
presslon  with  the  boundary  layer  terms  and  origin  shift  Included.  Figure  4 
presents  Che  mathematical  variation  of  Che  base  pressure  as  a  function  of  the 
upstream  boundary  layer  momentum  thickness  for  the  cylindrical  afterbody  con¬ 
figuration  of  Figure  3  for  a  pressure  ratio  equal  to  3.0.  Current  calcula¬ 
tions  for  both  the  standard  and  modified  criteria  are  compared  to  Che  same 
types  of  calculations  by  Wagner  and  White.  Again,  there  Is  excellent  agree¬ 
ment  between  Che  two  sets  of  calculations  while  both  types  of  calculations 
somewhat  overpredlct  the  experimental  data  point  extracted  from  Figure  3. 

One  Important  observation  from  this  figure  Is  the  relatively  small  gradient 
of  Che  base  pressure  with  Increasing  momentum  thickness  beyond  a  nominal 
thickness  value.  Therefore,  for  experimental  reports  which  do  not  include 
boundary  layer  data.  It  appears  to  be  reasonable  to  estimate  the  momentum 
thickness  using  engineering  methods. 

A  comparison  of  the  same  types  of  calculations  for  this  same  configuration 
as  shown  In  Figure  3  Is  shown  In  Figure  S,  but  the  free  stream  Mach  number  has 
been  Increased  to  3.27.  There  Is,  again,  excellent  agreement  of  the  results 
from  Che  two  programs  for  both  the  standard  and  modified  ONERA  criteria  and 
also  for  calculations  using  Che  Addy  recompression  coefficient.  Calculated 
results  using  Che  two  ONERA  criteria  agree  well  with  the  experimental  data  at 
low  pressure  ratios  but  a  *  ~ery  much  greater  at  the  higher  pressure  ratios. 
This  Is  inverse  to  Che  trend  of  Figure  3.  Figure  6  continues  this  series  of 
comparisons  of  calculations  and  experimental  data.  For  the  first  time,  there 
are  differences  between  Che  current  calculations  and  those  of  Wagner  and 
White.  The  configuration  for  Che  data  of  Figure  6  Is  a  six-degree  conical 
boattall  afterbody.  The  other  analytical  comparison  In  the  figure  utilizes  an 
empirical  recompression  coefficient  for  boattailed  bodies  which  was  developed 
by  Dr.  White  through  a  correlation  of  the  experimental  data  [19].  Since  there 
Is  extremely  good  agreement  between  these  two  solutions,  a  strong  Indication 
Chat  both  programs  are  performing  the  boattall  flow  field  solution  correctly, 
the  difference  between  Che  solutions  utilizing  Che  two  ONERA  criteria  Is 
puzzling.  A  final  comparison  of  current  calculations  between  the  Wagner  and 
White  calculations  and  the  experimental  data  [21]  Is  presented  In  Figures  7 
through  10.  Figure  7  presents  Che  base  drag  coefficients  for  all  three  con¬ 
figurations  for  which  data  Is  presented  In  Reference  21.  Since  Che  base  drag 
coefficient  Is  somewhat  Insensitive  Co  Che  base  pressure  ratio.  It  was  decided 
to  recast  Che  data  of  Figure  7  Into  the  base  pressure  ratio  and  this  Is  pre¬ 
sented  In  Figures  8  through  10  for  Che  three  different  nozzles.  There  are 
small  differences  between  Che  current  calculations  and  the  Wagner  and  White 
calculations.  However,  Che  differences  between  Che  calculations  are  greater 
for  the  largest  nozzle.  Two  comparisons  with  calculations  published  by 


Dr.  Wagner  [14]  are  presented  In  Figures  11  and  12.  Figure  11  compares 
calculated  results  from  the  two  programs  with  additional  data  from  Reference 
21.  The  model  configurations  for  the  data  presented  In  this  figure  a.'e  a 
cylindrical  afterbody  with  a  sonic  nozzle  and  a  Mach  2.0  nozzle  with  uniform 
parallel  flow  In  the  exit  plane.  Current  calculations  agree  well  with 
Wagner's  calculations  and  both  sets  predict  base  pressures  somewhat  higher 
than  the  experimental  data.  Figure  12  presents  a  comparison  of  results  from 
the  two  programs  for  two  configurations  for  which  test  data  Is  available  In 
Reference  22.  The  current  calculations  yield  greater  values  of  base  pressure 
than  do  Wagner's  calculations  for  both  the  ONERA  criterion  and  the  modified 
criterion.  This  disparity  Is  disturbing  because  It  seriously  affects  the  dif¬ 
ference  between  the  calculations  and  the  experimental  data.  However,  a  con¬ 
certed  analysis  of  the  current  calculations  failed  to  disclose  any  reason  for 
the  differences  between  the  two  sets  of  calculations. 

Overall,  the  comparisons  of  the  current  calculations  with  those  of  Wagner 
and  White  show  good  to  excellent  agreement.  This  Instills  confidence  that  the 
analyses  and  program  modifications  have  been  executed  correctly.  On  the  other 
hand,  comparisons  of  calculated  and  experimental  results  show  some  cases  of 
good  agreement  and  some  of  poor  agreement.  However,  more  comparisons  of 
calculated  and  experimental  data  need  to  be  made  before  drawing  any  conclu¬ 
sions  concerning  the  overall  prediction  capability  of  the  computer  program. 

4.2  Comparisons  with  Experimental  Data 

The  overall  effectiveness  of  the  base  pressure  computer  program  In 
predicting  base  pressure  %d.ll  be  Judged  by  a  broad  based  comparison  of  calcu¬ 
lated  and  experimental  data.  Addy  [3]  collected  a  wide  ranging  data  base  to 
develop  the  correlation  for  his  recompression  coefficient.  Since  Addy's 
report  was  published,  Agrell  and  White  published  Reference  19  which  contains 
considerable  additional  data.  Experimental  data  for  the  comparisons  to  be 
discussed  In  this  section  will  be  obtained  from  these  sources.  Figure  13 
completes  the  comparison  of  the  current  calculations  with  the  data  of  Reid  and 
Hastings  (see  Reference  21).  Figures  7  through  11  for  the  other  comparisons. 
Agrell  and  White  [19]  published  a  rather  comprehensive  set  of  data  for 
cylindrical  and  boattalled  models  at  two  free  stream  Mach  numbers  and  two  dif¬ 
ferent  nozzles.  Some  of  the  data  Indicates  that  flow  separation  occurs  for¬ 
ward  of  the  base  for  the  larger  boattall  angles  and  pressure  ratios.  This 
data  was  not  used  for  comparison  with  calculations  since  the  current  program 
does  not  have  flow  separation  prediction  capability.  Some  of  the  Agrell  and 
White  data  Is  presented  In  Figures  3,  5,  and  6.  The  rest  of  the  Agrell  and 
White  data  for  attached  flow  on  the  afterbody  are  compared  with  current  calcu¬ 
lations  for  both  the  standard  and  modified  ONERA  criteria  In  Figures  14 
through  20.  In  general,  the  calculated  base  pressures  are  larger  than  the 
experimental  values.  The  amount  of  the  overpredlctlon  depends  on  the  nozzle 
characteristics,  boattall  angle,  and  free  stream  Mach  number.  At  M*  -  2.01, 
the  basic  trends  of  both  calculated  and  experimental  data  are  quite  similar 
over  the  pressure  ratio  range.  However,  at  M<»  -  3.27,  the  variation  of  the 
experimental  data  over  the  pressure  ratio  range  Is  much  smaller  than  the 
variation  of  the  calculations. 
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Experimental  data  of  Bromm  and  O'Donnel  [22]  are  compared  with  calculated 
values  of  base  pressure  and  presented  In  Figure  13  and  In  Figures  21  through 
25.  Bromm  and  O'Donnell  obtained  data  at  three  free  stream  Mach  numbers,  two 
nozzle  sizes  and  three  nozzle  wall  angles.  Once  again,  the  calculated  data 
are.  In  general,  greater  than  the  experimental  data;  however,  the  amount  of 
the  overpredlctlon  varies  significantly  from  case  to  case.  It  Is  suspected 
that  some  of  this  variation  Is  due  to  experimental  uncertainties,  especially 
for  the  ten  degree  nozzle  exit  angle.  Qualitative  trends  are  well  predicted 
by  the  calculated  results  and  the  experimental  trends  remain  similar  for  all 
three  free  stream  Mach  numbers.  Baughman  and  Kochendorfer  [23]  published  data 
for  bodies  with  five  different  boattall  angles  with  the  base  to  nozzle  radius 
ratio  determined  by  the  boattall  length.  The  models  were  each  tested  with 
constant  size  nozzles  that  produced  different  exit  plane  Mach  numbers.  These 
nodel  and  nozzle  combinations  were  tested  at  free  stream  Mach  numbers  of  1.91 
and  3.12.  Comparisons  of  calculated  and  experimental  base  pressures  for  these 
conditions  are  presented  In  Figures  26  through  38.  Again,  the  calculated  data, 
In  general.  Is  greater  than  the  experimental  data.  Qualitative  trends  are 
well  predicted  for  the  Mach  1.91  data.  However,  the  trend  of  the  variation 
with  pressure  ratio  was  not  well  predicted  for  the  Mach  3.12  data  and  this 
same  result  was  observed  in  the  comparisons  with  the  Agrell  and  White  data. 

The  next  set  of  comparisons  presented  Is  with  the  data  of  Henderson  [24]. 
These  comparisons  are  limited  to  Henderson's  cylindrical  afterbody  configura¬ 
tion  and  they  are  presented  In  Figures  39  through  44.  All  of  the  calculated 
base  pressures  are  greater  than  the  corresponding  experimental  data.  Thus, 
the  trend  that  has  been  established  through  all  the  previous  comparisons  con¬ 
tinues  without  exception.  Comparisons  of  calculated  base  pressures  with  the 
experimental  data  of  Reference  25  are  presented  In  Figures  45  through  53.  The 
configuration  for  the  data  and  calculations  of  Figures  45  to  47  is  a  cylindri¬ 
cal  afterbody  with  a  small  nozzle  and  the  mismatch  of  data  with  calculations 
Is  probably  the  largest  yet  observed.  However,  this  follows  the  trend  of  the 
Reid  and  Hastings  data  presented  earlier*  The  remainder  of  the  data  Is  for 
nine-degree  boattall  models  and  significant  mismatches  have  been  observed  for 
similar  configurations  In  earlier  comparisons.  Figures  54  and  55  present  com¬ 
parisons  of  calculated  base  pressures  with  experimental  data  of  Cortrlght  and 
Schroeder  [26].  As  In  all  previous  comparisons,  the  calculated  values  are 
greater  than  the  experimental  values.  Data  from  Harries  [27]  are  compared  to 
the  calculated  base  pressures  In  Figures  56  through  60.  The  Isoenergetlc  data 
of  Reid  [28]  Is  compared  to  the  equivalent  calculated  values  In  Figures  61  and 
62.  Comparisons  of  calculated  base  drag  coefficients  with  experimental  data 
[29]  Is  presented  In  Figures  63  through  65.  The  data  presented  Is  for  the  no¬ 
bleed  test  cases.  Since  the  predicted  base  drag  coefficient  Is  lower  than  the 
test  data,  the  predicted  base  pressure  continues  to  be  greater  than  experimen¬ 
tal  values.  Also,  at  a  Mach  number  of  3.5,  the  trend  over  the  pressure  ratio 
range  Is  not  well  predicted.  This  continues  the  pattern  noted  In  comparisons 
which  were  previously  discussed.  Figures  66  through  69  present  comparisons  of 
calculated  base  pressures  with  single  jet  experimental  data  from  Reference  30. 
This  is  the  only  set  of  data  which  generally  agrees  with  the  calculations. 

The  geometric  configuration  Is  quite  similar  to  configurations  for  which  com¬ 
parisons  have  been  discussed  previously  where  calculations  and  data  did  not 
agree.  Therefore,  It  Is  felt  that  this  set  of  data  Is  not  consistent  with  the 
rest  of  the  data  base.  The  final  comparisons  to  be  presented  are  shown  in 
Figures  70  through  73  where  the  data  Is  taken  from  Reference  31.  These 
figures  show  the  usual  trend  of  the  calculated  values  being  greater  than  the 
experimental  data. 


COMPUTER  PLOT  LABEL  DEFINITIONS 


Ordinates 

PB/PE  Ratio  of  base  pressure  to  free  stream  static  pressure 

CD-B  Base  drag  coefficient 

CP-B  Base  pressure  coefficient 

Abscissae 

PI/PE  Ratio  of  nozzle  exit  plane  static  pressure  to  free  stream 

static  pressure 

THETA/R  Ratio  of  external  boundary  layer  momentum  thickness  to  maximum 

body  radius 

POI/PE  Ratio  of  nozzle  stagnation  pressure  to  free  stream  static 

pressure 

Parameters 

BETAe  Body  (external  flow)  angle  at  base 

BETA!  Nozzle  wall  angle  at  exit  plane 

Me  Free  stream  Mach  number 

Mi  Nozzle  exit  plane  Mach  number 

Ri/R  Ratio  of  nozzle  exit  radius  to  maximum  body  radius 

Ratio  of  external  boundary  layer  momentum  thickness  to  maximum 
body  radius 
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f  calculated  base  pressures  with  experimental 
ference  28,  (figure  13),  Isoenergetic  data. 
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son  of  calculated  base  drag  coefficients  with  experimental  data 
ference  29,  (figure  4),  no  bleed,  M»  =  2.5. 
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mparison  of  calculated  base  pressures  with  experimental 
ta  from  reference  30,  (figure  10a),  single  jet,  H,,  =  2. 
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GLOSSARY 


Definition 

General  bleed  coefficient 

Crocco  number  (U/U  ) 
max 

Two-dimensional  to  axlsymmetrlc  spread  rate  ratio 

Integral  defined  by  Equation  C-38  of  Reference  17 

Integral  defined  by  Equation  C-39  of  Reference  17 

Total  length  of  mixing  layer 

Mass  bleed 

Mach  number 

Pressure 

Temperature 

x-component  of  the  velocity  within  the  shear  layer  or 
boundary  layer 

Mixing  layer  origin  shift  due  to  upstream  boundary  layer 
Intrinsic  coordinates  in  the  two-dimensional  mixing  region 
Modifying  bleed  coefficient 
Ratio  of  the  specific  heats 

Stagnation  temperature  ratio  (see  Equation  32) 

Dimensionless  coordinates  in  the  mixing  region  [oy(x^  +  x) ] 

Critical  reattachment  angle  for  the  no  bleed  case 

Critical  reattachfflint  angle 

Density 

Spread  rate  or  mixing  parameter 
Boundary  layer  momentum  thickness 
PRANDTL-Meyer  function 


GLOSSARY 

Definition 

Axlsymmetrlc 
Adjacent  Invlsdd  flow 
Adjacent  quiescent  region 
Base  Region 

Discriminating  streamline 
Stagnation  conditions 
Reattachment  conditions 
Upper  edge  of  shear  layer 
Lower  edge  of  shear  layer 
Two-Dimensional 


The  current  version  of  the  program  Is  configured  to  either  (1)  read  part  of 
the  data  In  an  Interactive  mode  from  a  terminal  and  part  of  the  data  from  a 
Namelist  data  file  or  (2)  to  read  all  of  the  data  In  an  Interactive  mode 
from  a  terminal.  Data  file  Input  was  found  to  be  an  efficient  way  to  run 
the  program  when  multiple  cases  were  to  be  run  with  perhaps  only  one  (1) 
Input  variable  changed.  For  single  case  runs,  the  Interactive  mode  Is  the 
most  efficient  way  to  run  the  program.  The  Namelist  Input  file  Is  name 
DATA,  and  the  definition  statement  Is 

NAMELIST/DATA/NPRINT,  NSHAPE,  XIE,  RIE,  X2E,  R2E,  BETD2E,  EMNE,  XII, 
Rll,  BETDll,  GCI,  GAMMAI,  EMNII,  TROEI,  KPRESR 

Definitions  of  these  variables,  with  the  exception  of  NPRINT  Is  given  with 
the  appropriate  cue  below.  NPRINT  Is  an  output  control  parameter  defined 
as  follows: 

NPRINT  •  -1,  Input  data  and  base  pressure  solution  printed 
NPRINT  -  0,  Input  data,  iterations  and  solution  printed 

NPRINT  ■  +1,  Input  data.  Iterations  with  constant  pressure  boundary 
data,  and  solution  printed 
NOTE:  In  the  total  Interactive  mode,  NPRINT  ■  -1. 

Cues  which  appear  In  the  program  are  listed  and  the  Input  quantity  defined 
In  the  order  In  which  they  appear  In  the  program.  Where  data  Is  optional, 
depending  on  other  Input  data,  notes  are  Included  explaining  the  require¬ 
ment  or  option. 

CUE:  ENTER  INPUT  OPTION  CHOICE,  1  -  FILE  2  -  TERMINAL 
READ  (5,  72)  INOPT 
FORMAT  (II) 

INOPT  -  INPUT  DEVICE  SELECTOR 


CUE:  ENTER  NRCMP  RECOMPRESSION  CHOICE 

0  -  EMPIRICAL  RECOMP  COEFF  AFTER  ADDY  OR  WHITE 

1  -  KORST  APPROXIMATION  OF  PAGE  CRITERION 

2  -  CORRECT  PAGE  CRITERIA  CALCULATED  ITERATIVELY 

3  -  ONERA  ANGULAR  REATTACHMENT  CRITERION 

4  -  ONERA  CRITERIA  MODIFIED  FOR  2-STREAM  REATTACK 

5  -  AS  4,  SECOND  ORDER  CORRECTION  FOR  PHID 

6  =-  AS  4,  NONLINEAR  TREATMENT  (RECOMP  -  l.NEC.) 

READ  (5,  13)  NRCMP 

FORMAT  (II) 

NOTE:  Only  options  0,  3,  4,  and  5  are  operational.  It  Is  suggested 

that  If  the  Addy  or  White  empirical  recompression  coefficient  Is 
used  that  the  boundary  layer  momentum  thicknesses  be  set  equal 
to  zero. 

The  following  cue  appears  only  If  NRCMP  •  0  Is  chosen,  then 

CUE:  ENTER  RECOMP 

READ  (5,  6)  RECOMP 
FORMAT  (F6.4) 

RECOMP  -  SELECTOR  FOR  RECOMPRESSION  COEFFICIENT 

NOTE:  If  RECOMP  ■  0. ,  a  value  Is  calculated  for  RECOMP  from  Addy's  or 
White's  empirical  equations.  If  a  finite  value  of  RECOMP  Is 
Input,  it  is  used  for  the  calculation. 

CUE:  ENTER  ALPHANUMERIC  HEADING  -  20A4 
READ  (5,  10)  (A(I), 1-1,20) 

FORMAT  (20A4) 

A(I)  -  ANY  ALPHANUMERIC  PROBLEM  DESCRIPTION 


CUE:  ENTER  NSHIFT 


NOTE: 


CUE: 


MESSAGE; 


CUE: 


READ  (5,  9)  NSHIFT 
FORMAT  (II) 

NSHIFT  -  0,  ONERA  ORIGIN  SHIFT  (RECOMMENDED) 

NSHIFT  -  1,  SHEAR  LAYER  AND  BOUNDARY  LAYER  MOMENTUM  THICKNESSES 
MATCHED 

If  a  Namelist  data  file  is  to  be  used,  the  next  cue  to  appear  is 
the  cue  for  the  number  of  pressure  ratio  cases  to  be  run.  For 
the  interactive  mode,  the  following  cue  appears. 


AFTERBODY  SHAPE  PARAMETERS 
0  -  CYLINDRICAL  AFTERBODY 

1  -  OGIVE  BOATTAIL 

2  -  PARABOLIC  BOATTAIL 

3  -  CONICAL  BOATTAIL  OR  FLARE 
ENTER  AFTERBODY  SHAPE  PARAMETER  -  II 
READ  (5,  13)  NSHAPE 

FORMAT  (II) 

NSHAPE  -  PARAMETER  DEFINING  THE  AFTERBODY  SHAPE 

BODY  AND  NOZZLE  DIMENSIONS  ARE  RELATIVE 
THEY  CAN  BE  INCHES,  FEET,  OR  CALIBERS 
"X"  DIMENSIONS  ARE  POSITIVE  AFT 

ENTER  "X"  AT  BODY  BASE 
READ(5,  11)  XIE 
FORMAT  (4F10.4) 

XIE  •  LONGITUDINAL  COORDINATE  OF  POINT  WHERE  SEPARATION  OF  THE 
EXTERNAL  STEAM  OCCURS 


w  v’f 


CUE;  ENTER  RADIUS  AT  BODY  BASE 
READ(5,  11)  RIE 
FORMAT  (4F10.4) 

RIE  -  RADIAL  COORDINATE  OF  POINT  WHERE  SEPARATION  OF  THE 
EXTERNAL  STREAM  OCCURS 


The  following  three  cues  are  dependent  upon  the  value  of  NSHAPE.  If 
NSHAPE  ■  1,  2,  or  3,  these  cues  appear.  If  NSHAPE  ■  0,  they  do  not  appear 
and  the  cue  for  the  free  stream  Mach  number  appears. 


CUE;  ENTER  ”X"  AT  START  OF  BOATTAIL 
READ(5,  11)  X2E 
FORMAT  (4F10.4) 

X2E  -  INITIAL  LONGITUDINAL  COORDINATE  OF  THE  BOATTAIL 

CUE;  ENTER  RADIUS  AT  START  OF  BOATTAIL 
READ(5,  11)  R2E 
FORMAT  (4F10.4) 

R2E  -  INITIAL  RADIAL  COORDINATE  OF  THE  BOATTAIL 

CUE;  ENTER  SLOPE  AT  START  OF  BOATTAIL  -  DEG 
READ  (5,  11)  BETD2E 
FORMAT  (4F10.4) 

BETD2E  -  INITIAL  BOATTAIL  ANGLE  (IN  DEGREES)  AT  (X2E,  R2E) 
COUNTER-CLOCKWISE  FROM  X-AXIS  IS  POSITIVE. 

CUE;  ENTER  FREE  STREAM  MACH  NUMBER 


READ  (5,  11)  EMNE 
FORMAT  (4F10.4) 

EMNE  -  EXTERNAL  FREE  STREAM  MACH  NUMBER 


ENTER  "X"  AT  END  OF  NOZZLE 
READ  (5,  11)  XII 
FORMAT  (4F10.4) 

XII  -  LONGITUDINAL  COORDINATE  OF  POINT  WHERE  SEPARATION  OF  THE 
INTERNAL  STREAM  OCCURS 

ENTER  NOZZLE  EXIT  RADIUS 
READ  (5,  11)  Rll 
FORMAT  (4F10.4) 

*  RADIAL  COORDINATE  OF  POINT  WHERE  SEPARATION  OF  THE 
INTERNAL  STREAM  OCCURS 

ENTER  NOZZLE  EXIT  ANGLE  >  DEC 
READ  (5,  11)  BETDII 
FORMAT  (4F10.4) 

BETDII  -  FLOW  ANGLE  (IN  DEGREES)  AT  (XII,  Rll)  COUNTER-CLOCKWISE 
IS  POSITIVE 

ENTER  NOZZLE  GAS  CONSTANT  -  LBF/LBM  R 
53.34  FOR  AIR 
READ  (5,  11)  GCI 
FORMAT  (4F10.4) 

GCT  -  GAS  CONSTANT  FOR  THE  INTERNAL  STREAM 

ENTER  NOZZLE  GAMMA  -  1.4  FOR  AIR 
READ  (5,  11)  GAMMAI 
FORMAT  (4F10.4) 

GAMMAI  -  RATIO  OF  SPECIFIC  HEATS  FOR  THE  INTERNAL  GAS 

ENTER  NOZZLE  EXIT  MACH  NUMBER 
READ  (5,  ig  EMNII 
FORMAT  (4F10.4) 

EMNII  -  MACH  NUMBER  AT  (JClI,  Rll) 


CUE:  ENTER  EXTERNAL  TO  INTERNAL  STREAM  STAGNATION  TEMPERATURE  RATIO 
-  TOE/TOI 
READ  (5,  11)  TROEl 
FORMAT  (4F10.4) 

TROEl  -  STAGNATION  TEMPERATURE  RATIO  OF  STREAMS,  TOE/TOI 

CUE:  ENTER  NUMBER  OF  CASES  -  II 
READ  (5,  13)  NCASE 
FORMAT  (II) 

NCASE  -  NUMBER  OF  PRESSURE  RATIOS,  PlI/PE  OR  POI/PE,  FOR  WHICH 
BASE  PRESSURE  CALCULATIONS  ARE  TO  BE  MADE  FOR  A  GIVEN 
SET  OF  CONDITIONS  AND  GEOMETRY 

CUE:  ENTER  TYPE  OF  PRESSURE  RATIO  INPUT 

0  FOR  INTERNAL  STATIC/EXTERNAL  STATIC  PRESSURE 
1  FOR  INTERNAL  STAGNATION/EXT  STATIC  PRESSURE 
READ  (5,  13)  KPRESR 
FORMAT  (II) 

KPRESR  -  0,  PRIIE  (PlI/PE)  IS  INPUT,  AND  PROIE  IS  CALCULATED. 

1,  PROIE  (POI/PE)  IS  INPUT,  AND  PRIIE  IS  CALCULATED. 

The  follwlng  cues  depend  upon  the  value  of  KPRESR.  If  KPRESR  ■  0,  then 

CUE:  ENTER  PJ/PP 

READ  (5,11)  PRATIO 
FORMAT  (4P10.4) 

PRATIO  -  PRESSURE  RATIO  (Pll/PE  -  PRIIE)  FOR  CASE  1. 

If  KPRESR  -  1,  then 

CUE:  ENTER  POJ/PP' 

READ  (5,  11)  PRATIO 
FORMAT  (4F10.4) 

PRATIO  -  PRESSURE  RATIO  (POI/PE  -  PROIE)  FOR  CASE  1. 


CUE:  ENTER  BOUNDARY  LAYER  MOMENTUM  THICKNESS  AT  BODY  BASE 
READ  (5,  11)  BLMTE 
FORMAT  (4F10.4) 

BLMTE  -  EXTERNAL  FLOW  BOUNDARY  LAYER  MOMENTUM  THICKNESS 
AT  SEPARATION  POINT  -  SAME  UNITS  AS  BODY 
DIMENSIONS 

CUE;  ENTER  BOUNDARY  LAYER  MOMENTUM  THICKNESS  AT  NOZZLE  EXIT 
READ  (5,  11)  BLMTI 

BLMTI  =  INTERNAL  FLOW  BOUNDARY  LAYER  MOMENTUM  THICKNESS 

AT  NOZZLE  EXIT  PLANE  -  SAME  UNIT  AS  BODY  DIMENSIONS 
The  cues  PJ/PR  or  POJ/PJ  along  with  the  cues  for  BLMTE  and  BLMTI  reappear 
after  the  solution  for  each  pressure  ratio  Is  obtained  until  the  NCASE 
solutions  have  been  completed.  After  the  solution  for  the  last  pressure 
ratio  has  been  completed,  the  following  cue  appears. 

CUE:  DO  YOU  WISH  TO  MAKE  ANOTHER  RUN?  YES  -  1,  NO  -  0 
READ  (5,  *)  IRUN 
FREE  FORMAT 

IRUN  -  0,  NORMAL  PROGRAM  EXIT 
1,  RETURN  TO  FIRST  CUE 

Cues  are  typed  above  exactly  as  they  appear  on  the  screen  and  In  the  order 
In  which  they  appear.  For  the  Input  data  which  are  real  numbers,  a  format 
cue  In  the  form  of  XXXXX.XXXX  also  comes  In  the  screen.  Mnemonics  for  the 
Input  data  are  the  ones  used  In  Reference  4  and  the  definitions  are  basically 
the  same  as  those  presented  In  Reference  4. 


A-7/(A-8  blank) 


rionn  o  n  r*nn  n 


is*t*pp(  input, piiii'ur, nH.cn, T#^•^s«J^PUT,  I A('eo»«i<Tmi , 

1  rAPCTzpIlMCHt 

TaU-STRLAH  AAlSYNHfTPtC  HA5K 
PRrSSURF  PR  U  gran  ,  TSAHPP  -  2. 

AFTF.  RBOIII  OPTIUNPL  HtFUKF 
PnOGRAN  H*3  ROTH  InTIPNAL  AM>  STPtAN 

nOIINORT  LATFR  kFftCIS  IMCLUDFP  IN  THF  FURnAT  OF 
EFtClivr  NASS  HLFFI)  and  a  Nixmc  LATER  ORIGIN  SHIFT. 

****«MtSTEP  RHJUIRFS  INOuT,  UUTIP,  HOT2“,  ACPHS,  CROSS,  TJMIX, 

JTFR.  THF  VARIOUS  SUHHUUTINFS  CALL  UTMFRS. 


DlHfNSlUN  PNP(loO,5,2),  CMARl(5,J0>,  CHA RF ( S , 30 ) ,  PI(5I,  P2(5), 

I  PI(S),  A(20|,  1>ATA(I0,2),  •>PTI(5,30),  rtPTe(S,3U) 

CONNUN  pNB,  CHAHl,  CHARI,  PI,  P2,  PI 
CUNNOH  /ANC.LFS/THFTAl  ,  THFT  A2 , 1  HFT  AS  ,F  I  ,  F2 
CONNUN  /LRFVP/  PHt(lSU) 

COMHUN  /OATAtO/  GCl,r:ANNAl,F.NSll,Xtr,RII,HETAlf, 

1  CCF,GAhHAF..FNSlF  ,  X I F  ,  R 1 C  ,  ULT  A I E ,  PROl  OE, 

2  IRuEl  ,PRllC,PtCOap,A,ENNU,PRl<)l,EMNIF,PRiniE, 

J  NPK1NT,NCASI  ,NCASF,RLnRU,CNGRU,NF.  ,CNNF.,PRF.UE, 

4  NPUNCH,PR0E0t,PRIItC,P0tF01,NSHAPt,NPTSE,PRl  I  IE, 

i  NDEFLT,bLNTE,SLMri 

CUMMUN  /Cat/COHT,xriX,KI IX, NSCP,SLPRtS,XTEST,ECb,SFPR, DELTA 
CONNOR  /Ca2/X2C,R2C 

CONNUN  /CaJ/SIGNFl,StGNl2,NRCNP,NSIGNA,NHLECO,XPACE,NSHirT, 

I  THPSLP,RLCNPt ,NECNP2 

CONNUN  /1CRUS/1CRUSS,HPR 

F.NNNSF(CNS,GAaNA)sSURT(  (  (2,a«(ENS*«2n/(CANNA*1.0))/ 

I  (1.0-I(GANN4.I,0>/(CAnnA»I,0I)*(Cns**2))  ) 

aTFLNS  (tNS,GANMA)aSuRT  ( 2 . 0*GANNA/ (GANha « t . 0 ) I • 

I  (F.NS/(l.U>((GANNAa|.0)/(GANNAFt.O) )•(ENS••2)>) 

NOEFLT  «  1 
NCASFall 
N  NCASlBO 

10  IFIHCASI .EO.NCASL)  NCAStaU 

*»»»»RF.AO/aRITF.  BASE  PRESSURE  CASF  INPUT  nATA»»A** 

CALL  INOUT 

iriNCASF.FQ.U)  GO  TU  B 

••••LIHITINC  KAOll  FUR  (I)  ANU  (F)  STRFAHS  ARF  SPECIFIEU  HERE. 
RLiat.S«RlE 
RLE>U.S*RII 

•  ••««INITI  ALIZATION  OF  BASF.  PRESSURE  ITF.RATIUN  LOOP. 

OTKRUlB(l.n>TROEl  1/2. U 

RPRaO.SO 

HPRLaO.O 

«*««»;ap|P|CAL  separation  PRESSURE  RATIO  F:xPR|SSinN  FHOh-.> 

ZUKOSRI,  AIAA  JOURNAL.  UCTUBEN  INOT,  VUL.  S,  HO.  10,  PP. I  TAN. 1 75 J . 
PRSEP  a  1.0  *, 0.1bS«(NACH  NU.I. 

•  aaaap.xrERNAL/lNTERNAL  I  LOnS  SFI'ARATIUt>  PRFSSURF:  RATMS. 

PRSIF.  a  1,0  *  0. JbS*FHNlE 
PRSII  a  l.n  «  0.3b5*FNNll 
HPRR  a  PRSIC 

IF  (t<PNSll/PR51F.|»PRlf IFI  .LT,  |,0|  HURR  s  PRSIIVPRIIIF 

NOSOLNaO 

NOSNAX  a  in 

IHPRal 

IBpPNiialS 


B-1 


nn  r>  nn  n  o  n  non  r%  n  n  nn 


NT»Pt«l 

ir<»HS  (  NTrrFeJ 

20  ir  (1HP«  ,l.t.  IhPHP*)  r.o  TO  40 

i.HITt  (J,22)  Hl’PL,  RPH,  HPMO 

22  FOMMAT  (//,  ISX, 

1  SIN  NU.  UF  HASk.  PH»SS.  ITCRATIU'tS  EXCFKUFO*** 

2  ISXtlOH  ••♦RPRL  ■  ,FT.4,7X,7H  PPR  «  ,FT.4,2X, 

1  >H  bPRR  s  <F7.4,4H  ••*,/  ) 

1F(  C  AHS(HPR«HPRM)  .LF.  I  .UC-3  I  .HR  .  (  Ht>R  .C.T.  RPRN  )  )  4RI  TFt  3, 24) 

24  FUP»A1  (lbX,33H  •*«  PHUPARLF  F|.0«  RFPARATIO'^  FUR  , 

1  2om  specified  UAIA  *•*  ,/) 

WRITE  (1,26) 

2b  FORMAT  (I5X, 

GO  to  2b0 

«*«**CHECX  that  RPR  IS  IN  Tilt'  SUl.HTtUN  RAOGF,  ( MPRD ,  RRRR ) . 

40  IF  ((RPR  .GC.  RPRL)  .ANO.  (bPR  .LF.  HPHR))  CO  TO  SO 
bPR>(t>PHD4HPRR)/3.0 

•***«CALCIILATe  THb  EXPANSION  PRCSSURE  RATIOS  FOR  THE  RJUNOART  CALCS. 

SO  PRRIE  a  HPR 

PRROIF.  a  I)PR*PR1UIE 

PHROlaPABOIFaPOIFOI 

PRRIlapRBOl/PRlUl 

PNBEa(PRB01E*PR0l0F.)/PRtae 

PReutaPHBOte«PRUtaC 

CPa3.0«((RRRE-l .0)/ (CANNAE* (CNNE**2))) 

CO  a  -CP*((R1E**3-N|I«*2)/RE**2) 

*««**«RITC  THE  CURRENT  IRIAL  SOLUTION  DATA. 

CALL  OUriM(lBPR,A,ENN|l.PRIUI,PRROt,PRH|I,PROEOI,THJCI,PRIIE, 

1  CNN|E,PR|01C,PRBUIt,PRRlC,rNNC,PRE0C,PRR0E,PR01C, 

2  PRHt,NPRINT,BLDHU,F.KGRU,NSMAPE,NRCNP) 

•***«THr  internal  CONSTANT  PRESSURE  BNORT  IS  CALCULATEU  FUR  (PB/PUI). 
CALL  ACPHS(GANNAt,rNSIl,PRBUI,XIt,RII,RETAl(,RLt,lRPR,NPTSI, 

1  NPRlHT,l,LlHiri,RPTt,NSHAPEI 

•**««THC  external  CONSTANT  PRESSURE  RNURT  IS  CALCULATED  FOR  (PR/POIE). 
CALL  ACFBS(GAHHAE,ENSir.PRH0tr,XIC,RlE,8ET*lE,RLE,IBPR,NPTSE, 

1  NPR1HT,2,L1MITC,RPTC,NSNAPC) 

****«IF  IMPINCEbcNT  occurs,  THE  IwpiNtiENbNT  POINT  AND  THE  FLU* 

PROPERTIES  POWNSTREAM  OF  THF.  RECUPPRESSION  SHUCK  SISTEH  ARE  FUUNU. 

CALL  CROSS ( GANNA  1, BPTI, LIMIT  I, CAnMAE.RPTE, LIH I TE, 

1  N1C,NEC,NST0P,TJ“LI,TJNLF,PPSHI)K,MPRMT) 

IKRtCONPapRSMOK  .LT,  l.O  .ANO.  N5T0P  .KM.  1)  NST0Pa2 
CO  TO  (60,12,1)4),  NSTOP 

CALCULATE  THE  BOONORT  LAIER  woNFhTUN  THIKNKSS  AFTER 
THE  BXSE  EXPANSION  IS  COaPLETFO  USING  THF  RESULTS 
OF  NASH  PUBLISHEU  IN  R.b  H.  NO. 1 344 , 1 Nb }. 

80  BLNTCRaRLNTE*(EaNlE/)NNNSr(RPTF ( 3, NEC ) ,G4NM At ) ) * tl/PHRl  K 
I  •■•TrLNS(F.NSIt,GANNAE)/«TFLNS(RPTE(  )  ,  NEC ) ,  Gann  At  I 

RLNTIBb8LRTI*(ENN1  I/FNNNSF(BPTm,NIC)  ,«AWMAI  ) )  taX/HRMl  I 
1  «BTFLNS(ENSt I,CAMMA1)/WTFLRS(MPTI( 3,N1C),CAWNAI ) 


ITCaO 
CO  TO  NO 

***«*NI)  INVISClII  SOLUTIUN  trial  CASfS. 

btlWHI'H  0)  Nil  .SULiITIOn  fplAl.S  a  NIi»n*)i, 


uuuo  uuuuuu  u  uuu 


c*(«**Nn  soLuriUN— -;*o  ixpiotKorwi  hr  ii*»i'«i5sin>f:  shjck  S'hjii  ihm. 

82  HCHRSHPM 
cn  TU  88 

C***«*MU  SOt.UTI(JN«>-.SHijCK  &Y.STER  Ixtf.SRI  EXIST  » (JR  THML  »*bUE  UP  8PH , 
84  HPKL»«PR 
8b  HPRx(HPMU4hPPH)/2. 

NUSULN3N(lSn(.>U  I 

If  (NOsni,N.i.c.Nns<**xi  UO  111  2»> 

C»»»**H»XIHU*<  NUHhER  UP  IIU»S|iLUTIOR  TPIXLS  Fxrrtopo. 

c 

■NITP  (1.881 
88  P(tPi*T(//, 

1  l)X.4fH  ***i>XX1m(Im  ku,  up  NU  SIlUUTIliN  TRIXI.S  LXCP.LliP.D**  ,/. 

2  15X,4»H  •»»».»•.♦»*»•♦♦♦♦«••*♦••♦••••••»»•**»*»»»»»»**•  /) 

cn  TO  2bb 

c 

C*****ST«RT  HASE  pressure  and  TE'IPERATURE  PATin  irERAftilN  LUUPS. 
yn  TKHOlxTROCI 
IE«1 
HE»1 

100  TRHUt«TRPOI/TRnEl 
ITC»ITC*I 

C*«***CALCUI.AT10fi  AND  OUTPUT  UP  TURRUt.FI.'T  MIXING  HESUtTS. 

CALL  rjMIX(CAMMAl.GCI.UPTl(i.NIC),TRR01,TJNkl. 

1  GAMMAE.GCC.BPTEt  ),NF.C}.THHOf;,TJNI.E, 

2  Rll.CNSl 1 .BETAIl .HPTl(2,Nir).PRSH0R. 

3  PUIPUI.TROCI.RCCUMp.HLOR.PNCR, 

4  BLMTIB.BLNTEB.ITC) 

CALCULATE  THE  nilUNbRY  hLELO  RATES  RATtUO  TO 
THE  NUZZLE  MASS  FLOh  RATE 

nLOHOCsSLMTER/Rl I *HiP/Rll«PRRU4(1.0«C0S(RETAlI  ))• 

1  SOHT(GCI/GCC/rROEI  )*«TFI.*'S(HPTEt  3.MEC)  .CAMmAE)/ 

2  wrPLMS(E'NSll.CAMMAl) 
BL0RUlaBLMTlH/Rtl«PRBU«(l.0«CD8(HCTAlI)  I* 

1  NT»LHS(nPTI(3«NlC).GAMMAl)/Mrri.MS(PMSIl.CAMMAn 

HLDRU>HL0R0f  fRLUHOl 

PNGRaF.«HLURUC*CAHNAE«(GAMMAlol.O)*TRnEI*(CCE/GCI)/(GAMMAl* 

KGAMMAE'I.O)) 

ENGRUlaRLURQI 

CNGKUaENGRqttCNGHOI 

CALL  0UT2H(PRMe.PRRIl.PM0E(ll,TRi»Uf .TRBUI.THUEI.PHOIE  .PRlIk. 

1  BLUK.r.RGM.NPRINT.CP.CU.HLIlRU.PNCHO) 

»**»*SET-UP  ITERA110H  LOOPS  TO  PINO— 

NTYPEal  (NOblSnERPRGETlCl.  TPHOI  SO  THAT  ENCRbENGRU. 

NTYPEB2  (NONlSOfNCRCETIC}.  TRROl  SO  THAT  PLDRaHLDRU. 

NTYPPal  (ISUCNECETIC),  CONTjRtlP  TU  BASE  PRESSUHP  ITERATION  LOOP 

TO  PINO  BPP  S'l  That  Bl.nRNHLURO. 

GO  TO  (124,126,710),  NTVPC 

••«a4TRBOI  ITERATION  LOOPS  FOR  THE  NON» I SOFNKHGETIC  CASE. 

124  VARa(CHGM«eNCRO>ENGRU) 

CO  TO  IJO 

126  VARa(aLDM*HLORU>hLUR()) 

130  GO  TO  (140,142),  Nt 
140  UATA(IC,l)aTRBUI 
0ATA(lE,2)aVAR 

•  ••••iteration  for  IRBOI  such  that  PI.r.Rsf.NGRM  C«  HLiIxaifl.llRij, 

(NUIE  that  TRHOI  is  MESIHICTFU  Til  the  range  (TRUEI.I.O)  ) 

f 

142  CALL  lTERtTR»OI,0T».M(iI,1.0f.-4,1,*>,YAli,il.o,  I  .  Uf.-S ,  I  ,  (f. , 

1  THHII)  N,  YARN.  rRHOIP,  VAH'.'ISO  Y1  .I'SuNY/  ) 


rtrionn  ono 


IK(ThHUi>l.u)  isn,  isti,  |f>o 
IbO  CO  TO  (too,  I  00,2"*' I  <  Nt. 

•  ****FXrR*PUL*TIUN,  IF  NtCtSSAHY,  FOR  Te»Tf.HATIIRK  HAltO  THMUI 
SUCH  THtr  FHCRsKNCHll  OK  HLDKshI.ORO. 


M2 


lts|S-l 

IF(ABS  (UATAd  ,2)  )-AHS  (  0  ATA  (  I F  ,  2  )  )  )  17li,i;0,IM0 

l«l 

1IS2 

GO  TO  lOO 

l=lt-l 

II»1F 

NATlUs(OATA(lI,n>OATA(I,n)/(DA1A(ll,2)-l)ArA(l,2l) 

IKbUlaUAI  A(  1,  I  )-HATI0*IIATA(1,2) 

GO  TO  (202,2041,  NTTPF 

T»EO»Tt<eOI 

MYPK«2 

(;ii  TO  90 

TRHtisIHBOl 

NTYPCal 

•  F.NO  TRHOl  ITtRATlON  LliOPS. 

•CUNTlNUC  THF  BASt  PHCSSURL  RATIO  (BPR)  1TF.HATI04  LOUP  TU  FIND 
HPR  SUCH  THAT  OVARaO. 

«FUR  the  NON'lSOtNERCEriC  CASF. 

UVARa(TPEU-TMbU) 

CO  TO  214 

•FOR  THE  ISOENLhGETIC  CASE. 

OVARa(RLURU«aL0R0>HL0R) 

SIGNbUVRM/assiovAR) 

IFCSICNI  2111, 2IB, 222 

BPKRaBPR 

GO  TU  22R 

bPHL«BPR 

IF(lbPR«l)  210,220,234 
DUPRa(BPKH«aPRL|/2. 

GU  TO  2  It 
SIGNal.O 

ObPR««((BPR«epRn/(OVAR>OVARt  ))*()VAR 

HPRtBBPR 

OVARlaOVAR 

•ITFRATIUN  FOR  BPR  SUCH  THAT  DVAPSO. 

CALL  1  rCR(BPR,nBPH, I . nF<*4, SIGN, bVAR, 0.0,1 .OE-5,1 RPR, nbPR, 

1  HPRN,DYARN,bPRP,PVARP,l<SGNH| ,NSGNH2) • 

GO  TO  (20,20,242),  NBPR 
•SOLUTION  FOUND. 

GU  TO  (2bO, 250,254),  NTYPE 
•mrite  solutiun  data. 


250  aRITE  (3,252) 

252  FURNATI//,  20X,  J2H  •••NONoISUENEPRETIC  SOLUTION***  ,/, 

1  20A,  )2H  *******************************  ,//) 

GU  TO  25t 


254  MRITC  (3.256) 

256  FURNAT(//,  2T«,  28H  **«ISnttiCRRETiC  SOLUTION***  ,/ 

1  2TA,  28M  t*********************^**  /) 

C 

251  CALL  UUT3J(PRBE,PHHil,PRUCUI,TRPUt.TRH(M,TKilF(,PHI)IE  ,PNI1E, 
I  BLUR.ENCP,  1  ,CP,CII,bLI>RII,F‘<f;Rli) 

IP(NP)JNCN)  2R5,2R5,2T0 
C*****PUNCM  SOLUTIOm  .data. 

260  IF(NPUMCH)  2A5,2N5.265 


iOi  tl.iH)  I'RUlt.  i'PHh 

2b7  mh>««T(^n  1  ,4,5«  .  I  IMNM  SULUTliiN,  5X,  TH  pn/PR  s  PM.S) 

GO  TO  2HW 
C 

270  RlIFsHll/Ne 

C»»*»»CT*--T/0*  (  thrust  CflfFnClfuT). 

CT  *  (  <RU7,»*2)/CO.b»GAHHAK»(f.HNlf*»2)  )  )  •  (HMl  It«(  I  .I)*CAH»*A1» 

I  (fc-«<ll»*2j)-l.(>) 

C*»*»*HHf----JFl-TfJ-FHEKSTHtA!*  nOH^eiTUM  H.U*  HATIU. 

HHF  «  l(;AMPAl*(F;nNll*«7)*(KltE**2)*(’K1  lt)/(GA«^AF:*(KUML**2)) 

C 

HH1T7  (7,272)  HHOI t , PR  11) , HHHF ,CP , CO . H-F ,CT 

272  FIJRHAT(2Fi  I  .4,5F11.S) 

C 

280  IF  (NCASl  .KO.  NCASt)  HHITF  (T,2H2)  < A ( 1 ) , 1 ■ I , 20 ) 

282  FORMAT  (  20A4,/,80H4)*«**»«»*»»+«*4**»44**)***»****»*****««»***t 

144444444444444444444444444444444444) 

2H5  IF  (RCASl  .Gt.  NCASt)  CU  TO  290 
C 

C  »«»**C0  TO  NEXT  CASt4»»»» 

GO  TO  I U 
C 

290  CALL  bttPdO.S) 
oRITC  (b.lUO) 

300  FORMATC IHS, >00  TOO  »15N  TO  RAFF  ANOTHFR  RUM?  TtSsl.NOaO') 

RtAO  (S.*)  IRON 

IF  (IRON  .CO.  1)  CO  TO  10 

ENU 

C« t**)*********#**************************** •**•••»*•»* »•***•»*«*»*•• 

c 

SIIRROUTINC  I1UT1M(  1,A,)HNI1,RRI0|,(>RR0I,PRR11,PMUC')I  ,TR0E1  ,PR1  It, 

1  tMNlt,PR10tC,PRB0IF,PRt»lE,CMNC,PRt0b,PRB0t,l'R0lC, 

2  PRHC,NPRlNr,(lL0RO,t'NGRU,M5HAPE,NRCNP) 

C 

C»»»*»S0RR00T1MC  MRITCS  OUT  HEADINGS  AND  CORRFNT  DATA  USED  FOR  THE 
C  INVlSCin  FLOH  FIELD  CALCULATIONS. 

C 

C  AttVARlAHLCS*** 

C 

C  1  a  I>TH  VALUE  OF  THE  INPUT  HA5E  PRESSURE  RATIO. 

C  A  a  heading  card  DATA. 

C 

C  FOR  EITHER  STREAM  AT  (II),  (IC),  OR  ( F -‘•FREESTRt AN ) . 

C 

C  CNN  a  MACH  NIINHCR. 

C  PREO  a  PRESSURE  RATIO,  PE/PP. 

C  PRIO  a  PRESSURE  RATIO  AT  (I),  P)/P'«. 

C  PRRO  a  BASE  PRESSURE  RATIO,  PP/RO. 

C  PRBI  a  BASE  PRlSSURC  RATIO,  PB/CI. 

C 

C  PRIIC  a  INPUT  STATIC  PRESSURE  RATIO  PF  STREAMS,  PU/PE. 

C  THUEI  a  STAGNATION  rFMPERATURE  RATIO  OF  STREAMS.  TUF./TOI. 

C  PROCUla  STAGNATION  PH4:SSURE  RATIO  PE  STREAMS,  PUE/POI. 

C  HPNINTa  SEE  SUBROUTINE  *lNUUTa. 

C  RLDNU.ENGRO  a  SPtCIFlFU  VALUES  PE  THE  RI.EE  U  AnU  F.  4EKGT  RATIOS. 

C  NSHAPLa  0,  NO  HOATTAIL. 

C  a  1,7  OR  )-->UGIVt,  I'ARAPOLIC,  OR  CPMICAL  M'lAI  TAILS. 

C 

DIMENSION  A(20) 
lE(NpRINT)  I0T,|0T,N9 
C 

99  mHITE  0,100)  (A(J),.la|,2'»),N“r*’U,HHl|E:,rMiiF.l,PRIlE.0i, 

I  PRUIE  ,PLURO,E.NGRP,  I 


|(i0  rUHH*T(lH|,  S*,  70*4, J». 9m  mmCKP  *  ,J!,//, 

t  ISX.StH  ••♦•♦*»»»l»0-SIKt*M  P*S»  PHtS.SUKF  PMUGh * 4* •»••»»♦* ,  //, 
2  271, 25H  ♦♦•♦••CUMHtNT  ll»l**t*»»*,//, 

J  15X,UH  PJI/Pt  *  F9.4,17X,||H  TOF/TOl  *  F9,5,//, 

4  15X,HH  POt/POl  3  f9.5,l7X,l|M  POf/Pt  s  FH.J,//, 

5  |;x,9H  HLUPU  a  L12.i,  IPX,  4h  F'JGPl)  *  FI2.5,//, 

b  27X,]1H  1P[*L  nXSt  PKFSStIMC  HXTlH  un,  ,14,/, 

7  22X,JJH  »♦•»»  ••»♦♦»*»  •»»••  •••  ,//) 

C 

"HITE  (3,101)  EHHt  I  ,PF  mi  ,PHP(lt  ,PHUl  1  , 

1  EHMf  ,PP10|t,rPM0ir,l'HBlt 

101  FDR"*! (2HX,22H  •♦•lNTtHH*L  STMF*H*»*,  //, 

1  15X,tlH  LHNII  ■  F7.4,2SX,I0H  P11/P(iI  *  FR.b,//, 

2  1%X,9H  P8/PU1  a  rH.b,23X,9H  Ph/Pll  a  FH.b,//, 

3  2RX,22H  **«EX1FHHAL  5THF*H»*»,  //, 

4  15X,IIH  t.“HlC  a  F7.4,2bX,llH  Plt/PIUF  *  FR.b,//, 
b  lbX,9H  PB/POlF.a  r8.b,23X.9n  Pb/PIF  a  FR.b,//) 

C 

"HITE  (3,102)  tMNF,PRFOE,PRBOF,PHME 

102  FOHM*1(30X,17H  *»«FRtFSlHt*M»**  ,  //, 

1  lbX,7H  EMHE  a  F7,4,  23X,  9H  PF./POF  a  FB.b,  //, 

2  15X,  9H  i>n/PUE  a  F8.b,  20X,  9H  PB/PC  a  Ft.b,//) 

C 

IF(NSH*PE)  103,103,10s 
C 

lOJ  HHlTC  (3,104) 

104  )URh*T(21X,32H  44*  NO  HUATTAlt  BFFORE  BASF  *44  ,  /) 

HETUHN 

C 

105  aRITE  (3.10b)  NSHAPC 

106  FaRBAT(2SX,  27H  *4«  BUATTAIL  —  NSHAPE  a  11,  4M  444  ,/) 

C 

107  HF.rURN 
LNO 


SUBROUTINE  INUUT 

44444SUPRUUTINC  REAOS  IN  THE  INPUT  DATA  AND  THEN  CALCULATES  THE  INPUT 
DATA  FUR  the  MASTER  PROGRAM.  THE  IDENTIFICATION,  HEADINGS  AND 
INPUT  DATA  are  then  kHITTEN  OUT. 


444VAR1ABLES44* 


FOR  EITHER  THE  INTERNAL  (I)  UR  EXTERNAL  (E)  STREAM 


BETDI 

BET(>3E 

BLDRO 

BLMTE 

BLNII 

ENCKO 

ERNE 

EMSll 

GANNA 

GC 

INDPT 

kPRFSR 


a  FLO"  ANGLE  (IN  DEGREES)  AT  (XI, R|).  CC«  IS  POSTIVE. 

(  BETDll  IS  (4)  AND  BETOTE  IS  (*/•)  ). 
a  INITIAL  AFTERBODT  ANGLE  AT  (X2E,R2E). 
a  SPECiriEO  VALUE  OF  THE  BLEED  RATIO, 
a  EXTERNAL  STREAM  BUUNDRT  I.ATFR  MOMENTUM  THICKNESS 
a  Internal  stream  huundrv  laveh  momentum  thickness 
a  SPECIFIED  VALUE  OF  THF  ENERGY  RATIO, 
a  EXTERNAL  FREESTREAN  MACH  NUMBER. 

8  MAC)!  STAR  AT  (11). 
a  RATIO  OF  SPECIFIC  HEATS, 
a  GAS  CHNSTAnT  ( LBF-Fl /Lhm-R ) 
a  1,  input  is  RFAO  FROM  NAKELIST  DATA 

2,  INPUT-  IS  INPUT  ,RT  TF.HhINAL  UitPINr,  PROG.  EOCUTION 
a  II,  PRIIE  IS  INPUT,  AND  PhOIF  IS  rALCULATEU. 
a  I.  PHUIE  is  input.  A'O,  EuHE  IS  r  *  LCUIA  Tl-jii . 


r>r»onorinnnooooonnnoooo/»r;o^nnr>orto<^nor>nooor»r)nr»r»orir>r>onoonr>nor>r>nnnr»r> 


HC*SE  *  MUNHfB  Uf  I'WKSSUKK  «-*11US,  Pll/Pf  ,  f  UH  «HICH  H*St 

HPESSliHE  CAECHLATIOMS  »Nf  TO  PE  “*nE  fOU  *  GIVEN  SET  OE 
CUNOITtUNS  *NO  CEONtTPT. 

NDEFLT  »  0,  THE  VAH1»MLES  »Kt  HtStT  TO  THE  •t'Er»ULT  CONE  1 GUH AT I  UN* 
AFTER  THE  CASE  (SET  Of  ('RtSSOHE  HAFlUS)  IS  COHPLETFO. 

«  I,  the  VARIAHLES  >1LE  NOT  PK  FFSI'T  AT  UPON  CUMpEETION 
Of  THE  CASE. 

NOTE  CHANCING  THE  VALUE  tiF  •NtOFLI*  nILL  FIRST  AFFECT  THE 

CASE  SUCCEfPJNC  THE  CASE  IN  1.HICH  IT  IS  CHANGEO. 

NPONCH  a  0.  SUMHART  OUIPUT  DATA  nut  ruNCHFti. 

NfRINT  a  *1,  INPUT  UATA  ANU  HASP  PRFSSURF  SOLN  PRINIEU. 

NSHAPC  a  0,  NU  AFTERHUOV. 

a  1 ,  OGIVE  BOATIAIL. 

a  7,  parabolic  BOATTATL. 

a  3,  CONICAL  BUAIIAIL  UP  FLARE. 

PRIIF  a  STATIC  PRESSURE  PATIO  Of  STPFAnS,  PII/PE. 

PPUEUI  a  stagnation  PRFSSUPC  RATIO  OF  STREAMS,  PUC/Pi'l. 
ppuie  a  internal  stream  stagnation  pressure  TU  EATERNAL  strear 
STATIC  PRESSURE  RATIO  tNOZZLf  CHAMHER  TO  FHEESTREAN 
static  pressure  RATIO),  f‘OI/PE. 

HECUNP  a  RECOMpRESSlON  COEFFICIENT  aO.O 

NOTE  •—  IF  THE  INPUT  VALUE  OF  RECOHPaO.O  .AND. 

1)  NSHAPCaO,  THEN  KECONP  IS  CALCULATED  FRUR 

EMPIRICAL  COMATION  UN  CARD  NO,  1NOU2630. 

2)  NSHAPE3l,2,l,  then  RFCUMP  IS  CALCULATED  FROM 

mHITE'S  EMPIRICAL  COKRELATtON 
BASED  UN  Ff A  DATA 
DELTA  a  BOUNDARY  LAYER  THlCANtSS 
SlGMFl  a  FACTOR  FOR  SIGMA-1 

sicMr2  a  factor  fur  sigma-e 

NHCmP  a  0  empirical  RtCUMPHFSSlUN  CUEFFICIENTS 
mHITF.  RESPECIIVELT 

a  1  RCCOHPRESSION  CRITERIUN  AFTER  PAGE, APPRUAIMATED  AT 
KURST 

a  2  CORRECT  PAGE  CRITERION,  CALCULATED  ITERATIVELI 

(XPAGE  a  PHESCRieeO  FRACTION  OF  TJNLl  FUR  DEFINING 
RECONPHESSIUN  BEGIN) 

a  3  UNCRA  angular  REATTACNPCNT  CRITERlUN 
a  4  ONCRA  CPIT.  NUOIFIED  FOR  Tao  STNEAN  REATACHHENT 
a  S  AS  4,  BUT  SECOND  ORDER  CORRECTION  FUR  PHlu 
a  6  AS  4,  BUT  NUNLINEAR  TRCATNENT  (HCCOMPat.  NECESSARY) 
NSIGNA  a  0  SIGMA  AFTER  KORST  AND  1RIPP 
a  1  SIGMA  AFTER  CHANN APRACAOA 

NSHIfT  a  0  ORIGIN  SHIFT  AFTER  SOLIGNAC  ANO  SIHIEUA  (  NBLF.ED .  CT .  0  ) 
a  I  ORIGIN  SHIFT  AfTFR  HILL  (NPLEFD.CT.U) 
a  2  ORIGIN  SHIFT  SUPPPESStl' 

RDSICN  a  STING  DIAMETER  FUR  NOJET-CASF  (PRATIU«0| 

SEPRtS  a  ZUROSAI  CONSTANT 

TROCI  a  STAGNATION  TEMPERATURE  PATIO  OF  STREAMS,  TUF./TOI. 

XI, R1  a  CUUR0INA1LS  OF  POINT  «Hf.RF.  SEPARATION  UCCUM5. 

(RI'S  ARE  POSITIVE) 

X2K,N2Ea  INITIAL  CUURUlNATF.S  DF  THE  AFTEHAUDT. 


***PRrSSURL  RATIO  MUTATIUR*** 


PRIOI  a  P)]/P01,  PRllt  a  »’)J/PE,  PROlt  a  Pni/PE, 

PRIIIE  a  Pll/Pir,  PUltOI  a  POIE/POI.  PPIUIE  a  PIF./PUIF, 

PRUIUE  a  POIE/POE,  PREuE  a  PF/PUC,  PROtUJ  a  HUr/PIll, 

PRBIE  a  PH/PIF,  PRK(>|L  a  PR/IMilF,  I'RHil  a  HH/PIl, 

PRAOI  a  PH/PUI,  PRHE  ■  PR/f f ,  PRItE  a  Hll/Pf. 
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COHMON  CHAMI,  CHAk^,  Pt,  P?,  PA 

CUXXON  /UATAin/  GCI  .CiA'^wAl  ,F>*SI  I  ,  XI  I  ,kl  I  .litTA  1  t , 

I  GCt.GAk-AE.rusit.XlF.KlF.ilFrAlt.PkOlUE. 

i  TWIJFl  ,PR1  le.kfCOHP,  A,t  kkl  I  ,  PR  I D I  ,  F.  «lk  1 1 ,  PR  I  H  1  t  , 

J  NPP1i<T,NCASI  .hCASt.MLI'PlJ.mCHO.Rt.EkHF,  .PHFOF, 

4  NPUkCH.pPliFOI  ,HHI)IF,Pnit  J1  ,  N3HAPF.NPTSF,  PRI  lU  . 

5  «UtFLT,l*l,»TF,PLPTl 

COMMON  /C*  I  /COBT.XFl  X  ,PF|X,NStF  ,  StPW  1 S  ,  A  TFS  F ,  EF.E ,  St  PO  ,  OFL  f  A 
COMMON  /Cm2/X2E,K^E 

COMMON  /C«J/SIGMF1 .SIGmf  2,Nf<CMP,NSIGMA,NhLtED, APAGE.NShlFT.THPStP 
I  .RtCMPI ,HECMP2 

common  /1CP()S/1CH0SS,BPR 

niMENSlUN  PmB(IOO,S,2),  CHARKS.iO),  CHANElS.jO),  Pt(S),  P2($), 

I  P)(5).  A(20) 

FMnmSF<EmS,GAmma)«S<JRT(  (  «2.04(fMS»»2)  )/(CAMHA»l,i1) )/ 

1  ( 1 .0-< (GANMA-1 .0)/(Camma*1 .0) )*(fcMS»»l) )  ) 

tMSMNF(tMN,CAMMA)3SOPT(  (0.b»(GAMFA4l.0)»(F.-'N»*2)  )/ 

1  ( 1 .0*<i.5»(GAMNA«|.0)*(F-N»*2)  )  ) 

PRMNr(EMN,GAMMA)c( I .P»( (CAMNA-I .0)/2.(>)»(FMN»»2) ) •» 

I  (-CAmmA/(CAHMA>I .0)  ) 

IF  (NCAsi.Nc.o)  GO  TO  eu 

GCEsSJ.33 

CAMMAEmI .4 

KLCUMPaO.O 

NPUNCHan 

NPPlNTa-l 

ntOPOaO.U 

ENGPOaO.O 

SlGMElBl .0 

StG>F2a| .0 

NSHlETaO 

NAmEL1ST/DATA/NPH1nT,NSHAPE,XII,PIE,X2E,P2E.RCT02E,FMNE, 

I  XII  .Pll .BLTOtt ,GCI,GAMMAI,FMNll,THl)CI,KPPtSR 

INTtPACTIVC  INPUT  or  DATA 

NPITEIS.Il  ) 

71  rOPMATISX, 'ENTCP  INPUT  OPTION  CHOICE,  laflUE  2*TERmIMAL') 
BEAO(5,72)  INOPT 

72  FORMAT  (11) 

MRlTE(5,b9) 

b9  FOPMATdH  .•ENTFP  NPCMP  PLCOKPPF SS ION  C HOI CE '  ,  /  ,  5X , 

I'Oa  empirical  NLCOMP  COEFF  AfTEP  ADDY  OR  MR I TE ' , / , bX , 

2'1«  KOHST  APPROXIMATION  OF  PAGE  CR I TF R I  ON • , / , bX , 

J'2a  CORRECT  PACE  CRITERIA  CALCULATED  ITER A T I VEL Y • , / , bX , 

4'JB  UNFPA  angular  REATIACHMENT  CR 1 TER I  ON • , / , bX , 

b'4a  ONERA  CRITERIA  mOOIFIEO  FOB  2-STRFAM  PEATT ACH * , / , bX , 

N'ba  AS  4,SEC0Nn  ORDER  COHHECTIKN  FOR  PHIl)',/,bX, 

7*68  AS  4, NONLINEAR  TREATMENT  ( NLCUMPal . NF C . ) ' ) 

HEAD  (9,13)  NPCMP 
IE  (NRCMP.NE.O)  GO  TO  7 
NRITC  (9,9) 

FORMATCIH  ,'CN1ER  RICOMP') 

READ  (9,6)  RECOMP 

rORMAT(F6.4) 

rORMAT(20A4) 

FORMATIEF 10.4) 

FORMATdH  ,*XXXXX.XXXX') 

FORMATdl) 

FORMATdH  , 'ENTER  ALPRANUHtRIC  HFAUING  •  7'IA4') 

CONTINUE 

attirC(9,l4) 

RE4(>(9,I0)(A(|  },  la|,2'>) 


»  u  n  K  ( b , « ) 

.•f.'TCW 

MhAris.'JJ  N.sniFf 
Fiiw«*r(ii) 

IF  ilNUPI.Fg.n  c; 

•»»»»««*iNPUi  tjr 

FIIH^AfllH  .'AFIFH 
IH  ,  '  0 

M'P“*T(IM  ,•  I 

F0W«AT(1H  J 

FUBMATdH  3 

FORXATIIH  .'CMtH 
oRITtlS.  1%) 
wRlTtlS, 16) 
»RITH5,17) 
••RMECS.IU) 
•<R1TE(5,I9) 
kRl IF(6.2n) 
REA0(6,1 JjNSHAPe 
FORRAiriM  ,*POOI 
FCtRRATIlH  ,*THF.lf 
rdROATtlH  ."X*  0 
RRITL(b.21 I 
■R1TC(S.22) 
•IH1TC(6,2J) 
F0P“AT(1H  ,'ENfER 
«R1TF(6,24) 
i<RITE(S«l2) 
REAOCS, 1 t )XIC 
FORNATdH  ,'F**1ER 

HR|  rf:(S,25) 

RR1TF.{5,12) 

REA0(5,11)P1E 

IF  (NSHAPE.EO.O)  ( 

FURMATdH  , ‘ERICH 

■RITE(5,26) 

«HlTe(9,l2) 

HE*D(5,1I  )X2e 

FORMATdH  .'ENTER 

•<R1TE(5. 27) 

xRITEIS.d) 

HL*n(9.inR2C 

FORMATdH  , 'ENTER 

HR|rC(S.28) 

NKITC(5,I2) 

RE«0(5,U)RCTn2E 

FORNATdH  .'ENTER 

NRtTE(6.29) 

«RITF.(6.12) 

REAUIS.inFHNE 

FORNATdH  .'ENTER 
NRITEIb.JO) 
■RITC(b.i2) 
PCADtb.ldXII 
FORHATdH  .'ENTFR 
NRITCCS. 32) 
■R|TE(9,12) 
REAUlb.lDRII 
FORNATdH  .'ENTER 
«RITE(b. 13) 
NHITLIb. I  2) 
REAOCb. I t )HFTU1 I 


■  SHIF  1  •  ) 


g. ) )  gu  ru  70 

Ul  EXrPHAL  PARA-FTFR.*; 

•AFIFHHUl'I  SHAPE  PARAHETF  PS:  '  ) 

'  0  -  cri.ir.OHiCAL  Af  Tf  PH<tii»  • ) 

•  1  -  ncivt  HOAT7AIL') 

<  2  -  PAhAPOLlC  HUATTAIl') 

'  3  -  CONICAL  POAITAIL  OR  flARE') 

•ENTER  AFlERPOOr  SHAPE  PAHANFTER  -  ll‘) 


SHAPE 

'Pool  ANO  NOZZLE  PIHFHSIIIHS  ARE  RELATdF') 
'THF.T  CAN  BE  INCHES.  FEET.  OR  CALIBERS') 
'•X*  OINENSIONS  ARE  POSITIVE  AFf) 

'ENTER  ‘X*  AT  BOOT  BASE') 

IE 

'EN1ER  RADIUS  AT  BODY  BASE') 


EQ.O)  CO  TO  49 

•ENIEH  'X*  AT  START  OF  BOATTAIL') 


2C 

'ENTER  RADIUS  AT  START  OF  BflATTAlL') 


!E 

ENTER  SLOPE  AT  START  OF  HOATTAIL  • 


lCTn2C 

'ENTER  FNFL  STREAM  MACH  NUMBER') 


HNE 

INPUT  internal  PARANETFRS»»»*i 
'ENTER  'X*  AT  end  OF  NOZZLF') 


I  J 

'ENirH  NOZZLE  EXIT  HAhlUS*) 


•  B 

•enter  nozzle  tXFT  A'*OLt  -  PEH.') 


J 11 1,  ll« 


J'.'  ^L"!F  ' 


<4  FOHXATdH  ,'K*-rKh  laiZiLt  GAS  Cll'-STAHT  -  I  nr/l,it<*  «•) 

J5  FOKfArdH  ,•  53. J4  FUK  *IK‘) 

hHITt:(S.34) 

••NITt(5.35) 

«HITt(5,12) 
hbAUCS, 1 1 )GCI 

3^  FORMkTdH  .'FSTtH  NUZZLL  GAHk*  •  1.4  FUM  AIK') 

FRITb(5, 36) 

«UlTb(5,l2) 

HFAIXS.II  )GAmi4A1 

37  FOWMATdH  .•FNTEB  NOZZLE  EXIT  FACH  MIIMHKN') 

•I1»ITF(5,37) 

»R1TE(S,12) 

REArXS.  11  )rMNl  I 

38  FURNArdH  , 'ENTER  EXTERNAL  TO  INTERNAL  STREAM  STAGNATION') 

39  FORMATdH  ,  '  TENRER  A  TORE  RATIO  -  TOE/TO)') 

NRITECS.IS) 

••RITEIS,  39) 
nHITE(5.I2) 

REAU(5, 1 1 )TRUri 

40  FUNHATCIH  , 'ENTER  NUNbER  OF  CASES  •  11') 

7U  nRITCIS.AO) 

RCAD(S«1J)NCASC 
IF  (INOPT.EO.I)  HEAD(4,0ATA) 
ir  (INOPT.EO.I)  GO  to  79 

41  EURMATdH  ,'  ENTER  TYPE  OF  PRESSURE  RATIO  INPUT') 

42  FORNArdH  .'  )  FOR  INTEKNAl.  STAGNATIUN/EXT  STATIC  PRESSURE') 

43  FURNATCIH  ,'  0  FUR  INtERNAL  STATIC/CXTCRNAL  STATIC  PRESSURE') 

HH1TC(5,41  ) 

NRITe(S.42) 

>RITC(S,4I) 

REAUIS. 13)KPRESR 

44  FORNATdN  , 'ENTER  POJFPFd 

45  FORNATdH  , 'ENTER  PJ/PF') 

C*«**«CALCULATION  OF  PROGRAM  DATA. 

79  HETAII*  0.0174532«bCTDlI 

EMSII  a  EMSMNFIEMNII.GAMMAI) 

PRIOI  a  PRMNF(CMNII.CAMMAI) 

Oil  a  2.04R1I 

DIE  a  2.b«R)E 

XllOlt  a  XII/OIE 

EMSE  a  FNSMNF(emNE.GAMNAE) 

PPEOi  a  PRmnMKNNE.GAhmAE) 

RICl  a  Rtl/H|E 
IF(NSHAPE.NC.O)  GO  TO  56 

C«4*«*UNU0RM  EXTERNAL  FLOP  xlTHOUT  AN  AFTFRAOnV. 

RE  a  RIE 
ERNIE  a  erne 
ERSIE  a  ERSE 
PRIOIE  a  PREOE 
PKOlUE  a  1.0 
CO  TO  62 

C*AA«AAFTERBU0V  reforc  the  external  STRFAR'S  separation  point. 

56  BCTA2Ea0.01745I2*HETD2b 

CALL  ABTS(CAMNAe.EMSC.X2E.K2E,RETA2C,XIF,RiE,M5HAPE, 

1  1 .NPTSt.NERROR.COBTl 

C*««*4SET>UP  0A1A  FOR  EXTERNAL  STREAM'S  SEPARATION  POMT, 
XlEaCHAREd.l) 

RlFaCHARC(2,l  ) 

CSIL  a  CHAREO.D 

ERNIE  a  ERNRSFIFNSIF.GA^MAE) 

RFTAIE  a  CHARE! 4. 1 ) 

HETOIE  a  57 . 2R57 7 R5 • hi  | A  I E 
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PhlOlt  »  PB*«'.r 
PwniL'Fsl  .0 
PE  ■  H2l 
U2E  ■  2.<>*R2e 
X2EI>2L  «  X2E/U2E 
XIFD2F  a  XIE/n2C 
DPIE2E  »  tiir/02t 

C«»»»*KtCOMPPCSS10N  COEFUCIEPT  UfcTEP»lPXTlUN. 

62  ir(PtCU«P.CT.l.'it-OJ)  GU  TIi  an 
IF(6SHXPF .NC.O)  GO  TU  64 
C«*««*F0R  CFLINDRICXI-  XFTF RPOOI tS. 

RECOUP  «  .4H1  ♦  I.OahXRlCI  -  O.  R74  1 E 1  •  *  2  ♦  0  .  JO  }♦  R I F.  1 «  *  I 

GU  TU  RO 

C«««**FnR  BOXrTAlt.EO  XFrF.MhOniES. 

64  RXTIUa  ( R 1 1 /R I  F ) / ( ( JO . 0 -ULTn 1 1 ) / ( 20 . Of PK T02E > ) 

SPETAs  eLTDjepBF.TOl  I 
Caaa  -0, I H 7S*SHKTAf 0. bS 

IF  (SHETA.l.T.'R. )  EPPa-0.b25aSHFTA-2.bS 
IF  (SPCTA.CT.-4. )  tPPa-0.0J7SaSHETAfl.lS 
RECOPpa  EPM*HArinf O.S 
C 

80  MCASl  a  NCASl  *  I 
C(aa«*THAIiSFEH  UR  HEAD  NF.a  CASE  OATA. 

C 

IF(XPRE5R.GT.0)GO  TO  82 
CAUL  BEEPdO.b) 

«RITF.(5.4S) 

WR|TE(S>I2) 

REAPCS.I I tPRATlO 
GO  TO  88 

82  cage  HEtPdO.b) 
aRtTE(5«44) 
aRllE(5,t2) 

KCA0(5.1dPRATrU 

C 

88  tF(KPReSR.NC.O)  CO  10  90 
Caaaa«FOR  Pd/PE  <PK1IC)  INPUT. 

PRUEaPKATlO 
PROIEaPRlle/PRlUl 
GU  TO  92 

C«*a««FUR  POI/PE  (PRUIE)  INPUT. 

90  PRUICaPRATin 

PRllE  a  PMUlLaPRIOl 

CaaaaaCALCULATE  VARIOUS  PRESSURE  RATIOS  FRUP  NFN  CASE  OATA. 

92  PHUFUIaPRIoi/IPREOEaPRIlE) 

PUlFOlapHUrOiaPRUIOL 

PRI llEapPtoi/IPOIEOiaPRIOlCJ 

PPtCEaPR10IC*PR0inE/PRE0C 

aRITE(S,46} 

46  FORPATdH, 'ENTER  BOUNURV  LAYER  pUPENTUN  THICKNESS  AT  BODY  BASE') 
aRITE(5,l2) 

REAOIbflDHLMTE 

PR1TE(9.4T) 

47  FORHATCIN, 'ENTER  HUUNURY  LAYER  pOPENTUP  THICKNESS  AT  NUZZLE  EXIT') 
HRlTE(b,13) 

REAOIS.IIJBLPTI 
CaaaaapRiNT  CASE  DATA. 

PRITE  11,94)  (ACD.Ial.TO),  NHC-P, NCASl 

94  FURPATdHI  ,5X,20A4,2X,  ■  NRCPPa  •  ,  1 1 , 9X  ,  1 5HPRJHLE  4  fliPHER  IJ,//) 

C 

IF(NSHAPC.tU.O)  G"  TO  IRO 
GU  TO  dO0,t20,|4O),NSHAPE 
C 
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o  r»  o  o  r>  o  r>  r>  r»  o  r»  n  r»  o  n  o  o  on 


100 

NRlTf 

(3.110) 

1  IN 

FORMAT 

(29X,21H 

•  ••OCl 

VK  BOATTAIL***  //) 

c 

CO  TO 

IbO 

120 

•  RITE 

(3.130) 

130 

format 

(27X,2SH 

•  **PAR 

AHOLIC 

BOAITAIL*** 

//) 

c 

CO  TO 

IbU 

V. 

140 

•  RITE 

(3. ISO) 

ISO 

FORMAT 

(2I)X,23H 

***CON 

ICAL  BOATIAIL***  // 

) 

C 

ibO 

•ritf; 

(3.170) 

X2E. 

h2F.,  PETI)2E, 

EMNE, 

COBT,  PRILE 

170 

FORMAT 

(  ISX.bH 

X3F.a  ,( 

b.),7X 

,bH  P2Ea  ,F6. 

3.9X, 

I4H  HETA2E(DEG)a 

1 

F7,3,/, 

ISX.BH 

ERNE  a 

,F7.4,  4X,8H 

CUBT 

*  #  F  6  •  3  p 

2  7X,9H  Htt/Pt  «  P7.S.//) 

C 

IHO  kMlTtt  3,190)  NC*Sl,G«MMAt,<;CI,Xll.PII.BrTI)lI.bMi«l  I ,  bMS  1 1  ,  PH  ini , 

I  CAHH*b,GCb,XIL.H|F,HrTUie,eMNIR,EHSIE,PHlUlC 

190  FORMAr(IOX,4tH  * « • «Thu*S1 HKAH  BASF  PRESSURE  PMOCRAH****  ,SX, 

1  lOH  PHOb,  NO,  14,/,4X,31H  «*****1NPUT  OAT A •*•••*,/ , 

3  28X,22H  444INTERNAL  STREAM***,  //, 

3  I3X,9H  GANMAIa  FS.3,  3X.16H  CAS  CONSTANT  a  P7.2,I1H  LB-PT/LB-R, 

4  /  ,13X,6H  Xlla  F6.3,7X,6H  Rlla  rN.3,9X,t4N  BETAIKOCOa  FT. 3,/, 

5  ISX.tH  EMNir  arT.4,4X.SH  EMSll  BrT.4,l<X, lOH  Pll/PUt  aPT.S, 
b  /,2RX,22H  •**E(TEHNAL  STREAM***.  /, 

7  1SX,9H  CAMMAEa  rS.3,  9X.tbH  GAS  CONSTANT  a  F7.2.11H  LH-FT/LR-R. 
a  /  .ISX.bH  XlEa  Fb.3.7X.bH  Rita  F6.),9X,14H  BErAtC(OEG)a  FT. 3,/, 
9  ISX.BH  ERNIE  aF7,4,4X,8H  EMSIE  ar7.4,bX.llH  PIE/PUIE  arT.S/) 

C 

NRITC  (3.200)  PRIIE.  TROFl.  RLORO.  ENGRO.BLMTE.BLMTl 

200  FORMAT(2tX,3bH  ***«**bASE  PRESSURE  CASE  DATA******.  /, 

1  ISX.llH  PlI/PE  a  F9.4.t7X.l|N  TOE/TOI  a  FB.S,/, 

2  tSX,  9H  PEDRO  a  CI2.S.  tbX,  9H  ENGRO  a  FI2.S./ 

3  ISX.  9H  BLMTE  a  ei2.b.  IbX,  9H  BLNTI  a  E12.5./) 

C 

NHITE  (3,210)  HECUMP 

210  FORMAT!  liX.  12H  ***REC0mPRF,SS10N  COEFFICIENT  a  F5.),  3H***,  f, 

1  ISX.SIH  *•♦•*»*•«*«*•*♦•••*♦**•****************♦*•**•*♦**•  /) 

c 

HF.TURN 

STOP 

END 


SUBROUTINE  ACPBS(GAMHA,EMSl.PRATlD.XCO,Rri>.BETAO,REMT,NCALC.NPTS, 
1  NPMlNT.NFkON.NPPTS.PPTS.NSHAPC) 

M***AXlSyHNCTRIC  CONSTANT  PRESSURE  PUUNOART  SURPROCRAN  (ACPBS). 

internal  flow  (NFLUaal)  •••  UNIFORM  UR  CONICAL  SUPERSONIC  FLO*. 

CALCULATIONS  ARC  FOR  THE  *LOHER>HALt* 
OF  THE  FLOW  FIELU. 

EXTERNAL  FLO*  (NFLOHa2)  —  INITIALLY  IInTFORM  SUPERSONIC  FLOP. 

CALCULATIONS  ARE  FOR  THE  •UPPER-HALF* 
OF  THE  FLON  Fir.LO. 

NOTE  —  INPUT  ANO  OUTPUT  DATA  ARE  FOR  THE  *UPPF.R-HALF*  OF  FLO* 

FIELD.  THE  adjustment  OF  THESE  DATA  FOR  THE  CALCULATIONS 

IS  nape  internally. 

SUMPNOGRAN  REUUIRES-— OUTPUT.pMSHR.ilFl.or.CNFLOC.FPS,  APS.CPms. 

«CI>ATA,im'Th|'Y,T>  ST. 


non  nnn  n  n  n  rt  n  rt  n  n  n  r>  n  n  n  n  r>  n  n  r>  n  r>  n  r>  n  n  n  n  n  n  n  n  nn  n  rt  <->  rt 


»»»V*HI  AhLF.S»*» 


UAMti*  «  KATIO  Ur  THt  SPKClfiC  HFATS. 

F.xSi  •  INITIAL  **«CH  SIAM  AT  POINT  1. 

PRATIUa  tXPANSinN  PPtSSUSE  PATIO  iP/Pni, 

ten  ■  LUNC.ITIIOINAL  COORUInatI  oHKPb  FXPAI^SIUN  IS  CKNTFRen, 

HCU  a  RAUIAL  CUOROINATE  NHCRE  EXPANSION  IS  CENTKRFU,  PUSTlVt, 
hftao  8  FLOW  angla:,  radians,  at  (tcn,kC'i)  ruR  iNFAtHNAL  flon,  pos. 
PLNT  a  LlHiriNG  VALUE  OF  THE  RADIUS  FOR  TERnInATInG  CALCULATIONS 
(MAX.  N  for  1NTA>NAL  FLO*  AND  MIN.  K  FUR  FXTERNAI.  FLO») 
NCALC  a  CURRENT  CALCULATION  NUMRAR 

a  I,  THE  INITIAL  CHARACTERISTIC  DATA  IS  CALCULATED. 

.UT.I,  initial  char.  DATA  TAAFN  FROM  ONE  OF  THE  STORED  ARRAYS 
NPTS  a  HO.  OF  POINTS  OR  INCREMENTS  ON  INITIAL  1 1 -CHAR ACTE R I  ST IC . 
NPHlNTa  -I  UP  0,  C.P.E«.  DATA  NOT  PRINTED. 

Al.  C.P.P.  0A1A  PRlNIEO. 

NELOH  a  I,  internal  FLUM. 

2,  EXTERNAL  Fl.Oa. 

NHPTS  a  NIIMRF.R  OF  BOUNOARY  POINTS  CALCULATED. 

HPTS  a  BUUNOARY  PUINT  DATA  ARRAY,  Hal, LIMIT. 

PMR,  CHARI,  CHARE  a  ARRAYS  FOR  METHOD  OF  CH AR AC lERIST ICS, 

«A«nUTPUT  DATA  (IN  ORDER)*** 

INPUT  DATA  TO  ACPHS 

pHATlOa  EXPANSION  PRESSURE  RATIO  (P/PO). 

EMN2  a  HACH  NUMBER  ALONG  BOUNDARY  AFTER  EXPANSION. 

ENS2  a  MACH  STAR  ALONG  BOUNOAHY  AFTER  EXPANSION. 

X  a  longitudinal  coordinate  or  BOUNDARY  POINT. 

R  a  RADIAL  COORDINATE  OF  BOUNDARY  POINT. 

THFTA  a  LOCAL  FLOa  ANCLE  AT  BOUNDARY  POINT  (IN  DEGREES). 


DIMENSION  PNR(100,S,2),  CHARKS.SD),  CHAPC(S.)O),  Pt(S),  P2(S), 
I  P)(S),  PPTSIS.IO),  SfGN<S) 

CONNON  PMB,  CHARI,  CHARE,  PI,  P2,  P3 
••*«INPUT  DATA,  SOME  OUTPUT  DATA,  AND  COLUMN  HEADINGS  ARE  PRINTED. 

CALL  0UTPUT(CAMMA,EMSI,PRATI0,BETA0,NPKINT,HPLUM) 

**«*SET  INPUT  DATA  EOP  THE  FLOa  FIELD  CALCULATIONS. 

GO  TO  (2,4),  NELO* 

2  NCa-PCO 
XCaXCO 
BETAa-BETAU 
CO  TO  b 
4  HCaPCO 
XCaXCO 
BETAaBCTAO 
b  CONTINUE 

****SeT  SIGNS  FOP  CONVERTING  OUTPUT  OATA  TO  THE  *UPPER-HALF* 

OF  TNC  FLU*  riELO. 

DO  JO  Mai, 5 
CO  TO  (10,20),  NFLOH 
10  SlCN(H)a(-| .0|**(M,i) 

CO  TO  JO 
20  SlGN(M)al.n 
JO  CONTINUE 

***«tTHE  NAXIMUM  NUMNEN  of  family  I  CHAP ACTEM I  ST  ICS  FUR  .RICH 

CALCULATIONS  APE  m«0E  IS  SPECIFIFU  HEHL  (MAX.  LIMIT  IS  JO). 

LIMiTaJO 

C*«**«THE  iNiriAL  ll-CHAR.  IS  NUa  SUHCIVIDFl)  A'.U  1  HE  IVIllAL  CUAR. 
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oon  o  or*.  on 


hat*  CAbClII.AIf  I)  (  AAA.  »  II.  I1»  I»  (  i'l- •■*•  >1T.S  »  2'i). 

IF(f'CAl.C-H  ^n.bO.no 
sn  CU  TO  (»<n.<A>l),  itrLU* 

•  ••••rPR  iNfKHntL  FLU*  flKI.Ii. 

60  1F(AbS  (bF.1A)-l.ue«4)  7U,7U.M0 
•♦•••FOR  UNIFORM  FLO*. 

7U  CALL  UFL(lC((;AMM»,r:*SI  .XC.RC.Nt'TS, CHARI  .NFI.l)*) 

CD  TO  110 

»****FnR  conical  FLOfc. 

RO  CALL  CNFLOCICAKNA.tMSl.hHA.XC.PC.HPTSl 
CO  TO  110 

*****FUK  LXTLRRAL  FLOA  FIFLO. 

40  1F(NSHAPF|  96, 06, 100 

•  ••**FUR  uniform  FXTF.RNAL  FLO*  N11Hr>UT  A  bOATTAlL. 

46  CALL  UFLUC((:AMMA,L«SI,XC.RC,LlM11-l,rHARF,4FL0*) 

NPTS«LIR1T 
CO  TO  1 1 n 

•*»*FUR  ONIFURN  EXTERNAL  FLO*  RITH  A  BOATTAIL. 

100  LlRlTaNprS 

•«*«THF  PRANOTL'NEYLR  EXPANSION  AT  (XC.RC)  IS  NO*  SUHOIVIDFD. 
no  CALL  PmSRR(GAHMA,FNS1  ,PRAT10,l>F.1A,XC.RC,K) 

•*««RI  IS  NUHHF.R  of  FARILV  11  CHAR.  FOR  SUHDIVlUCO  EXPANSION. 

K1  a  K  A  1 

«*«*STORAGC  or  INITIAL  HOUNOART  POINT  DATA. 

NHPTSal 
00  120  Ral,4 

120  nPTS(M,l)aSIGN(R)«PNR(Rt  ,M.I) 

•*«ATRE  INITIAL  POUNOART  POINT  DATA  IS  PRINTED. 

CALL  ouraurd.NPNiNT.HpTSi 

A**«THt  FLO*  FIELD  CALCULATIONS  ARE  NON  HADE  ALONG  FARILY  1  CHARS. 
starting  from  THE  INPUT  POINTS  UN  THE  SURDIVIDED  INITIAL 
FAMILY  II  characteristics  TO  THE  bUHNDARY.  THIS  SEOUENCE  IS 
NOT  APPLICABLE  FOR  THE  FIRST  AND  SURSEOUENT  AXIS  POINTS. 

DO  ISO  NaT.MPTS 

AAAALOAO  INITIAL  FAMILY  II  CHARACTERISTIC  DATA. 

DU  ISO  Mb1,4 
GO  TO  (130,140),  NFLOM 
130  PNB(1,H,2)8CHARI(M,N) 

GO  TO  ISO 

140  PMB(I,M,2)8CHARE(M,N) 

ISO  CONTINUE 

*««*CALCULAT1UNS  ARC  FOR  THE  CURRENT  N^TH  POINT  UN  THE  INITIAL 
FAMILY  II  CHARACTERISTIC. 

no  160  L«I,K 

C(**«*CALCULATIUNS  ARE  FON  THt  CURRENT  L*TH  EXPANSION  INCREMENT. 
C«A***LOAO  DATA/  FIELD  POINT  CALCULATION/  STORE  DATA. 

CALL  MCOATA( I ,L,L«I .L1,KPTS) 

CALL  FPSIGAMMA,  PI,  P2,  P3«  NERROR) 

IF(NCRRUR)  2T0,|S4,lb4 
ISA  CALL  MCOATA(2,L1.L2,La|.KPTS) 

160  CONTINUE 

C*****ALL  FIELD  POINTS  ON  N*TH  FAMILY  1  CHAR.  HAYE  BEEN  CALCULATED. 
C«****LOAn  DATA/  boundary  POINT  CALCULATION/  STORE  DATA. 

CALL  MrOATA(4,RAl,6«l,Ll,RPTS7 
CALL  CPHSCCAMMA,  pi,  P3,  P3,  NERRUR) 

IF(NCRROR)  270,164,164 
164  call  MC0ATA(2,LI «L2,F*2,KPTS) 

NBPTSaNBPrSAl 
DO  170  Mat, 4 

170  RPTS(M,N|asiGN(H)*Pl(N) 


onn  nnnoo 


C»»*»»CH*RACfF»lSTJC5  0*f*  SMlfT. 

CAM,  MCOATA( 3,LI .L2,L1,Ka2) 

C***«*THf.  CURRtRT  POUROAR*  PUINT  RATA  IS  HUR  PRI-iTCD. 

CALL  UOrhOY (R, SPRINT, hPTS) 

CALL  TESMMLNT.nSImT  ,nh,0».N,«PTSI 
GO  TO  (IRQ, 260),  NSTNT 

(;y****Al)V«NCC  INLCX  FOR  NKXl  INPUT  POINT  ON  INITIAL  CH  AN  ACTFR I  ST  I C . 
180  KS|(>1 

GO  TO  (190,260),  NFLON 

*««*ATH1S  SbUUCNCC  APPLIES  ONLT  TO  THE  INTERNAL  FLON  ONKHF  THE  AXIS 
POINTS  ARE  CONSIUERCO. 

**««*THE  NUNHEH  OF  POINTS  TO  RE  C AI.CUL ATEU  ALONG  EACH  FANILV  1  CHAR. 
IS  NON  constant  ANU  GIVEN  RT  R|. 

190  KIsKAl 

KPTSsKlAl 

NsNPTS 

««*««TMF.  ELENCNTS  in  THE  N-TH  COLUMN  OF  THF.  PMH  ARRAY  AHE  SHIFTED 
DU«N  ONE  ROn  to  SET'UP  THE  CALCULATION  SEOUENCE. 

DO  210  Lal,Kl 
LI  a  Al-L«l 
DO  200  M«l,4 

20R  PNB(LlTl«N,l)apaB(LI,N,l) 

210  CONTINUE 

C««*««THF  CALCULATIONS  ARE  NOW  HADE  ALONG  THE  (N«1)-TH  FAMILY  I  CHAR. 
220  NaNyl 

C«**««LOAO  DATA/  AXIS  POINT  CALCULATION/  STORE  DATA. 

CALL  MC0ATA(l,|,2.L).kPTS) 

CALL  APS  (GAMMA,  P2,  P),  NENRUR) 

IF(NFRROR)  220,224,224 
224  CALL  MCUATA(2,L1.L2,1,NPTS) 

C««***CALCULAT1UN  of  HF.MAINOER  OF  FIELD  POINTS  ON  N-TH  FAMILY  I  CHAR. 
00  210  La2,Kl 

C**«««LOAO  DATA/  FIELD  POINT  CALCULATION/  STORE  DATA. 

CALL  MCOATAd.L-l.LAltLl.KPTS) 

CALL  FPSIGAmma,  PI,  P2,  PI,  NERROR) 

IF(NCRROH)  370,228,228 
238  CALL  hC0ATA(3,LI,L2,L,RPTS) 

210  CONTINUE 

C***««LOAO  DATA/  boundary  POINT  CALCULATION/  STORE  DATA. 

CALL  HC0ATA(I,KI,A1«1,L1.KPTS) 

CALL  CPBS(GAMMA,  pi,  P3,  PI,  NERROR) 

IF(NEHRON)  270,2)4,214 
214  CALL  NC0ATA(2,L1,L2,R1«1.RPTS) 

NBPTSaNBPTSTl 
DO  340  Mat, 4 

340  PPTS(N,N)asiGN(M)«Pl(N) 

C**«««CHARACTY.R1STICS  DATA  SHIFT. 

CALL  NCDATA( ),LI ,L2.L1,APTS) 

C*****THC  current  boundary  POINT  DATA  IS  PRINTFP. 

CALL  OUTBUY(N,npr1nt.BPTS) 

CALL  TESKRLNT.NSTmT.NFLOh.n.bptS) 

GO  TO  (390,360),  NSTaT 

C**«»«COMPARISON  alTH  LIMITING  NUMHFR  UF  FLOM  FIELD  CALCULA T J HNS. 

390  IF(N-LINIT)  220,2»0,26U 
C***««IF  NECATIVF,  continue  CALCULATIONS. 

C*«**aiF  ZERO  OR  PUSTlVt,  RETURN  To  masTFR. 

360  CONTINUE 
270  RETURN 
END 

SUHROUriNE  CRnS5(CAMM»I ,mPT1 ,LI“ITI ,GAMM»r ,HPTE,LImI Tl ,N|C,'<tr. 
I  NSTUP,  I  J«LI  .T.IMIf  .PRSMiiy  .NURliri 


B-15 


•  ••••THIS  SUBHUUTI'<t  C»LtUI*TfS  THE  I **»  I PHINf  'If  THE 
SUPfHSUNlC  I«iCEKN«L  (1)  AMO  LXtIKHAI,  CF)  STREAMS. 

SUhRUUTINE  RCUUI  RES-* -t'RSHA.  SLIP. 

•  **VARtAtlLES*** 


CAMmAI 

MPTI 

LIMITI 

CANMAE 

BPTK 

MHITE 

MIC 

NEC 

NSTOP 


TJMLl 

TJMLE 

ri.ra  « 

PKSHUK 

NPRlNf 


■  RATIO  Uf  1HE  SPFCIRIC  HEATS  mP  IriEJ  INTERmAL  STREAM. 

A  Internal  stream  RouNUARir  hata. 

■  NUMHCR  UE‘  INTERNAL  STREAM  HOlINtiART  PUINTS. 

*  RATIO  or  THE  SPECIPIC  HEATS  I  OH  THE  EXTERNAL  STREAM. 

■  EXTERNAL  SIRIAM  HtlllNOARY  BATA. 

>  NIIHBCH  UF  EXTERNAL  STREAM  HOIINDART  POINTS. 

*  LOCATION  NO.  Of  INTERNAL  STREAM  IMPINGEMENT  POINT. 

A  LOCATION  NO.  OE  EXTERNAL  STREAM  IMPINGEMENT  POINT. 

A  0,  SOLUTION  rOlIND. 

A  I  ,  NO  impingement. 

A  3,  NO  SHUCA  solution. 

■  3,  impingement  hefope  separation. 

A  internal  IURROLENT  jet  mixing  LENGTH. 

■  EXTERNAL  TUHHULENT  JET  MIXING  I.ENCTH. 

ABSCISSA  FUR  AXISTMMETRIC  CORRECTION  TEHMS  HSEJO  IN  UNERA 
angular  RFATACHMENT  CRITERION 

■  STATIC  PRESS.  RATIO  IRISE)  ACROSS  JHLIJUE  SHOCK  SYSTEM. 

■  SEE  SUBROUTINE  AINOIITA. 


HPTI(M.N)  AND  HPTE(M.N)  ARE  ROONOAMV  POINT  DATA  ARHATS  nHCRC 
Mat. 4  ANO  INOICATES  VARIARLC  AS  IN  PHR  ARRAY. 

NaI. LIMITI  UR  LIHITE  INOICATES  THE  HUUNUARY  POINT. 


OIMCNSION  XI( JO),RI(30),XE(}0)>HF( IU).RPT1(S»10),HP1C(S, 30) 
common  /AnglES/THETAI .THETAS.THFTAS.FI ,F3 
EMNMSF(eMS,CAMMA)sSURT(((a.0A(EMSAA2))/(CAMMAEl. >)/ 

I  (I .0«((GAMma«1 .0)/(CAMMAt|.u) )a(Emsaa2)) ) 

C«aaaaluA01nC  or  CONSTANT>PRESSURt  BOOROARY  POINT  DATA. 

00  10  N81.L1N1TI 

Xl(N)saPll(l,N) 

10  RUN)  a  HFT1(2,NI 
00  20  Na|, LIHITE 
XE(N|  a  BPTEd.N) 

20  RE(N)  a  BPTE(2.N) 

caaaaaset  initial  values. 

NSTOPal 

PPSHOKaO.O 

NtM«XaLlNITl-l 

NEMAXaLINITC'l 

CAAAaachcCK  for  IMRIHGEHENT  upstream  HF  THF  separation  Pi)|NTS. 
caaaaaeor  the  internal  strf:am. 

SFaO.O 

NEal 

00  30  Nlal.NIMAX 

SI  a  (Rl(Nt*l|  .  M|(Nin/(Xl(NI«I>  •  XIINII) 

IP( ABSCSC-SI )  .LT.  I.UE'OS)  CO  TO  30 

XIHP  a  (HI(NI)  •  RE(NL)  «  SFAAKNE)  •  SI  •  X  I  ( N I  ))  /  ( Sb  •  SI) 
IFTTXIMP.CE.XI (NI ) ).ANI>.CXImp.le.XI(R1*I ) I.anU. 

I  (XINP.LE.XC(HE)))  CO  TO  SO 
30  CONTINUE 

caaaaaeor  the  external  STHEAm. 

SlaO.O 

NIal 

DO  40  Nial.NEMAl 
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St  a  (HtiNr»i)  •  (N^ )  j/(  xt  1 )  -  xr(Nt)) 

IK(*bS(St«.SI)  ,LI.  l.<.K-Ob)  i;0  T(‘  40 

XIMp  a  (Rl(Nn  -  HK(Nn  *  St*XI(St)  •  SI  XXI  (141  )  )/(St  •  sn 

1F(  (Xl«*P.ar..XF(NC|).**'l>.  (  X1»P.I>  .XE(NI  ♦(  )  )  .4N0. 

I  (Xlap.l.t.Xl  (N1  n  )  GU  TU  70 
40  CUNTINue 
GO  TU  100 

Ctt***JK  IxplNGEurXI  UCCUBS, 

so  RJMP  a  (SL«St*(XF(NX.)«Xl(Nt))  ♦  StaPKNt)  »  S I  •*<£  ( N  t ) )  /  ( SK-S I  I 
C 

■RlTt  (3.b01  XIHP.RIPP 

bO  FUHMXTt  ISX,  4tlH  »•••♦*  I  •‘PlbCtPEHT  fIF  TMF  INTEPRtb  STKIXM  UCCUHS  / 

1  7IX.  47H  MEFURE  StPXRXTlUN  UF  THE  E:xTtRNAL  SThEAH* *  * • • •  ,  /, 

2  IbX,  27H  I*<PtMG7:MFAr  UCCUPS  AT  X  «  FI0.6,  SX,  9H  AHO  R  a  F10,6  /) 
GU  TU  90 

C 

TU  RIPP  a  (StbSl * (Xt(NE)«X I (Nl  )  )  *  SL*R1(NI)  •  S I • RE ( NE  )  ) / (SE-Sl ) 

C 

■HITE  (J.NOI  X1PP,RINP 

ao  FURMAT(  ISX,  4RH  1 mP t RGEMbPT  OF  THE  EXTERNAL  STREAM  OCCURS  / 

1  21X,  47H  HEFORT.  StPAKATlUN  OF  THE  INTERMAL  STRE  A 'X*  •  •  *  •  •  ,  /, 

3  IbX,  37H  impingement  OCCURS  AT  X  >  FlO.b.  SX,  9H  AND  R  ■  FlO.b  /) 
C 

90  NSTOPal 
GO  TU  210 

C***««CALCULAT1UN  OF  CUNSTAN f'pRESSURC  hOUNDAPlES  IMPINGEMENT  POINT. 

100  UU  120  Nlai.NlMAX 

SI  a  (RKNIbl)  .  Rl(Nl))/(Xl(Nl»n  •  ll(Nl)) 

DO  110  NEai.NEMAX 

SE  a  (RE(NE«n  •  RF  (NCn/(XC(Nt4l)  •  XC(NE)) 

IF(A(|S(SC>S1I  .LT.  l.OF'OS)  GU  TO  110 

X1«P  a  (RKNI)  •  NC(NF.)  *  SCaXFINE)  •  SI  «X  I  (N I  )  ) /(SE  •  SI) 
IF((XlMP.GC.Xl(Mn)  .AND.  (XIMP.Lf  .XI(NT«1))  .AND. 

1  (XlNp.CC.XCCNC))  .AND.  (XIMP.bE.XElNExl)))  GO  TO  140 
110  CONTINUE 
130  CUNTINUE 

C«44««F0R  NO  impingement  of  THE  STREAMS. 

WHITE  (1,110) 

110  FlIRWATdOX,  4tH  »  •  •  1  wp  1 NCCMENT  DOES  NOT  OCCUR  wITHIN  THE  ,/, 

1  I9X,  44H  RANGE  UF  CONST ANT-PRESSURF  BUMNOART  DATA***  /) 
NSTUPa3 
GU  TU  210 

C*****FUR  IMPINGEMENT  OF  THE  STREAMS. 

140  RIMp  a  (SC*Sl*(XE(Ne)>XI(Nl))  *  SE*RI(NI)  •  S I *HE ( NE ) ) / ( SE-Sl ) 
NICaNltl 
NFCaNE*! 

C**«**lNTERPObATION  FOR  THE  FLOW  VARIAHLES  AT  THE  IMPINGEMENT  POINT. 

00  ISO  M«l>4 

BPTI(M,nK)  a  BPTKM.NlC'l)  ♦  ((XINP  -  Xl(NIC-l))/ 

1  (XI(NIC)  -  XI(NIC-I) ))*(MPT)(M,NIC)  -  BPTKM, NIC-1)) 

ISO  BPTF.(H,nIC)  a  ilPTE(N,NEC-t)  ♦  ((XIMP  -  Xt(NEC-l))/ 

1  (XC(NEC)  -  XE(NtC-l)))*(MPTt(N,NtC)  -  BPTEI M , NEC- I ) ) 

C*****STORC  COUROINATES  UF  THE  IMPINGEMENT  POINT. 

BPTKI.NIC)  a  XIMP 
BPTI(2,HIC)  a  HIMP 
BPTLd.NEOa  XIMP 
BPTC(3,NLC)  a  RIMp 
THETAIa  BPTI(4,N|C) 

THFTA2a  PPrL(4,NEC) 

C**CAI,CIILATION  UF  THE  muinG  LENGTHS  A*  l>  OF  F  FUN  AXIST*.  iimFRA  CRIT. 
TJMLlaO.O 
Flao.O 

UU  IbO  Na2,NIC 


(>S>SMHT(  (M»>TI  (  I  1  (  »  ,u-1  )  )»»;?+(HPrH2,h)-i'Hril  2,'J-n  )**i) 

TJHLI«rjMLI*US 

160  H*Fl»(BPTI(2,N)*t»Pri(2,N-l))»r».S»l)S 

F  |sF|/(TJhl,l*KlMP) 

F2«0.0 
TJMLt'>0.0 
UO  no  Hs2,HKC 

OSaSUHri  (HFTF.(  I  ,  N  ) -BPl  t  (  I  ,  N- 1  )  )  •  •  2«  0»PTt  (  2  » ‘O  -  HPTE  (  2  ,  M- 1)  »  •  *  2  ) 
TJFLFBTJxLCtOS 

no  F2«F  2'»(BPTC(2,N)*hPlfc(2»N-i  )  )*0.!>*US 
F2»F2/tTJ"bt»R|aP) 

C*****OlirPUT  IM^tiNGCxeNT  POINT  DATA. 

tMNi  s  EMNMsr chpric j,nici,oamfai) 

THF.TOi  a  bl  .  BP  r  I  (  4  ,*i  IC  ) 

CHNt  a  ENhasr(HPre(i,aeCI.GAMMke) 

IHiTDE  a  57.20bn9S«BPTE(4,NtC) 

IF(NPHINT.I.T.U)  tiO  TO  200 
C 

aRlTt  (3,100) 

IRO  FUMMAT(  IHl  J 

c 

aPlTE  (3.140)  X  IMP,  R  lap,tMal  .THETOI  ,TJOt.l , 

1  XlRP,Rlap,FaoF,THFTl)E,TJRL« 

190  FORMAT!//, I•X,42H*««AT  lOTCRNAL  STREAM  IMPINCCMEMT  POINT***  ,//, 

1  5X.  9H  X  a  Flu. 6,  9X,  3H  R  a  FI0.6,  SX,  12H  MACH  HO.  a  F10.6,//, 

2  5X,  19H  THLTAIObG.)  a  F10.6,  5X.  ITH  MIXING  LENGTH  a  F10.6,//, 

3  lOX,  4)H  •••AT  EXTERNAL  STREAM  ImPinCEMFNT  POINT***  ,//, 

4  9X,  5H  X  a  Fto.6.  bX ,  SH  R  a  F10.6,  SX,  I2H  MACH  NO.  a  FlO.b,//, 

5  SX,  ISH  THLTAIOEG.)  a  F10.6,  SX,  ITH  MIXING  LENGTH  a  F10.6,//) 

C 

C**«**CALCULAT1UN  OF  THE  RECOMpRESStON  SHOCK  SYSTE*. 

C*****CALCULAT1UN  OF  THE  SLIPLINC  ANCLE. 

200  CALL  SL1P(«PTI(3,N1C),SPT1(4,N1C),GAmmAI, 

1  BPTE(3,NCC),ePTE(4,NFC),GANMAC, 

2  THETAS, NSTOP) 

C***«*OUES  THE  SOLUTION  FOR  THE  SLTPl.INE  ANGLE  (XIST. 

GO  TO  (210,230,230),  NSTOP 

C*«***CALCULAT1UN  OF  THE  STATIC  PRESSURE  RATIO  ACROSS  THE  SHUCK  SYSTEM. 
C  (NOTE  PHSHUKtapRSHOKCaPRSHOK.) 

C 

210  DELTA!  a  (BPTI(4,N1C1  •  THETAS) 

PHSHOK  a  PRSHK(BPTI(l,MlC),nFLTAl,GAMMRI) 

THFTAS  a  ST . 29STT9S*THCTAS 
IF(NPRINT.LT.O)  go  to  230 
C*****0UTPUT  or  SHOCK  SYSTEM  UATA. 

C 

NRITE  (3,220)  THETAS, PRSHOK 

220  FORMAT! ISX,  4«H  •••OBLIQUE  SHOCK  SYSTEM  AT  IMPINCEMCNT  POINT***//, 

1  SX,  23M  SLIPLINE  •NGLC(OCG.)  a  »|R.b, 

2  SX,  24M  STATIC  PHESSUHE  RATIO  a  FI0.6,//) 

C 

2)0  RETURN 
END 

SUBROUTINE  TJM1X(GAMNAI ,CCI .EMSI.TRBOl.TJMLl, 

1  CAMha3,CC2,FMS7,TRR02,TJML2, 

2  RN|  ,EMSN1  .HFTANI  .RIMP.PHSHilK, 

3  PR021  ,TPU2I  ,Rr:COMP,B|,t>R.ENCN,BLM  ri  ,HL"T2.  ITC) 

C 

C«a***tH19  SUBR0UT1MF.  CALCULA1ES  THF  UIMFNSllHiLFSS  BLEED  An|> 

C  energy  RATIOS  FOR 'the  laO'SHRFAM  T HTFRACT I'J.<  PMUHLEM. 

C 

C  SUBROUTINE  HC0UIHE5-— TEGRAL. 

c 


♦  ••VAHI  »l*LKS»*» 


FOK  F.1TH»H  S1He*«  1  UR  i 

gamma  a  KATIU  OF  SPECIFIC  hlATS. 

GC  ■  <i*S  C<iriSTAhT---(LbF-FT/l,B“-R>. 

t«S  »  maCM  Sr»R  A1  XMHlNGF.MtNT  »-UlMT. 

TH^TA  S  fl.O>  ANGUL  AX  IMPIUGFMC'^T  po|MT  (  I M  RAUIANS). 

TRHO  *  BASK  TO  FRtt-SIKtAM  STAGrATJUN  TKMP^MATIIM^  HAHU, 

TJML  ■  TURMUGKAT  JfcT  MIXING  IKAGTM, 

K  «  AHSCISSA  FUR  A X  I S < mmKT P IT  COPR EtT I DN  TtH^S  IN  UKKRA  AN' 

angular  RtAlTAC'MtMT  fHlTFRlClN 

RN|  a  NUZZLF  EXIT  PAUIUS  OF  STREAM  I  (fNrF.RNAL). 

EMSNl  a  nozzle  LX|T  maCH  STAR  UK  STREAM 

RETANI  a  nozzle  EXIT  Kl.Oa  ANCLE  AT  RNI  (IN  HAUIAtS). 

RIMP  a  RAUIAL  OUROINATF  OF  I MPI MGFMFKT  POINf. 

PKSHUA  a  STATIC  PRESSURE  RATIO  (RISE)  OF  OHLIUIJE  SHUCK  STSTEM. 

PPOIl  a  STAGNATIUN  pressure  RATIO,  PU2/POI. 

TROZI  a  RATIO  UF  STAGNATION  TKMPERATIIHF.S  UF  THE  T*U  STREAMS. 
RECUMP  a  RFCOMPRESSION  CULFFICIENT. 

RLOR  a  MASS  SLEF.D  RATIO  REKERtNCED  TO  FLO«  UF  STREAM  I, 

(G  HLCEO)/(C  NOZELLl). 

ENGP  a  ENERCr  BLCtU  RATIO,  (UMECAR)/((C  NOZZLE I ) aCP 1 *TUi ) . 

■HERE  OMEGAB  is  NEFERFNCEU  TU  TaO, 

CUHMON/ ANGLES/ THETA  1 ,THFT A3, THfTAS,f 1,F2 

COmmON/CnJ/SIGMFI ,S1CmF3,NNCHP,RSIGma,NBLEE0,XPAGE,NSH1FT,THRSLP, 

1  RECMPl ,RECMP3 

CH3MSF(EMS,GAmmA|  a  ( (GAMMA>| .0)/ (GAMm Ak I .U ) ) « ( EMS* *2 ) 
EHNMSF(EMS,GAMMA}asuRT  ( ( ( 2.0* (EMS«*3 ) )/ (CAMNAkI ,0 > ) / 

1  (1,(;«((GAMMA«1.U)/(GAMMAK|.0) )*(EMS**2))  ) 

EHSMRFCtMN, GAMMA)  8EMN •S0RT( O.S* (CAMM A K 1 . ) / ( I , tO , S* ( CAMM A> I . ) • 
1EMN**2)) 

CMSPRF(PR,GAMMA  laSuRT  (  (  (CAMmAkI  .D/CCAmma-I  .0)  )a 
1  (1.0-PR«*((CAMMA-l.n)/GAMMA))) 

HTFLMS  (EMS,GAMMA)aSURT  ( 2 .0*CAMM A/ (C AMM AK 1 . 0 ) ) a 
t  (tMS/(l.O-((CAMMA-I.«)/lCAMMA»| ,0) »•(tMS••2)) ) 

PHMSF(eMS,GAMMA |8( I ,0>( (GAMMA«t.O)/(GAHMAtl .0) )«EMS*»3)** 

J  (CAmmA/IGAHhA-I.U)) 

PHIDF(CNR,TRBO)  a  CNR* (O.S«CNRa( I .O'TRRU)  ♦ 

I  SORT  ((CNR*a2)*((U.S*(I.O«TRBO))**2)KTRBO)) 

m*««pranotl>mcyf.r>fiimction 

0MCGAF(EMH,GAMM A )8SURT((CAMMAmI.)/(GAMMA>|.)) FA  TAN (SURTf (GAMMA-1.) 
l/(GAMMA«l.)*(EMN«*2>l.))}>ACOS(l./(MN) 

••••■••page  CRITERION, APPRUXIMATFO  BT  KORST 
XKAPROIEMNJ  a  0.02 J J J«ERN*0.6B)3) 

••«*«0R1GIRAL  PAGE  CRITERION 

XRAPPA(PMIO)  a  O.S*(l.*CUS().l4l!>9*(PHID-n.OI))) 

•«*a*ONLRA  ANGULAR  CRITERION 

RKATTACHNCNT  ANGLE  FOR  COaO  IN  PLANE  FLOH 
P810(ENN)  a  0.O442HISK0.ITI94*KMN-O.ntT9S?*F.HN*A2 

P810(FMM,PHIJ,GAMMA)a0MtC4F(FMN,CAHMA )-0MFG4F(E"N*SURT( 1 .-PHI J»»2) 
1 .GANNA) 

CHANCE  Of  P810  «ITH  CU 

PSIPH0(F,NM,PHIJ,GAMMA,TRH0)aSORTItMN»»2*( I .-PHlJ**2/( IRHUkI I .-TRBU 
I )*PHIJ))-I. )/(l .♦'l.b*(GAMMA-l,)AKMN*»2*(l .-PHIJ**2/(TRHU*( l.-TRMII 
2*PMlJ)))*PHIJ/S0RT((l,-PMia)*(TI.RO»l*HIJ»*(  I  .-PM1.I»»2/(TRH0K(1.- 
ITRRU)*PMI  J  J  )  )A(THH0K(  )  .-TRBIi)*FMIJ/2.  )/(TRM'j*l  I  . -TRHU  )  •PH)  J  I  •  •  1  .  b 
PMICO(PHlJ,KTAJ,CS«0,TNHiJ)a(lwM|.«(  |  .- fh  nii  )  •PHI  )-PH  IJ  •  •  /  •CSUn  )  /  (  t  .  • 


2CSui)  •O.SM  I 

C  FXCrOB  FUH  Si-CHSU  (iHlitK  (>t  f'MIi'  i.ITr<  CO 

P2COS>  (PHlJ.tTXJ.CSUf,  IPHU)*(  tlHHO*  l-HI.J»»2»CSijr>)* 

10,5t>4J<»*F.XPt-tT»J**2)*2,«rrMHO*«  »  ,-TPHii)*PHlJ-PHI.)*»2»CS‘)D)*f  l*J» 
2PHIJJ/PHIJ»»2/(! .-CSUD) 

C  ClIKKXCriON  Ft)H  AXlALtV  StMHKTPlC 

DPSl  *X  (  AKi:  )  s  11.114  M  24«0.06  IOh7*  (  AHG«*2*1  )  /  (  AWG*  *7*  1  ) 

APC(F)  «  tXP(r»b.4-4.4) 
c  coP»tcrio*<  run  ‘•om  hact  hpikip  shift 

CdRUSIEH**  J«2,'*OH(i-l  .OS4»tHN*O.I5bl  2*t*'K**2 
C»»»»»BSI.K()<S  CDRWKCIIUH  for  LIHITINC  AMKI.FS  I.IMI 

EPSIl-FIF«i»)«0.2<»li7b-0. 1  7  73JS*tMP*0.04  1  321  5»ENN*»7-y  ,00  Jl  1  Jy  I  * 

IFM  4«»1 
C 

FPS«I .t-i 

IF  inpchp.lt. 3)  emu  b 

PKUOlHl . 

THtrtSaTHETAS  /57.2S577P5 
PR003ct . 

C«**««UTARH  accounts  for  a  correction  of  RfeATTACHHENT  A'lCI.ES 
C  OUF  TO  THE  CONFLUENCE  PkObI.EM  OF  TWO  SHEAR  LAYERS 

OTARRaO. 

C«««A*bPStL  ACCOUNTS  FUR  ONEHA  LIMITING  ANGLES  IN  CUMHESPONDENCE  TO 
C  THOSE  or  KORST 

EPS1L1bEPSILF(F.mnmsF(EmS1  .GAHMAI  t) 
EPSIL2«CPSILF(EMNMSF(F:mS2,CAHha2)) 

ITaUSTal 

NTtPlal 

NJUHPaNRCMP«3 

DTa«0.U2 

OPSlaOPSlAX ( ARGIFI ) ) 
nPS2aOPSI  AX(ARG(F2n 
IF  (NRCMP.LT.b)  GOTO  7 
22  OTlaTHETAl«THETeS«OPSI-urARR 
DT2aTHF:TES<rHETA2*0PS24UTAR8 
IF  (OTl.GT.O.)  GOTO  24 
0TARH80Tl4nTARR>t.l>b 
GOTO  22 

24  IF  (DT2.GT.0.I  GOTU  2S 
0TARHa.(UT2>0TAPe)4l .(-6 
GOTO  22 

25  nCLTHlaOTI/ITHETAl-THETESI 
DrLTH2aoT2/(THCTES- theta?) 

GOTO  3 

b  IPAGElal 
1PAGE28I 
NTTPblal 
NTYPETal 

OXtaO.OS 

0X2*0.03 

IF  (NRCMP.NC.2)  GOTO  5 
TJNLlaTJMLl-XPAGE*TJML2 
TJNL2arjHL2*( l.-XPAGE) 

3  IF  (NRCMP>1)  1,2>2 

C**«**CALCULATION  OF  DISCRIMINATING  STHEAMEINC  VELOCITY  RATIOS 
C  BASED  ON  THE  RCCOMPRESSION  COEFFICIENT.  SINCE  THE  PRF.SSItRE  RATIO 

C  ACROSS  THE  OBLIQUE  SHUCF  SYSTEM  Is  EQUAL  FOR  SIHEAHS  I  Anl  2, 

C  THE  DISCRIMINATING  STREAMLINE  STAGNATION  PRESSURE  PATIO,  P/POl), 

C  IS  ALSU  THE  SAME. 


I  PROntal  1 .0/(RLCOMP*l'RSHli|l  A) 
PROOTaPHOUl 
RF.CMPlaRECOMP 


B-20 


oo  o  norio  nnon 


OtiTii  7 

•  •»*»c»i,rui,»i  i0'4  fiF  nisct^i**-  iiim:  si**f  axi.jji  vvuutiTr  w*tio  haskh  m* 

PACK'S  CKITtKin-J,  AlphUXlMATFO  »•»  tillPST,  irKPATl^V^  FUM  MOKE  EXACT 
EVALUAIIUN  or  FACE'S  Ch  i  TF.I.  I  I'EP  (»/Pul>)ll  IS  SlKFKKt><T  FU«  STKFA"  I 
ANO  STNKAM  2, 

2  ntLTHlsXKAPKU(KM(,«SF  (F,»'SI  .CAKMAl  )  ) 

DtLTH2  =  XAAPxn(t'<'»‘<SF  (  fc*(S2  ,CA>'m*2  )  ) 

J  PPOOJaPKSHX  (FXS!  , -(DELThI •(THE TAl -THF TKS ) ) .CAKMA I ) 

PHUn2*PRSHA  (E«S2, -(l)€LIH2» { IHf TA2-THF  f ES ) ) , CAX" A 2 ) 

RECxPlsPROOt /PRSHUK 
RFCPP^=PRO02/PRSHOR 
PROUlal ./PWOni 
PhO()2b|  ,/PRUU2 

•«**«FnR  STREAP  J. 

7  CS'JUlU=CR2''Sr<F>’SPRf  (  PRUO I  ,  CAMH  A 1  )  ,CAK«AI  ) 

CSUlJl  =  CR2RSF(F‘'S1,GAHRA1) 

C1=SIJRT  (CSuni) 

CNRI  *  SORT  (CSOOIO/CSUni ) 

PHlUl  -  PH1UF(CNR| .TRbul  ) 

•***«KUR  STREAM  2. 

CSQ(l2D  a  CH2MSr(EMSPRF(PR002,CAMMA2),GAMMA2) 

CSOn2  a  CR2MSr(EMS2,GAMMA2) 

C2aSURT  (CSa02) 

CNR2  a  SORT  ( CS0n2P/CSU02 ) 

PH1P2  a  Pt<IPF(CNR2.TRPU2) 

•  ataaiNTEHATION  FOR  MORE  EXACT  PAGL'-rRlTERlUN 

IF  (NRCmP.NC.2)  GUTU  4 
PELlaOFLIHI 
PKL2SDCL1H2 
OELTRlaXEAPPACPHIOt ) 
lteLTH2aXKAPPA(PHI02) 

VARlaDCLTHl^DfLI 

VAR2aUEtTH2«OEL2 

IF  (NTTPEl.NF.l)  S IGN 1  a V AR t / ABS ( V AR | } 

IF  (NTrPE2.hC.3)  SIGR2aVAM2/AnS(VAR2) 

IF  (NTYPEl.NE.J)  CALC  ITER  ( DECTHI , DX 1 , 1 .C>4 ,SI G4 I , V ARi , 0. 0 . EPS, 
lIPAGEl .NTYPEl ,XN1, (N| ,XPI,yPt .NSICAI.NSIGRI ) 

IF  (NTYPE2.NE.3)  CALL  ITER  (OrLTH2,OX2. I  .F>4 , SIGN2 , Y AR2 , 0 . 0 .EPS, 
UPACE2,NTiPF3,XN2,rN3,XP2,YP2.NSl(;A2,NSlGB2) 

IF  (PELTHl.LE.O.)  nELTHIa0.54F.PS 
IF  (OCLTH2.LE.O.)  UELTH2aO.S*EPS 
IF  (DELTHl.CT.I.)  UELTHIal. 

IP  (UELTH2.GT.1.)  UELTH2BI. 

IF  (NTyPll.EQ.3.ANO.NTYPC2.E0.3)  GOTO  4 
IF  (IPAGEI .GT.3n.OR, 1PAGE2.CT.20)  GOTO  49 
GOTO  3 

A««aaEVALUATION  OF  BLEEU  AND  ENERGY  RATIOS. 

4  SIGMAla|3.0*3.TSHaKMNMSr(EMSI ,CAMma| ) 

SIGN A3  al3.n43.7SR«FMNMSF(eFS3,GAnMA3) 

IP  INRCmP'6)  37,3«,3T 
39  EllUlaEPSILI/(l.-CSUUl) 

C11038CPSIL3/(t  .'-CS003) 

GUTU  311 
37  CllOlaO. 

E1103a0. 

21  CALL  TEGRALIPRIOl.CSOni .TRHUt ,EIIJl,FIIt>i , E 1 3 J I , E I 3U 1 , 7 TA J I , 

1  PMrjI.ElSJII 

CALL  TECRAL(PH|02,CSUD2,TRb02,FIl J3,F:lin2,EI3J2,EI3n2,ETAJ2, 

I  PHIJ3,C12J3I 

•  ««*aAPPLlCATIUM  PF  UNF.RA  ANGULAR  rRIlERION 

ACCOUNT  FOP  origin  SHIFT 


IK  (ITC.I.K..2)  Gliro  9 
IK  (NSMIKT.nC.O)  GUTO  400 

TJMLlaTJ»'GlKSlG«**l‘Hl,(<Tl/(tI2JM(l.0-CSi)Pl  ) ) /CUROS  ( t»tNHSF 
2(F.MS1  (CAMMAI  )  I 

TjaL2«TJ*‘l.2«SICMA2*HL»'T2/<tl2J2*<  I  .0-CSl)t)2  )  ) /CURDS!  EXiiMSK 
2(r><S2,CARRA2) ) 

400  IF(NSHIFT.Nt. I )  GOTO  9 

TJRLlaTJRLI  KSIGRAi  ARL^Tl/dl  1.11  •<  1  .-CSOni  )  ) 
TJMl.2«rj“L2*SIC><A2«HI,>‘T2/(KI  1J2»(  1  .-CSU02)  ) 

9  IF  (NKCRP.tT.l)  GOTO  8 
FHN|stMNRSF(FRSI,GAMM*l) 

K*‘N2«tRN*»SF«F*«.S2,GARKA21 
IF  (NRC"P.FO.»»)  GOTO  21 

PHICOIaRHICOlPHlJI ,tTAJ| .CSUni .IRBOI ) 
PH]C02=PHICO(PHrJ2,KIAJ2,CSUO2,TRBn2) 

1 1  OTlx-THtTEStTHKTAl-OPSI AX( ARGiri ) ) -PSI  (J  ( F" ’M  , PH IJ 1 , GAR M A 1 ) -UTARB 
0T2«  THETES-THK:  rA2-0PSIAAI  APC(K2))“P510(KRfi2,PMIJ2,GARMA2)*0TARB 
CU|a(OTI/PSlPHU(ERNl ,KH1JI ,c:amhA1 ,TRB01 ) / PH  ICO 1 -EPS  1  LI  ) 
C02«(UT2/PSIPHO(ERi42,PHIJ2,CARKA2>TRBn2)/PHtCU2-EPSIL2l 
EIlDtaEIIJlKC01/( l.-CSOUt) 

En02aFIIJ2»Ca2/(  I.-CSU02) 

EH1P|a-COI*rjRLl/SIGPAI 

FHT«2=-CQ2*TJRL2/SICRA2 

FACt«0. 

FAC2*0, 

IF  (HHCRP.E0.5J  FAC1«P2Cu5KIPM1J1,ETAJ1 ,CS0D1 ,TRB01 ) 

IK  (NRCHP.EU.S)  KAC2aP2COSK(PHIJ2.LTAJ2,CSOD2. rRB02) 
PHIOlaPHIJl«C0I«PHlC01«(l,-COl«FACl/2. ) 

PHI02«PH1J2kCU2»PHIC02*( 1 .-CQ24KAC2/2.) 

IF  (PHIOI.LE.l.E-3)  GO  TO  lUOO 

CALL  TECOSLCPHIDI ,CSOUI .TRBOI.FIIOI ,EI3DI ) 

GO  TO  loot 
1000  EllOlaO.O 
CUUlaO.O 

toot  IK  (PH102.LF.I.F-3)  GO  TO  1002 

CALL  TECDSL(PHI02,CSQ03.TRBU2.E1I02,CI302) 

GO  TO  21 

1002  Ln02a0.0 

E1303a0.0 

21  PRRlapRRSr (L"SI ,GAMMA1 l/PHMSK ( t MSHNF ( ERN 1 (SORT ( 1 . -PH  1 01 4 • 2 ) > GARM A I 
I  1  ,GAM"AI 1 

PRP2*PRRSK(F.RS2.GARRA2  l/PRPSFdRSBNF  tEHR2»S0RT(  I  .-PHID2»*2)  ,GARRA2 
|),GAMnA2) 

IF  <NRCHP.LT.4)  GUTO  10 
VARaPRRl-PRR2 
IF  (VAR.EU.O.)  GOTO  10 
SlCHaVAR/ARSIVAR) 

IF  I  ITrUST.nE.  I  .  1  OT8-(OTARt-OTARBOI/I»AR-(4ARO)»VAH 

DTaSlGN«DT 

VAROSVAR 

UTARBUaOTARH 

CALL  ITER  (urARR,0T,l.t-4,SIGN,VAP,U.0,F.PS. ITaOSr.RTTPK.XN.YN.XP, 

I yP,NSlCA,NSIGB) 

IF  (HRCRP.EU.6.AR0, (PRRl .LT.PRR2.AN0.0ri.Lr.2.t-6.0R.PRRI.CT.RRR2. 
1AND.UT3.lt. 2. F-Rl)  GUTU  10 

IF  (NTIPE.Nt.I.ANH.llRUST.LT.l” 1  GOTO  (  11  ,  1 1 . 22  1 . NJliRP 
10  lHR8LPB(TriC'US*UTAHh)«S7.29b» 

RFCRPIaPKRl/PRSHUA 

HFCRP2aPRR2/PRSHUK 

■  PRBNlapRHSFCLPSl.CAMRAI )/PRRSK I K K 5HI .Ga R« A  I  1 

CUtFf  1»(  (  I  .04CU5  'IBKTANI  )  IFSlr.RA  1  )  ♦  I  R  I  "P/H  R  I  )  •  (  I  JRl.  1 /RN  1  )  •  (  PHMN  1  )  « 
1  SORT  (2.oaG4R"AI/IGARRAl-l  ,«  I  1*(  I  .n/»Tf  LR.S  (  E  RSr  I  ,  GAR- A  I  )  J 
CnEFF2a  (  TJHI,2/TJ«L  1  )  »  1  .Sir.-A  | /S  IGRA?  )  ».S dR  I  1  (  I  .  D/ TR 1  > 2 1  )  • 


r>  o  n  n  o  r>  o 


I  ((;*!>  ■'*7/0 A 1)  •  (GCl  /GC2  )•(((;*’<>'*  1  -1 .0  )/(G*'<'‘*/-l  .11)  J ) 

CD)  FF  J*(S1G'‘»I  /SlGf  *7)*(IJ«L2/TJH,1  )»S1)HT  (  (GC2/GCI  i»(lP(i2n  >• 
I  ((« (;*'‘«»2/G*hP*l)»((GAi'**»l-l.<i)/(<;»««<*2-1.0  )))••», S) 

Pl,t)R*-t0f.f  FI  •  (r  1  •  u.  1  im  -tl  I  Jl  )*CLiF  FF2»C2»(FilD2-t;itJ2)) 
f  NGP»-ctiFf  f  I  •  (Cl  •! 1 1  jui -rRKrn*Fi  iJi  )*rnt:i  f  i*C2* 

1  ( 1 1 J02- rRBU2*Ll I J2)  ) 

RtTUHN 
99  STIJP 


tNn 

SUPXnullNt  0IIT2'<<PRUF.  I'RHt  1  ,F‘NIiF  01  ,TPnilF,l  RHUl  ,  THI)F  1  ,PRU1K. 

I  PR  J  lt,HL0R,f.NGR,M>Hl91  ,rp,Ct>,BLl)Rl),F'({:R0) 

C 

C****«0NT2R  »RJTFS  OUT  THt  CAUCULATKD  MIXING  RFSULTS  A'lD  CURRtKT  UAf*. 


r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


100 


c 


»«*VARl*BLtS**» 

PRB  *  HASt  PRFSS'JRC  RATIO,  PH/PF. 

PRHIl  3  BASE  PReSSURC  RATIO,  PB/Pll. 

PHIJLUl  3  STAGNATION  PRESSl'RC  RATIO,  POE/POl. 

TRHOE  3  base  TLMPERATURF  RATIO,  TB/TOE. 

TRBOl  3  BASF  temperature  RATIO,  TR/TUI. 

THOEl  3  STACNAItUN  TEMPfRATURF  RATIO,  TOF/TOl. 

PROIL  3  INTERNAL  STAGNATION  TO  LxT.  STATIC  PRESS.  RATIO,  PUI/PE. 
PRlie  3  input  STATIC  PRESSURE  RATIO,  PlI/PE, 

BLORiLNCR  3  SEE  SUBROUTINE  'TUHIX*  FOR  PEFINITIONS. 

NPHINT  3  SEE  SUHROUTINE  *IN0UT*. 

CP  3  BASE  PRESSURE  COEFFICIENT. 

CU  3  BASE  DRAG  COEFFICIENT. 

RLORU.ENCMO  3  SMtClFlEU  VALUES  OF  THE  BLEED  AND  ENERCT  RATIOS. 


IF  (NPRiNT.LT.O)  GO  TO  lOJ 


nHITE  (1,100)  PRlIE,TROEI,PR0Eni,PR0IE,RLDR0,ENCR0 

FORNAI(19X,  AlH  ***»3*TURBULCNT  JET  NIXING  RESULTS* *»»•* ,  //, 


1  iOX,  19H  ♦•♦CURRENT  0ATA*»»  ,//, 

2  MX,  IIH  PII/PC  3  F«.5,  17X,  IlH  TDE/TOl  s  FI, 5,//, 

1  MX,  lIH  POE/POl  3  F8.5,  ITX,  IIH  PUl/PF  «  F8.I  ,//, 
4  t4X,  9R  HLORO  «  EI2.S,  ISX,  9H  ENGRO  3  EI2.5,//) 


♦RITE  (1,101)  BLDR,E9GR 

10)  EORnATOOX,  IBH  ♦♦♦NIXING  DATA***,//, 

1  MX,  8H  HLDR  3  02. S,  IBX,  BH  CRGH  3  02.S,/A) 

NRITL  (1,102)  TRBOC. TRBOl ,PRME  ,PRBII,CP,CD 

)02  FORNAT  (1BX,4SR  ••♦NASt  PRESSURE  AND  TENPENATURE  RESULTS***  ,//, 


1 

MX, 

lOH 

TB/TOt 

m 

rt.s. 

l«X, 

ION 

TB/TOI 

« 

f*,5,  // 

2 

I4X, 

tUH 

PB/PE 

m 

FB.5, 

IBX, 

lOB 

PB/PIt 

m 

Fi.S,  // 

1 

MX, 

lOH 

CP-B 

m 

Fi.5, 

IBX, 

lOB 

CD-B 

m 

FB,5  ,// 

4  20X,  40M  ♦•♦♦•♦END  Of  CURRENT  CASE  RESULTS* ♦•♦ ♦•  /, 

5  20X,  40M  ♦♦•**♦♦♦♦♦•♦••••••••*•••••♦♦♦•♦♦••♦♦♦♦•,  //) 

C 

101  RETURN 
END 


SUBROUTINE  1TEN(X,DX,ERRaRA,SIGN,Y,YGIVEN,ERRUHT,NI r,NTYPE, 
I  XNEG,YNEG,APnS,XPOS,NSIGN| ,NSIGn2> 


♦  ♦♦♦•SUBROUTINE  frERfOR“S  AN  |Tf  RATION  TO  FIH|>  X  SUCH  THAT  THE  AMSOl.iITt 
VALUE  OF  (T'tGlVENI  IS  LESS  THAN  OP  EUOAL  10  FKHORV  OH  THE 
ABSOLUTE  VALUE  UF  (  X  (  I  ♦  I  ) -X  (  I  )  )  IS  LESS  THAN  OP  E'JUAL  TO  FHRORA. 


♦♦♦VARIAHLFS*** 
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C  *  »  lui  EPr.MUKNT  V»H1A»*1.F. 

C  Ul  «  lM.Rt.‘'F*«r  |N  M)KNT  VAPIAHLF. 

C  ERhUR*  «  RAXI-ti-  VALUf  OF  AH<;  f  x  (  I  ♦  I  ) -X  (  t  )  )  FUH  .SOLliriON. 

C  SIGN  »  Al.c  mK  -1.0,  OtHRFS  INtHFRtUlIwG  FRDM  X  INITIAL. 

C  Y  >  IKPENDCNT  VAUIahLF.. 

C  YGIFCM  a  GIVEN  VALUF  OF  OtI  LtlotM  VARIAOLF. 

C  EHRORY  a  «AXIMUN  VALUE  OF  ABS ( Y- YGI VFN ) . 

C  NIT  a  InCNF*CNT  NONBFR, 

C  NTYPt  a  1,  increment. 

c  a  3,  Interpolation. 

C  a  3,  SOLUriON. 

c 

OYaY-YGIVEN 

IF(AHS  (1>Y)-FRR0HY)  90,90,10 
10  IF(DY)  30,90,30 
30  NSir.N3a-l 
XNFGax 
YNEGaY 
GO  TU  40 
30  NSIGN3a«t 
XPOSaX 
YPtlSaY 

40  GO  TU  (50,80),  MTYPE 
SO  IF(NlT-l)  70,70,00 
60  NSIGNsnSIGNI  «Ni!IGN3 
IF(NSIGN)  80,80,70 
70  NSlGNlaNSlGN? 

NlTaNlTAl 

Ct***a3NCRENENr  TO  FIND  SOLUTION  INTERNAL. 

XaXAStGN«OX 
GO  TU  100 
80  NTYPLB3 

C»*«»»INTENP0L*7I0N  FOR  SOLUTION. 

NITaNlI*! 

xsAvesx 

RA1IUa(XPUS-XNEG)/( YPOS-YNtG) 

XaXNCGARATIU«(YGI VEN-YNEG) 

C««***ACCELERAT10N  OF  CONVERGENCE  OF  ITERATION— REF,  aECSTEIM,  NBS. 

A  a  1.0/RATlU 
IF(A-l.O)  82,88,82 
82  0  a  A/(A-1.0) 

XaGSTN  a  U«XSAVC  t  (l.U-Q)*X 
IF(XNEG-XrGSTN)  84,86,88 
84  1F(X«GSTN-XP0S)  86,86,88 
86  XaXaGSTN 

88  IF(ABS(X-XSAVE)  -  ERRURX)  90,90,100 
90  NTYPEa3 
100  RETURN 
END 

SURHUUTINE  ABTS(GANNA,tNSI ,XRTI ,RHT1,AnGHT1 ,XHT2,KUT3,NSHAPE, 

1  NPR1NT,NUCPTS,NERR0R,CD) 

C 

C 

C  AXlSTHurTRlC  BOAITAII,  SU  BPROCHAH 

C 

c********»***$** »*,»****$* t»****t**%9***»*******» *»******»»*••***•*  »»•*» 

c 

C  WHITTEN  BY  —  A,  L.  AUUY 

C 

C  NOIE  — -  INPUT  AN’u  OUTPUT  HATA  ARI  FUR  1HK  •UPPEH-MALI*  'IF  FI.J- 

C  FlF.LP.  the  adjustment  IIF  TMFSF  IIATA  FUN  THI.  CALCUI.A  n  ii  46 

C  is  WAIIE  l«rFN'4ALLl. 


nno  o  o  r>  r>  n  r>  o  r>  r>  o  r>  o  r>  n  o  o  r>  r>  n  rj  o  n  o  o  r»  r»  r»  o  o 


SIlhPHUKH*"  KtOHIHtS  0UiBTI,l>.SP**,KPS,H[HPS,'«Cr)Ar*,llUTBT2, 

BTITf H. 


««*V*RIAHttS**« 

GAPMA  a  RATID  OF  SPACIFIC  HF.ATS. 

Ff>SI  s  iRITlAt  FRFtSTRtAM  MACH  STAR  AT  STATION  1. 

XBTI.HBTt  «  COUPUlfiATtS  OF  URST  lOTRT  OH  ROAtrAIL. 

ANCRTl  a  INITIAL  bOAMAlL  ANCLE  A)  STATION  I. 

NEGATIVE  AND  IN  RAUIANS. 

XB12,RBT2  a  FINAL  POINT  UN  BQATTAIL. 

NSHAPE  a  SEb  SUBROUTINE  *BTCNST*. 

NPRINT  a  •!  OR  0,  H.H.P.  UATA  NOT  PRINTED, 

*1,  H.B.P.  DATA  PNINTFD. 

N0CP1S  a  f4UNBER  OF  11-CHAR.  POINTS  CALCULATED  ON  CHAR.  THROUGH  (2) 
NERkOR  a  SEE  SUBROUTINE  *BT1TER*. 

*»»0U1PUT  DATA  (IN  URDf.R)*** 

INPUT  DATA  TO  ARTS 

X  a  I.ONGirUUlNAL  COORDINATE  UF  HOUNDARV  POINT. 

R  a  radial  coordinate  OF  BOIINPANT  POINT. 

THETA  a  LOCAL  FLOa  ANCLE  AT  BOUNDARY  POINT  (IN  DECREES). 

note  — •  THE  11-char,  data  THROUGH  (XHT2.Rbr2)  IS  TRANSMITTED  TO 
THE  MASTER  PROGRAM  THROUGH  aCOMHON*  IN  THE  ARRAY  CHARE. 


dimension  PMB(10U.S,2),  CHARKS.JD),  CHARE(5.)0),  P|(S),  P2(S), 

1  P)(S).  CII0(5) 

COMMON  PMH, CHARI, chare, PI, P2. PS 

PRMSF(EMS,GAMMA|a(1.0-((CAMMA-|.O)/(CAMMA«1.0))*eMS**2)** 

1  (GAMMA/(GAMhA-1.0)) 

EMNMSF(EMS,GA«MA)«50Rr  ( ( ( 2.0* (EHS** 2 > ) / ( CAMM Aa I .0 > )/ 

1  (1.0-((eAMMA-|,0)/<CAMMAa|.0))*(EMS**2))  ) 

CALL  HTCNST (XBri,RBTl,ANGBTl,XHT2,RRT2, NSHAPE, Cl, C2,C)> 
•«***INPUT  DATA,  SOME  0U1PUT  DATA,  ANU  COLUMN  HEADINGS  ARE  PRINTED. 
CALL  OUTBTI (GAMMA, EMS I ,XBT1 ,RHTI,ANGBT1 ,XHT2,RBT2, NSHAPE, 

1  CI,C2,C3,NPRINT1 

•A«««SET  INITIAL  VALUES. 

NGOTOal 

NOCPTSal 

NIal 

PRin|aPRMSF(EMSI .GAMMA) 

EMN|aCMNMSF(EMSI .GAMMA) 

M«««rUJMBER  OF  POINTS  CALCULATED  UN  THE  I 1 -CH AR ACTEKI SYI C  ORIGINATING 
AT  (XBr2,RBT2)  IS  SPECIFIED  HERF.  (LlMITE  MAX.  a  )a). 

LlNlTEaJO 

C****«FOR  UNIFORM  FLON. 

EMSaEMSI 

DRB0.U2*RbTI 

OXaSUHT  (EMN1*«2-I .0)«DR 
CO  TO  10 

C«*«**LOAD  INITIAL  VALUES  AT  (Xbri.HBTl)  INTO  THE  PMb  ARRAY. 

JO  PNB(l,l,I)aXBTl 
PMB(I,2,1)8RBTI 

PMMd,  J,l)aFMS 

PMB( 1,4,1 laO.O 

1F(AHS  (ANCPTn-r<Ul-3)  40»4O,5O 
C****«BUATTAIL  MITH  ZERU  INITIAL  TURNING  ANCU.F. 

40  KaO 


uuuuu  vuu 


i 


GU  TO  7n 

C*****rOH  AN  AFTLNBOOy  «ITH  INITIAL  TUNNING  ANCLE. 

SO  IF  (ANGHri )  S2.S3.S4 
C**«**FOR  A  HOATTAIL  (UtTA2C.  NEGATIVE). 

S2  K«(AhS  (57.2'»57a»ANGPTI  )♦!  .0) 

GO  TO  S6 

C«****AP)>RaXlMATE  ANAL7SIS  FuR  A  FLARE  (HFTA2K  POSITIVE). 

S4  K  >  I 
5b  FK  ■  K 

UTAsANGBTl/FR 

C**«**CALCULAT1UN  OF  CHAR.  ARRAY  OATA  fOR  POINTS  LnI.KlI  AND  ns). 

UU  bU  Lai. A 

PNb(L4l , I . 1 )3pH8(L, I .  I  ) 

PMP(LFt,2,l)aPMR(L.2,l) 

PHb(L«t ,4, 1 )aPNB(L.4. 1  )  ♦  UTA 

bO  PNH(L4t.i,l)aEMSPN(PNb(l.).l),PPB(l,4,l).P4b(L*1.4.l).GAMNA) 
C**«**R1  IS  number  of  (AMILY  II  CHAR.  FUR  SURDIVIOED  EXPANSION. 

70  KlsKtl 

C*****THE  INITIAL  RUUNUARY  POINT  DATA  IS  PPINTED. 

DO  80  Mai. 4 
80  P)(M)sPMM(Kl.N.l) 

CALL  OlJTBT2(CANMA.EMSl.EMNl.PRltJ|,P),Nt  .NCUTO.NPRInT.CD) 

**«*«THC  FLO*  FIELD  CALCULATIONS  ARE  NON  MADE  ALONG  FAMILY  1 

characteristics  STARTING  FROM  THE  INPUT  POINTS  ON  THE  SUBDIVIDED 
INITIAL  FAMILY  FOR  THE  FIRST  AND  SURSEOUENT  AXIS  POINTS. 

83  NlaNlFl 

*44««CALCULAT1UN  OF  THE  INITIAL  I I'CHARACTERISTIC  DATA  POINT. 
•A«««LOAO  CURRENT  II-CHARACTERISTICS  DATA  POINTS  INTO  PMB  ARRAY. 
PMR(|.1,3)8PMB(1.1,1)  ♦  OX 
PMB(l,2.3)aPMB(l,3.1)  4  OR 
PMRC1,4,2)8PMB(I ,4.1) 

PMB(1.3,2)aEMS 
GO  TO  (90,90.100),  NCOTO 
90  DO  92  Mai. 4 

93  CllD(M)aPMB(l.M.2) 

CO  TO  91 

94  00  9b  Mai. 4 

9b  PMb(l.M.2)aCI10(M) 

98  IF(K)  140,140,100 

***««CALCULATIUNS  ARC  FOR  THE  CURRENT  b>TH  POINT  ON  THE  INITIAL 
FAMILY  II  CHARACTERISTIC. 

too  no  no  lmi.k 

C*****CALCULAT1UNS  ARC  FUR  THE  CURRENT  L'TH  EXPANSION  MCHEmENT. 
C««*««lUAO  data/  field  POINT  CALCULATION/  STORE  UATA. 

CALL  MCDATA(l,L,L4l,L3.MPTS) 

CALL  FPSCCAMMA.PI .P3.PJ,NCRR0R) 

IF(NCRROR)  200.110.110 
no  CALL  MCDATA(3,Ll,L2.L4l ,KPTS) 

C««4*«ALL  FIELU  POINTS  ON  N-TH  FANILY  I  CHAR.  HAVE  BEEN  CALCULATED. 

GO  TO  (140,140.1201,  nguTO 
C*4«4«STORC  BOATTAIL  I 1 •CHARACTERISTIC  DATA. 

120  NUCPTSbN0CPTS«1 
DO  no  Mai, 4 

DO  CHARE  (M,NUCPTS)aPi(M) 

C*«**4CHARACTCH1ST1CS  OATA  SHIFT. 

CALL  NCDATA(J,LI,L2.L}.K*I) 

1F(R0CPTS>LIMITEI  82,200,200 
C«**«*LOAD  DATA/  BOUNDARY  POINT  CALCULATION/  ST'lHC  UA1A. 

140  CALL  MCDATAl 1,F«1 .N4l .L3.MPTS) 

CALL  BTBPS(CAMha,PI .PT.PIT.MSHAPF ,CI ,C2 ,C) • mCRHOH I 
1F(N».RR0R)  700,144,144 

E«****CUNT1NUE  BOATTAIL  CAICiaAIION,  flPHATF  FOR  I -CH  AR  \<- 1  r  H  1  S  1 1 C 
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r»  rt  r>  r>  on  r>  r>  r>  r»  o  ft  n  r>  rj  o  r>  r>  r»  r>  o  n  rt  n  n  ftftft 


IHHUIIGM  IHF  nfiAriAIL  FNU  PdlUT  <XMT2,MHT2),  IJH  CALCUL»Tf  TMt 
I  I'CHAfttCrttf  ISTIC  OHlGl»*Tir<r.  A1  THr  POINT  (XHI'2,KHr2). 

144  CAUL  BTirLR(XHTI,XHT2.PJ,Clllt,KGOTU.htPKi)H) 

IF(NbPROH)  200, Mb, 146 
146  GO  TO  (170,94, ISOl,  MGOfO 
«*«**LIIAD  F1KS1  POATTAIU  I  I •CHAKACTf 61ST 1 C  POINT, 

ISO  00  160  MSI, 4 
IbO  CMAHF.  (n,  1  )sPI(M) 

170  CALL  NCOArA(2.Ul,U2,K«2,KPTSI 

CUKRbNT  nUIOinAMY  POINT  DATA  IS  N(l«  PPlNTrO. 

CAUL  nuTeT2(CAMMA,»;-Sl  ,r.NM,PRt(i|,P7,Nl  ,Hi:urO,NPRlHT,CI<) 
M»»*CHARACTEPIST|CS  data  SHIPT. 

CALL  MCDATA( 1,UI ,U2,U3.A42) 

•  •••AAOVAOCb  INtir.X  FOR  nlXT  INPUT  POINT  (IN  INMIAL  CH  IP  ACTLH 1  ST  1 C . 
KsKaI 
GO  TO  82 
200  PETURb 
END 

SUbROUTlNb  (1ICNST(XHT1  ,KbTl  ,AN(;pT1  ,XnT2,PHT2,NSHAPF.,Cl  ,C2,CJ) 


B 


5 


XRTl  a 
RBTI  a 
ANG8T1  a 
XRT2  a 
RBT2  a 
NSHAPC  a 


ci,r2,cj 


Initial  longitudinal  koattail  cuoroinatc. 

INITIAL  RADIAL  OOATTAIL  CUOROINATE. 

INITIAL  UUATTAIL  TURNING  ANGLE,  RADIANS,  CCa(«). 
TERMINAL  LUNCilTUOINAL  HUATTAIL  COORDINATE. 
TERMINAL  RADIAL  BOATTAIL  COORDINATE. 

1,  OGIVE  BOATTAIL. 

2,  parabolic  BOATTAIL. 

J,  CONICAL  BOATTAIL. 

a  COEFEICIENTS  IN  THE  BOATTAIL  PROFILE  EDUATIONS. 


SLOPEls  TAN  (ANCBTl) 

GO  TO  (10,20,341,  NSHAPE 
M«*«0G1VC  BOATTAIL  (NSHAPCal). 

10  Cls(0.S)*(  (XBT3«XBT1 IMAT-T.'lbSLOPEI FRHTl •(XBT2>XBT1 MKBT2**2 

1  -RPT1**2)  /  (RBT2-RBTI-l,0»SLnpri*(XHT2-XBTl)  ) 

C2a  X6T1  ♦  SLUPL1*(RBT1«C1) 

C3a  (XBTI«C2)*«2  ♦  (R8ll-Cl)*42 
GO  TO  40 

l***«PXRAeOLIC  BOATTAIL  (NSRAPta2), 

20  Cla(  RBT2«RBTt«SLOPLl«(XBT2«XBTI I  )  / 

1  (  XBrM«34XBT7««2  -2.0*XI*TiaXBT2  ) 

C2aSLUPCl  •2.0FCMABTI 

C3SRBTI  •  (  C3«XBTI  «  Ci*(XHTI**2)  ) 

CO  TO  40 

laFFACONICAL  BOATTAIL  (NSRAPtal). 

30  ClapBTI 
C2aSLOPCl 
C3aXBTI 

RMT2>RBT1«SL0PEI*(XBT2>XBT1 I 

40  NETIIRN 
END 

SUBROUTINE  pUTBTl(CANNA,EMSI,XBTl,RnTI  ,ANr,BTl  ,Xttr2,MMT2,NSHAPF, 
I  C1,C2,C3,NPRINT) 

MaasTRIS  SUBROUTINE  PRINTS  INPUT  DATA,  SliMF  dUTPUT  DATA,  AND 
HEADINGS  Fur  the- boa (TAIL, calculations. 

PPN5F(EMS,CAaMA  )ai  I  .  »•(  (  GA«m«.  |  )  2  (  GA  «  »  a  4  1  . '• )  I  FF  I  j  •  •  2  )  *  * 


B-27 


nnnnonnnonoonnnnnrjrinnnnnrin 


1  ((;*Mi4A/((;AHhA-i  .0) ) 

KM.4MSF(e*S,0AMHA)«SilkT  ({{2.0»tr«'.S*»2))/('JAMM**|.i)))/ 

1  (i.n-((<;*H«*-i.o)/(f:****<**j.o))»(Ki«s»»2j)  ) 

IF(SPRINI)  10,10,1(10 
100  t,»<N|stMNMSf  (F.MSI  ,(;aMH*) 

PHim*»>hK5F(e“S!  .GAA«A) 

Bt:TAD«S7.20&ri*AN«.HI  I 
C 

ANlTt  0,11  GAM(‘A,FFM,PHtni 

1  FOAMATC  Ihl  ,///,21  X  ,2  IH  AXlSFkHfTFlC  HOATTAIl,  /, 

1  I5X,30H  WITH  IlhlFOMM  SUPEMSOMr  FMih  //, 

2  21X,iUH  IkPUr  0A1A  •••  //, 

3  7X,0H  GA«***  s  FS.3,3X,i2H  kACH  k(t.  *  FS.J,3X,  8H  P/F(l  ■  Fk.4//) 
C 

GU  TO  (2,4,6),  NSHAPF 
C 

2  kOITK  (3.3) 

3  FOPkATdM  ,|9X.27H  •  OGIXE  HOATIAII.  pRIlFlt.b  •) 

GU  TO  • 

C 

4  write  (),b) 

5  FORMATOH  ,19X,32H  «  PARABOLIC  PUATTAIL  PROFILE  •  ) 

Gt'  TO  e 

C 

6  fcRlTE  (3.7) 

7  FURkATdH  ,19X.30H  «  CONICAL  HOATTAIL  PHCFILE  •  ) 

C 

•  -RITE  (3,9)  XBrt,RttTl,HETAD.XnT2,RB12,CI,C2,C3 

9  FORMATdH  ,//,7X.  8H  XBT)  >  F6.3,3X,  RH  RHTl  ■  F6.3, 

1  4X,10H  AkCRTl  a  F 6 . 3// . 7X , BH  XBT2  a  F6.3,3X,IH  RnT2  a  F6.3//, 

2  7X,IM  Cl  a  r7.3,2A.RH  C2  a  FT.3,3X,10H  C3  a  F7,3///, 

3  20X.37H  •**  BOAfTAlL  SURFACE  OUTPUT  DATA  •••  //, 

4  13X,lHX,14X,lHR,lnX,RHNACH  NO. , 9X , 4HP/P1 , 9X , 9HCP( LOCAL )  //) 

C 

10  RETURN 
CNU 

SUBROUTINE  RTHPS(GAMMA,P1,P2,P3,NSHAPE,C) ,C2,C3,NCRR0R) 

BOATTAIL  BOUNDART  POINT 
SUBROUTINE  (BTBPS). 

THIS  SUBROUTINE  CALCULATES  A  POINT  P3  UN  THE  HUATTAIL  nALL 
GIVEN  THE  PROPERTIES  01  A  POINT  PI  IN  THE  FLOn  FIELD. 

•VVVARIABLES*** 

GANNA  a  RATIO  OF  SPECIl IC  HEATS. 

PI(J)  a  J>TH  FLUM  VARIABLE  AT  THE  POINT  1  aHERE  la), 2, OR  3. 
PI(J)  AND  P2(J),Jat,S  a  FLOW  VARIAHLES  AT  KNONN  POINTS  I  AND  2. 
P3(J),Jal,9  a  FLUB  VARIABLES  AT  THE  UNKNOWN  POINT  3. 

THE  J  SUHSCHIPT  INDICATES  THE  FOLLOnlNC  VARIABLES  — 

Jal  CURRESPUNUS  TO  X. 

Ja2  CURRESPUNUS  TO  R. 

Jal  CONRESPONUS  TO  naCH  STAR  (FNS). 

J«4  CORRESPUNnS  TO  THETA  IN  RADIANS  (IHET). 

NSHAPC  a  SFE  BELUa. 

C1,C2,C3  a  CONSTANTS  IN  THE  HOATTAIL  PROFILE  COUATIIINS. 

NFRROR  a  A  CVINTRUL  VARIABLE  FOP  CHFCRINr,  THE  PUSSlHlLlTt  (HAT 

the  current  CHARACTF.RISI  IC  kISSFS  THE  h'lAIIAIL  ANii  AN 

iteration  IS  required. 

NERNUR  a-1  ...  ERROR. Ik  C ALCULAT I UN . 
nERROR  a-  u  ...  NU.ITtPATION  REUUIREO. 

NF.RRUR  a  I  ...  AN  ITFpATiUk  |  S  HFQIlIRKU, 


POINTS  t  AhP  i  ARt  ASSU»>EO  CONNECTED  BY  f  AMltY  I  4HbRE 
PUiNT  }  IS  ON  TME  ■ALL. 

THE  HOATTAIL  PROf  lLf.  IS  SPECjritl*  BY  fOUATt'JNS  OF  THE  FOHH--- 

1.  IF  NSHAPE>|  OGIVE 

R»  Cl  ♦  SuPTl  Cl  -  (X-C?)»*2  ) 

2.  IF  NSMAPEB2  PANAP'ILIC 

R»  Cl  ♦  C2**  ♦  C1»(X**2) 

I.  IF  NSHAPCal  conical 

H«  Cl  ♦  C2*(X-C3) 

taUfcPt  Cl,C2,ANn  Cl  HAVE  BtfN  CALCULATE  FNOa  THE  INPUT  DATA 
IN  SUHRUU1INE  tHTCNST*. 


OlNbNSlUN  PI(S).  P2(b)«  PKS) 

ALPHAF(LNSTAR,GANMA)aATAN  (SURKd.O  -  ( ( G  AHN  A- 1 . 0 )  /  (  CAHh  A  ♦  I  .  0  )  ) 

1  *(EMSTAH«*2n/(LNSTARaa2«1.0))) 

AVGKA.B)  ■  (A  «  n)/2.U 

PCOEFFIEMSTAH. ALPHA)alMSTAH*TAN  (ALPHA) 

OCUtFF (NPOI NT, RADIUS, ENSTAH, THETA, ALPHA )8(( EAST AH /RADIUS)* 

1  (IAN  (ALPHA)«*2)*TAN  (THt:TA))/(TAN  (THETA)  ♦  ( ( - 1 . 0 ) • •NPOl NT ) • 

2  TAN  (ALPHA)) 

HUCUEF  ( RAOIUS,CNSTAH, THETA, ALPHA )a((EMSTAR/RAUIiJS)*T AN  (ALPHA)* 

1  SIN  (ALPHA)*SIN  (THETA)) 

C***«*NP01NT  IN  OCUEFFO  INDICATES  THE  KNOmn  POINT  HEINU  USED— 1  UK  2. 
C*****ENBOB  FLAO  SET. 

NFRRORao 

NCOONTaO 

nctnaxbis 

EMSMAXaSQHT  ( (UAHNA* 1 .0 )/(CAMHA«1 .4 ) ) 

c«****xNUaN  Input  data  frun  points  i  and  2. 

XlaPKl) 

RlaPl(2) 

ENSlaPI ( 1) 

THETIaPKA) 

H2aP2(2) 

CNS2aP2( J) 

THCT2aP2(4) 

CaaaaaFUR  AN  INITIAL  ESTIMATE  OF  THE  VALUES  AT  PUINT  1. 

R3aAVCF(Rl,R2} 

EFSlaAVCriEMSl ,CNS2) 

THET3aAVCF(THETl,THET2) 

GO  TO  17 

C***«*1TERATI0N  FOR  VARIABLES  AT  POINT  3. 

C*****IF  HSHAPE  a  ],  OGIVF. 

I  Aal.U  *  (TAN  (D1FF|I))**2 

B«2.0*(R1-CI)*T*N  (UIFF13)  -2.0*C2-2.0*X I  * ( T AN  (DtEF13)  )**2 
Cm  C2**2  -  CJ  ♦  (  ( R 1 -Cl  )-X  I*T*N  (DIFF13)  )**2 

D1SCR«B**2-4.0***C 
IF(DISCN)  14,14,3 

3  X1b(-H-SUBT,  (B**2-4,0*A*C))/(2.0*A) 

R34RI*(X3-X|  )*TAN  (IIIFFI3) 

THET3aATAN  (  (C2-X 3  ) / ( K l-C 1 )  ) 

GO  TO  10 

C*****1F  NSHAPE  a  2,  PARAHOLIC# 

4  AaCI 

haC2-TAN  (OlFFII) 


uuuuuuuuuuuuu 


l)ISChs|:>«*2«4.n*««C 
ir(OlSCR)  I9.l9,k 

6  X3a  (-B*StiRT  (B**2-4.<i***C)  )/(2.U«*) 

BJaRMtXJ-xn^TAN  (Oimj) 

THIT3S41A»  (C3«2.0«CI*Xi) 

(iO  TU  to 

C»»»»4If  NSM4t>t  «  3,  C04ICXL. 

7  X3a  <CI-Ht-C2«C3»Xt*T*N  (Ult»'l3)  )  /  CTAn  (DfFKlJ)  -  C2  » 
M3aRI«(X3-Xt  )*TAN  tl>lFF13) 

IFtBtl  19,10,0 
9  THtT3«ArAA  (C2I 

C**«*«TkST  AnU  'f.VAUUATIUK  FOH  HUBtZUtilAt  1  •CHAW  »C  ICR  IS  7  ICS  . 

10  IFIABS  (UlFFl  l)-l  .UF.-3)  11,11,12 
C««*«*FOP  I  HOHIZONTAL. 

11  PHODllaHOCuer  (Rl3,tHSll,THtT13.ALPH|3)*(X3>Xi ) 
r.U  TU  13 

C***«»FOR  1-CHARACTtKlSriC,  O.R, 

12  PRUD13»iCueFr( t ,H13,tMSl l.THtTI 3,ALPH1 3)*tR3-Rl ) 
c**««4CAi.cuLAriuN  ur  fluo  vahiablils  at  poimt  i. 

13  F.MS3«CHS1-Pt  3«(THeT3-THbTl  IwPROOli 

oirrNSa(k.HS3-sAvti  i/sAvei 

IFdrMSI.LT.l.OI  .UR.  (eHSS.GT.F.HSMAXt)  GO  TO  20 
iriABS  (UIFFMS)  .Lb.  1.0F.*41  GO  TO  IB 
17  NCOUNraOCOUNT*! 

IFINCOUNT  .GT.  nCTMAX)  GO  TO  IR 
SAVCl  a  EMSI 
R13aAVCr(N| ,R3I 
fcMSlSaAVCriEHSl ,eMS3) 

THCT13aA¥Gr(TlieTl,THeT3) 

A|,PH13aA|,PHAr(EMS13,CAMMA) 

OIFFllaTHtTll-ALPHlI 
P13«PC0Gf P(eMSI3<AtPHll) 

GO  TO  (1.4,7),  NSHAPC 
IB  P3(l)  a  X3 
P3(2)aR) 

P)( llaEMSl 
P3(4)aTHET3 

ir(RCOU«T  ,GT.  NClPAX)  aRITC  O.IRO)  t»COUNT,OirroS 
180  FORMAT!/,  SX,37H  «•*  CONVCRCENCb  KRRUR  IN  *BTnPS*,  (  ,I3,2H  ,  , 
1  E10.3,0H  I  /) 

RETURN 

19  NCRRUNa.l 
RETURN 

20  NCRRORa.t 
NRITC  (3,21) 

21  rURMAT(//.23X, )2H  4**  error  I«  •HTBPS*  CAtC.  •••  //) 

RETURN 

CNU 

SUBROUTINE  0UTHt2(GANRA,ENSl .ENNl.PRint ,P3,N| ,NGUru,NPHlNT,CU) 

aaaaaTHlS  SUBROUTINE  PRINTS  THE  CAt.CUl>ATF.U  ROATTAlt,  SUMFACE  DATA 
AT  THE  LOCATION,  Na  NUBPTS,  IN  THE  RPTS(M,N>  ANRAT. 

•4«vARlA8LCS*4a 

GAMMA  a  RlflO  OF  SPECIFIC  HEATS. 

EMSl  a  FREESTREAh  m«CH  STAR. 

EMNt  a  FREESTREAH  RACH  NUMHER. 

PHIUI  a  FREESTREAH  SI  A  nC'Tli^STAGNAT  ION  PRESSURE  RATIU. 

P)(J)  a  BOATTAIL  BUUNUARt  POIM  DATA. 

THE  J  SUBSCRIPT  1NUICA1FS  THF.  FULLOnIHG  VANIAHLFS-*- 
Jal  CURRFSrUNUS  TO  X. 


B-30 


uuuuouuuuuo  uuo  uwu  u  u u o u u u u u u u u u u u u 


js2  c■JMnts^Ml^rs  iu  fc. 

J«J  C'iHKKSP'INUS  TU  ***fH  ST»H  (F»*.S), 

Jb4  CUPKFSPnNdS  TU  TMtTA  IN  MAl>l*NS  (fHer*). 

Ni  ■  1,  ...  LuCATLS  the  nuuNOAHr  PMiNf  UN  rnt  BOArTAit 

SIJKFACC. 

NtiUTO  s  I,  NOhAAL  HUATIAIU  CAl.rULAT10N. 

B  i,  ITCPATIun  FUN  l-CMAHACTEHlsriC  rHt<UUr:H  (XR12,PNT2I. 

B  1,  CAtCUUAIlON  UF  11-CHAHACTKHISTIC  THHUJf;H  (XRTJ.KHT?) 
NFRINT  B  SbF  .SUBPUIJIINL  «AHTS*. 


DlNFNSIUh  P}(S) 

PH“SF(tNS,GABNA)B(  l.0-(  (CANNA-1  .0  >  / 1  GANN  A  ♦  1 . 0 )  )  •EXS  •  •  2  )  •  • 

1  (GAPMA/(CAMMa-| .0)) 

tNNFSf  (EXS.GANMAIB  SUf)T((  (2.o*(tH.S**2>)/('iXNMAAl  .n)  )/ 

1  ( 1 .0-1 ICANNA-I.O J/CGAMmaaI .0)  )»(txS»»2) )  ) 

IF(ltPHINr)  HO, HO. to 
10  XBP1,(|) 

Hst i(21 
F.NSBPi(3) 

tPNsCHNM&FCLNS. GAMMA) 

PPOnuiBl .0 

PRBlB(PRMSF(CMS,GAMMA)/PHiai )*PRUI)OI 

LUCAL  PRESSURE  CUEIFICIENT  IS  CALCULATED.  CP  IS  BASED  UN 
THE  FREESTRCAH  MACH  NliMBCR  AND  PRESSURE. 

CPb(PRB|*1.0)/(P.S*CAMMA*(EMN|**2)) 

HRITE  (3.20)  X.R.EMN.PHBI.CP 

20  FllRMAT(7X,riO.S.SX.ri0.5.SX.F10.S.SX.FIO.S.SX.riO.S) 

•A««*THC  BOAITAIL  DRAG  CUEFFICIENT  IS  CALCULATED.  CD  IS  REFERENCED 
TO  THE  FKEESTREAM  PRESSURE  AND  MACH  NUMRER  CUNDITIUNS. 

IF(N1«1)  30.30.40 
•44441NIT1ALIZE  CO  CALCULATION. 

30  cnao.o 

OeNOMBO.S*GAM«A«(CMNl**3)*(RBB2) 

GU  TO  SO 

40  AVGPKB(0. S«(PRMSF(F.NSL.CAMMA)«PRMSF(CHS. GAMMA) )*PRUR01 )/PR101 
CnBC04(((1.0-AVGPRl*(RL*A2-R**2)»F0EN0M) 

SO  RLbR 

CMSLbCMS 

GO  TO  (H0.R0.60).  NGOTO 
60  bHITC  (3>70)  CO 

70  FOHMAT(/,2SX.2aH  DRAG  COEFFICIENT,  CD  b  F«.S,3H*B4  ,  //) 

80  RETURN 
END 

SUBROUTINE  NTI TER (XBTl .XBT2.P3,C1ID,NGUTU, NERRUR ) 

•••««SUBROUTINt  CONTROLS  BOATTAll.  ITERATION  FUR  I-CHARACTERISTIC 
PASSING  THROUGH  (XBT2.RBT2). 

•••VARIABLES*** 

XBT2  B  LUNGITUDINAL  CUURU.  OF  TERMINAL  POINT  OF  THE  ROATTAIL. 
P3  ■  CURRENT  BOUNOAKT  POINT  FROM  SURRDUriHE  •HTbPS*. 

CUD  B  CURRENT  INITIAL  1 1'CHARACTERl STIC  UATA  POINT. 

HCOTO  >  1,'BOATTAIL  CALCULATION. 

a  2,  ITERATIUN  FOR  I'CHAM ACTFH I  ST IT  THHUUGH  (XHT2,Rur2). 

B  3,  CALCULATION  or  II-CMARACTERI.STIC  THROUGH  (XHr2.HHr3) 
NERHOR  a  •! ,  ERRUH  IN  IICrATIOR.  GO  TU  NEXT  CASE, 
a  0,  HOUNOAHT  POINf  CALCULATION  n.A. 

a  I.  ERROR  IN  MUUNUARY  POINT  CALCULA rillN ,  S1AHT  ITERATIUN 


uuuuwuuuuuuuuuuu 


OIMFNSIUN  PJ(S)f  S«VtL(S)>  i>AVl.H(b).  C1II<(S) 

XHT  ■  (XI»T2>XB1I) 

C«****eKRI)R  OR  ITCRAIIUN  ritlECTlOH. 

GO  TO  (I0.6O),  NGUTU 
10  IF(NFRHOM)  20,20,bU 
30  IF(XH12-P3(I ))  SO,190,iU 
30  ITENal 

DU  in  ••>1,4 
4U  SAVEL(H)«CI  I0(**) 

Rl TURN 

C»**»*1TEBATI0N  SFO'JFmCE, 

50  NG0T0>2 

00  ir(NERROR)  TO.TOtllO 

70  lF(*llS((XbT3-Pl(l))/XhT)-I.UE*4)  19u,i9n,R0 
•0  lF(XHT3«Pi(l))  110,190,90 
90  on  100  M>1,4 
100  S*VEL(H)aCllD(N) 

GO  TO  1 30 
110  no  120  Nal,4 
130  SAVFHlNlaCl  10(>t) 

130  IFIITF.R-ISI  100,160,140 
140  NERHORa«l 

NRITC  (3,150) 

150  FURN*T(//,5X,67H  *•«  MAX.  NO.  ITERATIONS  FXCEEOEO  IN  SHR.  BTITER. 
1  GO  TO  NEXT  CASE.  //) 

RETURN 

160  1F(ABS  ((SAVEL(l)*SAVER(l))/XBT)*1.0F-4)  190,190,170 
170  ITERalTEK*! 

CaaataiMTERVAl,  H«t,ve  FOR  VALUES  ON  INITIAL  1 1  •CHAN  ACTER ISTIC. 

DO  ISO  N8l,4 

ISO  CIlD(M)aO.S«(SAVEL(N)«SAVER(N)) 

RETURN 

CvaaaasOLUTION  FOUND. 

190  NGUTOS3 
RETURN 
END 

SUBROUTINE  UFLOC(GANNA,ENS, AC.RC,N1 ,CHAR , NFLON ) 

C 

C**a«aTHlS  SUBROUTINE  SUHOIVIUES  THE  INITIAL  FANILT  II  CHARACTERISTIC 
ANU  CALCULATES  THE  INRUT  DATA  FOR  POINTS  ON  THIS  CHARACTERISTIC 
FOR  UNIFORM  FLOH. 

aaavARlABLES*** 

GAMMA  a  RATIO  OF  THE  SPECIFIC  HEATS. 

ENS  a  approach  MACH  STAR. 

XC  a  LONGITUDINAL  COORDINATE  mhENE  EXPANSION  IS  CENTERED. 

RC  a  radial  coordinate  «hf:re  expansion  is  centered. 

NEGATIVE  FOR  INTERNAL  FLOa  AND  POSITIVE  FOR  EXTERNAL  FLOa 
N1  a  number  of  INCREMENTS  OF  INITIAL  CHAR.  (MAX.  IS  29) 

CHAR  a  initial  CHARACTERISTIC  DATA  ARRAT. 

NFLON  a  1,  internal  FLUH. 

a  2,  EXTERNAL  FLO*. 

DIMENSION  CHAR(9,  )0) 

EMNMSF(CMS,V;AMMA)aSUHT(((3.0a(EMS*a3)  )/(r.AMMAA|.i>|  1/ 

1  (|.0-((CA*N*-l,il)/lCA‘*MA»1.0)l*(ENS»a2))  ) 

GO  TO  (10,20),  NFLOa 
CataaaFOR  INTERNAL  FLQm. 

10  N|a|5 

FN I aN I 


I 


nnnnnnnn  n  nnnn  r%  n  n  n  n  n  n  n  n  n  n  n  n  n  n  n  n  n  n 


t)M3*BS  (XO/rNt 
cu  to  lu 

C*****rni>  LXrtRNAL  FLUti. 

20  DH30.0J**BS  (RC) 

3U  |)XsbR*S0RT(  (tMNRSFttHS.GAM***)  )*»2'*l  .0) 

NPTS*N|*1 

PU  40  NsItNPTS  I 

FW3N-1 

CHAP  (1,N)  ■  XC  ♦  fN*l)X  - 

CHAR  (7,N)  ■  PC  *  FN«PR 

CHAM  (3,N)  X  RMS 

40  CHAR  (4,N)  m  n.n 

MFIUHR 
tNO 

SIJBROUriNr  CMFLOC<CARHA,fcHS,RFTA,XC,Rf ,H1 ) 

*«**«F0R  internal  CilNICAL  FLOP,  THIS  SUBROHTINF  SIIHOIVIPES  THE 

NUN'CHAPACTFHISTIC  UNIFORM  FLOa  CURVE  THROUCH  THE  PUlNT  (XC,NC) 
AND  THEN  CALCULATES  THE  INPUT  UATA  ALONG  THE  FAMILY  II 
CHARACTERISTIC  aHJCH  ORIGINATES  AT  THIS  POINT. 

SUBROUTINE  REUUl RES*--FPS , APS. 

444VARIARLES444 

GAMMA  a  RATIO  OF  THE  SPECIFIC  HEATS. 

EMS  a  APPROACH  MACH  STAR. 

BETA  a  ELOH  ANGLE.  NEGATIVE,  IlN  RADIANS),  AT  (XC.RC). 

XC  a  LONGITUDINAL  COORDINATE  NHERE  EXPANSION  IS  CENTERED. 

RC  a  RADIAL  COORDINATE  NHERE  EXPANSION  IS  CENTERED. 

N1  a  number  or  increments  op  INITIAL  CHAR.  (MAX.  IS  29) 

dimension  PMB(100,S,2),  CHARMS, )0),  CHAREIS.IO),  PH5),  P2(S), 
i  PJ(S) 

COMMON  PMB,  CHARI,  CHARE,  PI,  P2,  P} 

RCONEbrc/SIN  (BETA) 

••444SUB01VIS1UN  or  THE  NON-CHARACTFRISTIC  CURVE  INTO  N2  INCREMENTS. 
(N|a2*N2).  TO  change  THE  NUMBER  OF  INCREMENTS  CHANCE  ONLY  N2. 
(MAXIMUM  N2  IS  14). 

N2a|0 
rN2aN2 
Nla2«N2 

•••**STnRE  INITIAL  DATA  POINT. 

PMBd.l.lJaXC 
PMe(1,2,t)aRC 
PMH(l,3,naEMS 
PMB(l,4,l)aBCTA 
DU  to  Mai, 4 

to  CHARl(M,l)aPHB(l,M,l) 

•••4«THC  PLUM  FIELD  CALCULATIONS  ARF  non  maDF  ALONG  FAMILY  I 

CHARACTERISTICS  STARTING  FROM  THE  POINTS  ON  THE  SUHDIVIUEO 
NON.CHARACTERISMCS  CURVE.  THIS  SEDUFNCE  IS  NUT  APPLICABLE  FOR 
CALCULATIONS  INVOLVING  OTHER  THAN  THE  FIRST  AXIS  POINT. 

***«aTHF  CALCULATE!'  FLOm  FIELD  DATA  FUR  TRF  (NMl)  PJIvrs  ON  1  Ht 

FAMILY  11  CHARACIFRISTIC  ORIGINATING  AT  (XC,HC)  all.L  HF  SlURFD  Al 
CHARKN.M),  NHERE  N8|,N1«). 

DU  40  Ns|,n2 

C«*«**CALCULATL  DATA  ON  THE  NON-CHAN At TtRI ST) C  INPUT  CURVF . 


uuu  u  u  u  u  u u u u o u u u o u u u u u u u 


*NG(.eRsHK1«*(  I  .0*f  N/F  Iii2  ) 

PMp(N«l ,  1 ,7)sxC«HCUN»  •(Cns(*'ilGI.KK)>COS(nKT*) ) 

HHHI N*1 ,2,2)sHrONb«SlN( ANGLf R) 

PRe(H*l,4>2)s**<GLEK 
RPr5»N*l 
DU  20  Ia|,N 
LSNoltl 

C***«*LIJ*0  U*T»/  CRI.CULAIf.  rULb  PUINT/  STOPR  OAF*. 

CAI.L  MC0ATA(l,L*ltL,L.3,KPT&) 

CALL  FPS(GAMi4A,P1  ,P2,P).Nr»<K0R) 

CALL  HCUATA(2,LI •L2,L,KP1S) 

20  CUNTlNUt 

C**«**STaRL  INiriAL  CHANACTLHISllCS  DATA. 

[10  30  W>l.« 

30  CHARl(M,K«i)=PMa(l,M,2) 

C*****SHIFT  MbTHOL  nr  CHAHACteRISTirS  DATA. 
call  HCUArA(3.Lt«L2,L3.KPTS) 

40  CURTINUL 

«$«««THC  CALCULATION  SCUUkNCE  IS  NON  MiJbll  IF.D  FUR  SUBSRUUbNT  AXIS 
AND  netO  POINT  CALCULAIIUWS. 

00  90  Nal,N3 

MI«N24N 

LaN2*l*N 

***«*LnAD  DAIA/  CALCULATF.  FIGLO  POINT/  STORK  DATA. 

CALL  NC0ATA(1,L.L,L3.APTS) 

CALL  APS  (GAMMA, P3,P3,NERRnR) 

CALL  MC0ATA( 2,L1 ,L2,L,AP1S) 

1P(N1>n1)  70,10,50 
SO  NIIaL't 
LII«L 

00  90  Ial,NII 

«««*«L0A0  UATA/  CALCULATE  FIELO  POINT/  STORE  DATA. 

CALL  MC(>ArA(l,LII,Lll>l,L3,RPrS) 

CALL  FPS(CAMMA,Pi,P2,P3,NCRNOR) 

CALL  MC0ATA(3,L1,L2,L11«1,APTS) 

90  LtlaLlI-l 

•4444STORE  INITIAL  CHARACTERISTICS  UATA. 

70  DO  80  Mai, 4 
80  CHARI(M,Nl«|)apMH(l,H,2) 

•♦•••SHIFT  METHOD  OF  CHARACTERISTICS  DATA. 

CALL  HCOATA(3,Lt,L2,L3,L) 

90  CONTINUE 
RETURN 
END 

SUBROUTINE  PMSRR ( GAMMA , EMSTAR ,PRATIO, BETA , XC , RC , K ) 

••*»(THIS  subroutine  SUBDJVIUCS  THE  INITIAL  PR ANDTL'METER  EXPANSION 
(NAVES  or  family  IIJ  INtU  APPROXIMATELY  I  UECREC  ImCHEMLNTS. 
INPUT  DATA  IS  THEN  CALCULATED  FUR  THE  METHOD  OF  CHAKACTENISTICS 
NET  AT  THE  POINT  nnEHE  THE  EXPANSION  IS  CENTERED. 

SUBROUTINE  REUUIRES--«FMSPM. 

•••vahiablcsa** 

I 

gamma  a  RATIO  OF  SPECIFIC  HEATS. 

EMSTAR  a  APPRUACH  MACH  STAN. 

PRATIO  a  EXPANSION  PRESSURE  RATIU  (P/PU). 

BETA  a  Initial  FLOm  angle  In  radians. 

XC  a  longitudinal  CnUPDINATF  MHFPF  EXPANSION  IS  CENIEHKI). 

RC  a  RADIAL  COuHUInATE  mmenF  FXPANSinN  (S  CEK^mlx. 


'^cnrinnooonrtnon  n  n  n  n  rinr>r>orir>r>r»r>nr»r>r»n 


*  NiiHBtR  ur  iNChtxfT.s  nr  int  thuxi^c  *'<Gbt. 
a  «  i-n|MF.‘«SlUli«L  *MKAy,  H*<f<  ( t ,  a ,  m  )  ,  OF  UATA  FOM  THt 
MLTHOU  OF  CHAHACILRISTIC?  a(.T.  THC  SUttSCK  IPTS  L.A.A 
HAVE  THE  FULLUxlNb  HANf.tS  AMO  MtAVInGS--* 

L8|,K*1  AND  COHHtSPONUS  TO  THE  l.-TH  PUIVT  OF  THt 
SHBOtVIUEO  PHAhOn.'-HrYEK  EXPANSION. 

Mai  CnpHLSPONUS  TO  X. 

Maa  COHPfcSPtINOS  to  H. 

Ma3  COHHtSPONOS  TH  "ACM  STAR  tF»S). 

Ma4  CORKESPONOS  TO  THtTA  IN  RABIANS  (THETA). 

Nat. 2  CURRF.SPUNl;S  TO  THE  PREVIOUS  JH  CURRENT  l-CHAR. 
L.Nal  at  point  mHERF  THE  INITIAI,  FLO*  CONDITIONS  ARE 
SrtClFlEU  AUD  THE  P-M  EXPANSION  IS  CENTERED. 


dimension  pMR( |0U,S.2).  CHARKS.JU),  CHARE(S,10),  Pl(b),  P2(S) 
1  r3(S) 

COMMON  PNB,  CHANT,  CHARE.  PI.  F2.  PI 

OMCCAF( A,il)aSURT( (B«l .0)/(H-1.0))*ATAN  ( SUNT ( ( A  * * 2* I . 0 ) / 

1  ((e»l.O)/(B-|.0)-A*»2)))-AIAN  (SURT(((B*I .0)/(H-t.O) )• 

2  ( (Aa*2-1 .())/( CBAl .0)/(B-I,0)-A»*2)))) 
F.MSPRF(A,a)aSURT(((B*I.0)/(B«l.O))a(l.0-A**((B-1.0)/B))) 

EMSiaEMSTAR 

EHS2aEMSPRF(PRAT10,CAMMA) 

MA*«roR  mavcs  of  family  11. 

AriGLERa-(UMECAF(LMS2.GAMMAI  -  OMEG AF ( EMS I . GANN A ) > 

IF  (ANGLEH)IO, IU,20 
10  Aa(ABS  (S7.29ST8«ANGLEM)«1.0) 

GO  TO  JO 
20  Rat 
JO  FKaK 

OELTAaANGLER/FR 

MAAARnOmN  INITIAL  INPUT  DATA  FUR  PNB  ARRAY. 

PMBCl.t.DalC 
PMB(1.2.1)aRC 
PMBCI , J.l laEMSI 
PMBd.A.UahCTA 

••••tCALCULATION  OF  ARRAY  DATA  FOR  POINTS  Lal.Ftl  AND  Na|. 
no  1  Lal.R 

PMBILtl , 1 • I )aPMH(L, I ,t) 

PMB(L*l,2,t)apaR(L,2.l) 
pMp(L*l,4,|)8PNR(L,4,t)  *  UFLTA 
I  PmB(L«I , J.l )aEHSPM(ENSl.PMP( t.4,l).PMB(LA| ,4,1),GAMMA> 

RETURN 

END 

FUNCTION  ENSPM ( ENST AR. THETA I, THFTA2, GAMMA) 

•«*a*THIS  FUNCTION  CALCULATES  THE  FINAL  MACH  STAR  AFTER  A 

PRANDTL-MCVER  expansion  or  compression  GIVEN  INITIAL  M* 

AND  THE  turning  ANGLE  It*  RADIANS. 


••a  VARIABLES  ••• 

ENSPM  a  FINAL  MACH  STAR  AFTFR  THt  TURN  OF  (THE1A2  • 
EMSTAR  a  A«>PR0ACR  maCH  SIAR. 

THETA!  a  APPROACH  I LOo  ANCLE  (IN  RAdlANS). 

THF.TA2  a  FINAL  FLOm  ANCLE  (In  HAOIANS), 

GAMMA  a  RATIO  OF  SPECIFIC  HEATS. 

THF  SIGN  CONVENTION  FOR  ANGLFS  IS  C»(-)  »HI<  CC*(*). 


THETA  I ) . 


IJMEUAKA.Hla  S'Jt- r  (  («♦  1  .'•)  )  »»TAN  (  SuKtC  (*»»2-l  .0)/ 

1  (tB4l,0)/(B-I.Ul-A«A2)))-ATA*l  I  SOHT(((t.4l,n)/(H.l.O))» 

2  ((A»*2-».0)/((H4l,())/(h-l.l))-A»*2)))) 

c*»»»*stT  iNiriAi,  values. 

M  i  r  »  0 

>iir>*AA  a  20 
•*nPEa| 

C«4***l<T>Pr=l  .  iNTtHVAL  HALVE.  NT»PEa2,  1  MT  tPrf'L  A  TE  . 

NATIIlaO.) 

AN(;LE«(THLrA2-rHE  r  A1  ) 

IFIANLLE)  20,20,10 
C****»rt)P  A  HEVEPSIHLE  COHPRESSIUN. 

10  LHSHal.o 

OHECAHaii.O  . 

EHSPataSTAH 
GO  TO  JU 

C***«»EUB  A  HFVtRSIHbt  EXPAHSIOH. 

20  EHSNsEMSrAR 

Oa(GANa(lMk.GAF(EHSH  ,GAMHA  ) 

EMSPa  SORT ( ( G Ahm A4 1 . O ) / (GAHM A«1 .0) ) 

CIVAAAEVALUATE  UMELA  FUNCflllN  EUH  CONIilTlOM  *2*. 

30  UHF.CA 2a (UHEGAFICMSIAR, gamma )-ALGLE) 

C««**«I)UES  THE  SULUTKIN  EXIST. 

1F(UMCGA2)  40,60,70 
40  MHITE  (i,SO) 

50  FORMAT!//, |0X,25H  »♦*  ERROR  IN  -EFSPM-  /) 

RETURN 

60  CMSPMal.O 
RETURN 

INI TI AtLT  INTCHVAL  HALVE  AND  THEN  INTFRPULAfE. 

70  NIT  a  NIT  ♦  1 

IF(N1T  ,CT.  NIIMAX)  GO  TO  140 
EMST»tNSN4R4T10*(EM5P-EMSNl 
UNFGATaOMCGAF(CMST,GA«HAl 
01FFU8(0HtGAT«UMEGA2)/OMEGA2 
1F(ABS  (OlFFOI'l .0C«4)  140,140,00 

80  IFIUIFFU)  90,140,100 
90  ENSNaFMST 

UMEGANaUMEGAT 
GU  TU  110 
100  ENSParasr 

(IMEGAPaUMCGAT 

hT»PE«2 

110  OIFFMS  a  (C«SP<-ENSN)/LMSN 

irCAHSlOlFFMS)  •  1.0E«4)  140,140,120 

120  GO  TO  (70,130),  NTIPC 
C*a***iNTt:RPaLATC  FOR  THE  SOLUTION. 

1  JO  NATI0a(UHEGA2>0MCGAN)/ (OMEGAP-OMEGAN) 

GU  TO  70 

C«aaaasULUT10N  FOUNO. 

140  ENSPNaEMST 

iriNll  ,CT.  NlTMAXI  6N1TE  (3,150)  NlT.DlFFU 
150  rURNAT(/,5X ,34H  • * (CON VERCFNCE  EHNOK  IN  FmSPm,  (  ,  13,  2H  ,  , 
1  ElO.J,  6H  )  •a*  /j 
RETURN 
END 

SURNUUTINE  OurHlII  (H,NPKINT,I»PTS) 

C 

C«****SUHRnuTINE  PRINTS  THF  CURRENT  CALCULATED  HJJNPART  POINT  DATA. 

c 

c  a«*vAR)AbLFS«*« 

C 

Cm  a  NMMAF.R  OF  CDRhenT  m«iiin(>arv  point. 
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n  r>  o  o  o  n  o  n  o  n  on  no  nnnonnnnn  on  n  nnnnnnnnn 


Nk'xlM  *  -I  III.  o,  C.P.H.  l/ATA  liUl  IklNItH. 

♦I,  C.P.B.  PAT*  PkiurMi, 

HkTSCK.M)  s  CUkPtNT  oUUNUAHY  fiAIA. 

M=|  niHHtSkUNUS  TO  X. 

MS?  CURPtSkOkOb  TO  R. 

MsJ  CUI.R(  SPUNt'S  TO  HACH  star  fKMS), 

Ms4  ClIMktSPUNU.S  TO  IHtTA  tA  RADIARS  (THKTA). 


UlMtRSIUS  HPTS(%,iO) 

IF(NPRINT)  2,2,1 

1  XsHPiS(l,N) 

RsRPTS(2,N) 

THe  rA«S7.?9&7  8*BPrS(4.N) 

WRITL  (1.10)  X,  R,  THtlA 

10  KORMAHFlS.6,  F2'».b,  FJO.bJ 

2  RFTURN 
CNO 

SUBRUUTlAt:  MC0ATA(KiUP,Ll  ,  t.2 , 1. 1 , K PTS ) 

l**««SllHRUUTlNt.  LOADS,  SlUHtS,  OR  SHIFTS 
MtTHOD  OF  CHARACTLHISTICS  DATA. 

NOP  ■  1,  LOADS  PNH  DATA  IN  PI,P2. 

■  2,  SrUHFS  Pi  DATA  IN  PNB. 

■  3,  SHIFTS  PNB  DATA  FHON  1-2  TO  1-1. 


OINCNSION  PMa(IOO,S,2],  CHARKS.IO),  CHARk(5,30),  PUS),  P2(5) 
1  P](9) 

CONNON  PNB,  CHARI,  CHARC,  PI,  P2,  PJ 

GO  TO  (to, 10, SO),  NOP 

10  DO  20  MSI, 4 

PI  (n)spmh(i,|  ,N,2) 

20  P2|N)sPNB(L7,i4, 1  ) 

RLTURN 

JO  PU  40  MSI, 4 
40  PMB(L},M,  2)sP3(M) 

RtTURN 

SO  no  70  RlIsl.KPTS 
DO  60  MSI, 4 

60  PNB(A11,N,  DshmbIRII ,H,  2) 

70  CONTINUE 
RETURN 

END 

subroutine  FPS(GAHMA,P| ,P2,Pl,NFNRnR} 

>*«**AXISTNMCTNIC  FIFLO  PUINT  SUBROUTINE  (FPS) 

SS*VARIABLFS*«« 

GANNA  s  RATIO  OR  SPECIFIC  HEA1.S. 

PI(J)  s  J-TM  FLU»  FARIAMLE  AT  THE  POINT  I  •MtHE  |s|,2,nR  I, 
PI(J)  and  P2(J),Jst,4  s  FLHs  VARIANLFS  AT  ANUS «  POINTS  I  AnO  2 
PI(J),JB|,4  B  Fl.ua  VAHIABLCS  AT  Till  MNXIui-a  PUIMT  1. 


norjonnrjnooorj 


THE  J  SUbSCHlFT  IHlMCAftS  THE  EliLl.n^IHG  VamIABLES--- 
Ja|  CDAKESHUnUS  TU  K. 

Js2  COHHESPOf.US  TU  H, 

jBj  CKRHrSPUNOS  TU  haCH  iiTAH  (EKS), 

J«4  CXRKE.SPUNO.S  TU  THETA  IN  RADIANS  (THET), 

riERROK  a  -1,  ERRUR  IN  CALCULATION, 
a  0,  CALCULATION  U.K. 

POINTS  I  and  1  ARE  ASSUHP.L  CUNNFCTEO  BY  fAAILY  I. 
POINTS  2  ANU  3  ANK  ASSUNCIJ  CONNECIEU  BY  FAMILY  11. 


DIBENSION  Pl(5),  P2(5),  P3(5l 

ALPHAI  (|:mSTAR,CAMHA)zATAN  (SbRTdl.O  •  (  (  C  AMM  A»  1 . 0 )  /  (  GaMm  A«  1  .  ii )  ) 

1  *(EM,STAR*A2I  )/(LMSTAN*«2-l  .0})) 

AVGF(A,B)  a  (A  A  Bl/2.a 

PCUY.FFCCMSTAR, ALPHA )aCMSTAK*TAN  (ALPHA) 

OCUFFI  (NPUINT.RADIUS.CMSTAH.THETA,  ALPHA  )s(  UMSTAR/KADIIIS)* 

1  (TAN  (ALPHA)AA2)«TAN  (THETA)  )/(TAN  (THE'TA)  ♦  (  (  -  1 . 0  )  •  ANPUI  NT  )  » 

2  TAN  (ALPHA)) 

HOCUCF  ( RAUl US. CMSTAR, Theta. ALPHA )al (emSTAH/HAUIUS) (TAN  (ALPHA)* 

1  SIN  (ALPHA)*.S1N  (IHCTA)) 

(;***««nPU1NT  if  OCUEFPO  INDICATES  THE  KNOaN  PUINT  BEING  USCD->1  OR  2. 
C*»***tRROR  FLAG  SET. 

NCUUNIaO 

NCTNAXatS 

NERRORaO 

EHSHAXaSOHT  ( (CAHN A« 1 . U ) / (CANHA-I .0) ) 
c*««««knumn  Input  data  fhoh  points  i  anu  2. 

XlaPKI) 

RI=PI (2) 

F.MSiaPK  J) 

THtTIaPl(4) 

C 

X2aP2(l ) 

R2sP2(2) 

EHS2*P2( i) 

THET2ap2l4) 

C«*««*fOR  initial  estimate  of  average  values  BET-EEN  points  I-J  ANO  2-3, 

RisAVGFIRt ,R2) 

ChS3**VCF(EMS1 .ENS2) 

THET3=AVGr(THETI,THEI2) 

C(t  TO  11 

C*****lTERATtUN  fur  VARIABLES  AT  PUlNT  3. 

1  X3a(R2  •  Rl  ♦  XI«TAN  (UlFFt3)  -  X2«TAN  (5UM23))/ 

1  (TAN  (01FFI3)  -  VAN  (SUB2))) 

R3«(RI  ♦  (X3  -  Xl)*TAN  (01FFI3)) 

C««***TEST  AND  EVALUATION  FOR  HORIZONTAL  I  OR  II  CHAR AC Tf R 1  ST  1 CS . 

IF(ABS  (OIFF13)-I.OE-3I  2,2,3 
C«***«FOR  1  HORIZONTAL. 

2  PROOI3aHUCOtF  ( Rl 3 . ENS 1 3 , IHET 1 3 , ALPH I  3 ) • ( X 3-X I  ) 

GO  TO  4 

3  PROOlIaOCOEFFd.Rl J,tHS13,THETI3.ALPH13)»(R)-Rl ) 

4  IFIAHS  (SUB23)-1.UC-3)  S.S.N 
C**«**FUi>  II  HURIZUNTAL. 

b  PR0033sHUCnEF  (R23 , EHS23 ,THLT23. ALPH23 )* ( X 3-X2 ) 

GO  TO  7 

6  PROD2iaOCUEFr(2,R23,t'’523, THET2J,*LPH73)*(M  I-R2) 

C«**«*CALCULAT1UN  UF  FLOr  variables  A1  POINT  ). 

7  ThETJa(P13*THF.TF  ♦  P23*1HEI2  ♦  PROOIJ  -  PROUZJ  ♦  E^SI  •  t-S2)/ 

1  (PI34P23) 

EHS3aFBSI  -  PI  1*( THEI )-THETI  )  ♦  PPUPl  ) 
niFFMS  a  (EMS)-SAVLI )/SAVFt 


B-38 


IK  (f'Si.LT.  I  .III  .Ut..  )  )  uii  iU  li 

1K»MS  IUiF»«S)  .Lt.  I.UF-4)  W  TM  I? 

C 

11  NC'iuNiB.vcuui<r*i 
IKNC0liNT.GT.NC1‘<*X  )  GtJ  fO  17 
SAVri  a  {.MSi 

HI  IsXVUKKI  ,»j  J) 
n73a*V(.KH7,Ri) 
t<*Sl  3«*VGf  (►JXSI  ,FHSJ) 
t.MS7  3sA«(iKF*«S7,L»'Sl) 

THFTl  JsAVGK  IMEl  I  .THfeTJ) 

IHKriiaAKGKTMrra.TMtl  J) 

ALPHI 3aALPHAF(e  xbl I.CAMxa) 

ALPH7  3sAl,PHAr  (K<S3  I.GAMmA) 

PI  3aPCUFIK(ePSn,ALPHI  3) 

P23aPCUEFr(fMS73, ALPH21) 

OlFFJ  3aTHtTI 3-AtPMl J 
SIJM2i8TH».T7  3*ALPH2i 
r.O  TO  I 
C 

12  P3n)  >  XI 
P3(2)»R3 
P3( 3)«t«SJ 
P3{4)aTHCT3 

IFIMCOUNT  .GT.  NCTPAA)  aPITt  (1,920)  NCOUNT , D I FF9S 
920  FUPMATC/,  SX.ISH  CONVePGtNCI  CPttUR  IN  *FPS*,  (  ,I3,2H  ,  , 

1  EI0.J,6M  )  •••  /) 

RtTUHN 

C 

II  NEHHUNa-l 

•RITE  (3.14) 

14  rURMAT(//,33X,29H  *•*  ERROR  IN  AFPS*  CALC.  *•*  //) 

HF.TORN 

END 

SdbPOUTlNE  APS  (GAMMA, P3.P3,NCRRUR) 

C 

CfAAAAAXISTNMCTHIC  AXIS  POINT  SUbKOUTINF  (APS) 

C 

C  FOR  THIS  SUHROHTINE,  THE  UNKNUMN  POINT  3  IS  UN  THE  AXIS. 

C  THE  KNOaN  POINT  2  ANb  THE  UNRNUaN  POINT  3  ARE  ALONG  FAMILY  II. 

C 

C  AAAVARIABLES*** 

C 

C  gamma  a  ratio  of  SPECIFIC  HEATS. 

C  PI(J)  >  J'TH  FLOm  VAHIABLR  AT  THF  POINT  I  MHERE  (a|,2,UH  i. 

C  P2(J),J*t,4  a  ri.O*  YARIABLES  AT  KNOhN  POINT  2. 

C  Pl(J),Jat,4  a  FLO*  VARIABLES  AT  THE  UNKNOWN  POINT  3. 

C  THE  J  SUBSCRIPT  INDICATES  THE  FULLUalNG  VAHIAHLES  — 

C  Ja|  COPRCSPONUS  TO  X. 

C  ja2  CURNESPONOS  TO  R. 

C  JaJ  CORHESPUNOS  TO  maCH  STAR  (EmS). 

C  Ja4  CQRHESPONOS  TO  THETA  IN  RADIANS  (THED. 

C  NERROR  a  •!,  ERROR  IN  C ALCUI.ATIOn. 

C  a  U,  CALCULATION  O.E. 

C 

C 

DINRHSION  P3(S),  PKS) 

ALPHAr(EMSTAR,GAMMA)aATAN  (SUHTdl.O  •  ( ( GAMMA- I . 0 ) / ( GAHM A f  I  .  0 ) ) 
1  *(EMSTAR**2)  )/(r.MSTAH*»2«I.O))) 

AVGF(A,B)  a  (A  »  H)/2.0 

PCOE»F(EHS1AR,ALPHA)aEMS1 AH«TA«  (ALPHA ) 

OCUFFF(NPO|NT,RA1)IiiS,EmS^AH,TMH  A,  alpha  )a(  ( KNS  (  AM /M  AH  I  IIS  )  • 

1  (TAN  (ALPHA)AA3)A1  AN  (THETA  ))/(  TAN  (IhKAI  ♦  (  (  -  1  . 'i )  A  •  NP'i  I  •- I  I 


nnnnnnnnnnnn 


2  1*N  (At.PH*)) 

C**«**NPUINT  in  DCt'KFr  (  )  INOICATFS  iHf:  KNllxN  POINT  ItKiNC  USk:0>>l  ON 
C*****tNRUN  FL*0  SbT. 

NCUUNT  >  0 

NCrxAXalS 

NbRRORaO 

CNSNAXaSURT  (  ( G* nm * ♦  1 . 0 ) / (r;A-N*-|  . n ) ) 

C**«*«KNriwN  input  data  fur  POINTS  2  AND  1. 

X23P2(1  ) 

R2aP2(2) 

LMS2aP2(3) 

THt.T2aP2(4) 

PJaO.O 

THETJall.O 

C*«««*fUR  initial  FSTINATt  UF  AVF.RAGF.  VALUFS  HFTWbEN  POINTS  2  AND  1. 
EHSJaLRS2 
R2Js*VCF(R2,R3) 

TRFT23aAVGF(THEr2,THbT3) 

GO  TO  5 

C««***ITCRAtIUN  FUR  VARIAHLCS  AT  POINT  3. 

1  X3aX2  •  (N2/rAN  (SUM23)) 

EHS3aEHS2  •  P23«THCT2  -  a23*R2 
DIFFNS  a  (eaSI-SAVED/SAVei 

IFdCMSI.LT.l.O)  .OK.  (EHS3.GT.FRSNAX))  CO  TO  7 
IFIAbSlOlFFNS)  .LC.  l.UE«4)  GO  TO  b 

S  NCUUNTaNCUUNT^l 

IF(NCUONT.GT.NCTNAX)  GO  TO  b 
SAVC|aeMS3 

EaS23aAVGF(LNS2.EMS3l 

ALPH23aALPHAF(ENS23,GAMNA) 

SUN23aTHCT23«ALPH23 

P23aPCueFF(eNS23.ALPH23) 

Q23aOCUCFF(2,R23,iMS23,THeT23.ALPH23) 

GO  TO  I 

b  PI(naX3 
P3(2)aR3 
P3(3)«F.mS3 
P3(4)aTHET3 

IFINCOUNT  .GT.  NCINAX)  MRITF  (3,bO>  NCUUNT , D1 FFMS 
«0  FOMMATI/,  SX,3SH  CONVERGENCE  ERROR  IN  *AP5«r  (  ,I3»2H  ,  , 

1  CI0.3,bM  J  •••  /) 

RETURN 

7  NERRORa-t 
mRITE  (3,8) 

8  FORMAT(//,23A,29H  ERROR  IN  aAPSa  CALC.  //) 

RETURN 

END 

SUeROUTlNE  CPBSCGAMMA,  pi,  P2,  P3,  NFHROR) 

**«*«AX1STNMETR1C  CONSTANT  PRESSURE  BOUNnART  SUHROUTINE  (CPBS) 

POINTS  2  and  3  ARE  ON  THE  SAME  CONSTANT  PRESSURE  rtilUNUART. 
POINTS  I  and  3  ARE  ASSUMED  COtlNECTFO  NT  fAMlLT  1. 

a**VARIAHLES««* 

GAMMA  a  RATIO. UF  SPF.CIFIC  HEATS. 

Pl(J)  a  J.TH  FLO*  VARIaVlC  AT  THE  POINT  I  nHENF  la|,2,OR  ). 
Pl(J)  AND  P2(JI,Ja|,4  a  FLOa  VARIAHLF.S  AT  KNUaN  POINTS  I  AND  2. 
p)(«l),ja|,4  a  FLua  VARlAHLbS  AT  IMF  lINKNOitN  POINT  ). 


THt  J  SUHSfMIPT  I'U'ICAILS  TH»  FOLt.O^lNG  V  AM  1  AHLF.S'-* 
Jsi  COMMb&PUAOS  TU  X. 

Js2  cuhmksponus  to  P. 

Js3  COKRtSPUNUS  TU  HACH  STAR  (F.HS). 

Ja4  CURMKSPUNUS  TO  THF.1A  14  RADI  AMS  (THCT). 

NfRRUR  X  -I,  FRRON  IN  CALCUl.ATIUN. 

X  U,  CAI.CIILATIUN  O.K. 


OIRFNSIUN  PH5),  P3(S).  PJ(5) 

ALPHAF(E.xSTAR,GAR.'iA)aATAN  (SORTCd.O  -  (  ( G  ARM  A- 1 . 0  )  /  (  GAMM  A  4 1 . 0  )  ) 

1  •(e«STAh«»2})/(bMS1AH*«2>1.0) )) 

AVGF(A,HI  a  (A  «  h)/2.0 

PCUerr ( FRST ah , alpha )xtaSTAH*TAN  ( alpha ) 

HOCUFF  (RADII'S, (NSrAR,THFtA,ALPHA)a( ( f mstA R/R ADI US ) *T AN  (ALPHA)* 

1  SIN  (ALPHA)*S1N  (THETA)) 

UCULFFl NPUI NT, RADIUS, LmSTAK.THLTA, ALPHA )x( (KHSTAR/RADIUS)* 

1  (TAN  ( ALPHA)**2)*TAN  (THF:TA))/(TAN  (THETA)  ♦  ((-I.0)**NPU1NT)» 

2  TAN  (ALPHA)) 

C««*«aHPUlNr  IN  UCOF.FFO  INDICATES  THE  KNUHN  POINT  BEING  USED— 1  ON  2 
C*****CRNnH  FLAG  SET. 

NCOUNTaO 

NCTMAXalS 

NCRNURaO 

CTAttXKNUNN  INPUT  DATA  FROM  POINTS  I  AND  2. 

Xlxpi (I  ) 

Rlxp|(2) 

EPSIaPl ( J) 

THET1xp1(4) 

C 

X2«P2(1  ) 

N2aP2(2) 

EMS2aP2( J) 

THbr2«P2(4) 

C<««««FUH  INITIAL  LSTIRATE  OF  AVERAGE  VALUES  BF.TnECN  POINTS  !•)  AND  2>A 
R)aAVGF(HI ,R2) 

ThCT38AVCF(THCTI,THCT2) 

C«*«**SINCE  points  2  AND  J  ARC  UN  THE  SANE  CONSTANT  PRESSURE  BOUNDARY, 
EHS3>LNS2 

LRSl3aAVGr(F,RSI,EHS3) 

ALPHl3xALPHAF(eNS13,GANMA) 

P13xPCOEFF(CNSI3,ALPH13) 

CO  TU  t 

C****«1TERATI0N  for  VARIABLES  AT  POINT  3. 

1  X3«(R1  -  H2  ♦  X2*IAN  (THCT23)  -  Xl»TAN  (DIFF13))/ 

I  (TAN  (IHbT23)  •  TAN  (D1FFI3)) 

R3*(R1  *  (X)  •  XI)«TAN  (OIFFt3)) 

SIGN  a  Ri«SAVCl 

C««***1F  SIGN  IS  NEGATIVE  OR  AERO,  AN  ERROR  HAS  OCCURRED. 

IF(SIGN)  I,*, 2 

C«***«TCST  and  EVALUATION  FUR  HORIZONTAL  I-CHARACTCRI 5  TIC. 

C**«a«FOR  1  HORIZONTAL. 

2  IFIABS  (UlrF13)«I.Oe«3)  3,3,4 

3  PH0013aH0Cacr  (RI3,bNS13, THETIS, ALPH|3)*(X3-XI) 

GO  TU  S 

4  PROOI3auCUKFF( 1,H| ),CMSt 3,THFT13,ALPHt  3)»(H )-RI ) 

$  THET3a(THETl  •  ( ( EMSS'EMSt'PRODI 3) FPI 3 ) ) 

UlFrTa(THCT3>SAVC3)/SAVE2 
iriABSIDlFFT)  .LE.  l.Uf4)  GO  TU  7 
C 

B  NCOUNTaNCUUNT*'! 

IFINCOUNT.GT.ncTnax)  go  to  7 


u  u  c;  u 


SkVt  i«THk.l  1 
HI  JsAVUFOII  ,R)) 

THtri 3sAVGr(THLTl .THtT)) 
blKFliaTHtrtJ-ALHHI J 

OUsUCUtIFC  I  ,RI  J.LHSI  ),tHF.TI},ALI>Ht  )) 

THf:T2i*AVGI  (THr:T2.1Hk.Ti) 

GO  TO  1 

T  t>J(l)3X3 
P3(2)>HJ 
R3( 3)«eHS3 
P3(4)arHk;f  3 

iriNCUUNT  .GT,  nCTHAX)  •HITF  (3,70)  WCniKi T , OIFF T 
70  FORMATC/,  SX,3bH  **•  CUNVtRCCNCE  ERROR  IN  «CPBS*,  (  ,I3,2H  ,  , 

1  ElO.J.bH  J  «•*  /) 

return 

H  NERHOKS'I 
1.P1TE  (3,<0 

9  FUKNAT(//,23X, JUH  ERROR  IR  ACPRS*  CALC,  //) 

RETURN 

ENU 

SURROUriNE  0UTPUT(GAMMA,CNS1 ,PRAT10,RLTA,NPK1NT,NPLUM) 

*«««*SUHRUUTINE  prints  input  and  sure  OUTPUT  DATA,  AND  COL.  HEADINGS 
FU)^'''  THE  AXISTHNETRIC  CONSTANT  PRESSURE  BnilNUART  SUBPROGRAM. 

EMNPRF(PH,GAHHA)aSuHT((2.0/(GAHMA>1.0))* 

1  (PR»*<-(GANMA»1.0)/CAMMA)-|.n)) 

EHSNNF(ENN,GANNA)aSaRT((0,S«(CANNA»t.0)*(EMN**2) )/ 

1  {1.0F0.5*(CAMNA-),0)*rENN»*2))  ) 

CMNMSF<CNS,GAMMA)aSORT(((2.0*(e»>S**2))/(CAMMAiI.O) )/ 

I  (1.U-((CAMMA»I,0)/(CAHMA*I,0))»(EHS**2))  ) 

C 

IF(NPRINT)  70,70,10 
10  HETA0aS7.29S779S*HE1 A 

EMNl  a  eMNMSF(EMSl , GAMMA) 

EMN2aCMNPRF(PRATIO,GAMHA) 

E«S2aCMSMNF(eMN2, GAMMA) 

GO  TO  (20, SO),  NFLOM 
C 

20  1F(AHS  (BETA)>t.UE-4)  30,30,40 
C 

30  hRITE  (3,100)  gamma,  BETAD,  CMNI,  PKATIU, 

1  PRATtO,  FMN2,  EMS2 

100  FORMAT! Ihl,  tn,  21X,  31H  CUNSTANT  PRESSURE  JET  BUUNDART  /, 

1  I9X,  IbH  FOR  INITIALLY  UNIFORM  AX1>SYMMETR 1C  /, 

2  24X,  2SH  SUPERSONIC  INTERNAL  FLUa  //, 

3  2«X,  17H  •••INPUT  OATAP^^  //, 

4  7X,  9H  GAMMA  a  FS.3,  24X,  ISH  HETA  (PEC.)  a  FlO.b  //, 

5  7X,  12H  MACH  NO.  a  ER.b,  ITX,  *M  P/Pll  a  F8.6  H, 

6  22X.  27H  •••BOUNDARY  OUTPUT  OATAP^A  //. 

7  7X,BH  P/PO  a  FR.B,3X,1|M  MACH  M'.  aF9.b,3X,)3H  MACH  STAR  aF9.6//, 
B  7X,  2H  X,  27X,  2H  R,  23X,  I 3H  THETA  (DEG.)  /) 

C 

GO  TO  70 
C 

40  MRITE  (3,101)  gamma.  PETAP,  EMN|,  PHAriO, 

I  PHATIU,  feMN2,  EM57 

101  FURMATOHl,  ///,  21X,  31H  CONSTANT  PRESSURE  JE|  R.liir.UAKI  /, 

1  19X,  3bH  FOR  INITIALLY  (!UN|CA|.  AX T'SYMMVTR  1C  /, 

2  24X,  2SH  SUPEhSI'NIC  INTERNAL  FL<>«-  //, 

3  2RX,  17H  •••input  UATA*^^  //, 


U  O  (J  u  u  u  u  u  u  u  u  u  u  u  u  u  u  o  u  u  u  u  u  u  u  u  u 


4  7X,  911  (;»•««*  =  f'l.i,  74X,  till  I-MA  (ttft:.)  =  KlU.n  //, 

b  7X,  I2M  maCM  Nil.  a  fU.e,  I7l,  NH  H/r<i  a  FH.fc  //, 

h  i2t,  2;ri  ***r*iiUNOXH«  UUfl'llf  //, 

7  7X,I(H  P/FO  a  FM.b,lX,llH  NACH  NO.  3K<I.K,U.I2H  4ACH  SfAH  afU.A// 
H  7X,  2n  X,  27X,  2H  H,  2JX,  1  )H  IHUA  llifr..}  /) 

on  TO  »i) 

50  nKITP  (3,1U2)  OAmnA,  bFIAIi,  PHATM, 

1  PRAIIIl.  »PN2.  LMS2 

102  FURMATdHI,  ///,  2IX.  31H  CONSTANT  PHISSUPC  Jtf  OHMOAKy  /, 

I  19X,  ihH  FUR  INiriAbLy  UNIFOHM  A X I -S YNNFTR I C  t, 
i  24X,  25H  SUPFRSONIC  LXTCRNAL  M.UN  //. 

J  2f'X,  17H  •**irjpllT  0ATA«44  //, 

4  7X,  9H  OAMVA  a  F5.3.  24X.  15H  (iFTA  (IFG.)  s  FIO.s  //, 

5  7X,  12H  MACH  NO.  a  F9.b,  17X,  NH  P/PO  a  FO.b  //. 

b  22X<  27H  **«HUUNUARy  OUTPUT  oATA*«*  //, 

7  7X,FH  P/PO  a  F8.b,3X,llH  MACH  NO.  BF9.b,lX,l2H  MACH  STAH  aF9.b// 
H  7X.  2H  X,  27X,  2H  K,  23X,  iJH  THFTA  tOFC.)  /) 

70  KFTUHN 
KNn 

SDPHOUTINF  TtSr(PLMr.NSIMr«NFtOb,N,HPTS) 


XSUHMOUTINe  STOPS  CALCULATIONS  ANO  PCTIIRNS  TO  THE  MASTER  IP  ••• 

1.  THF.  INTEHNAL  BOUNUART  HAUIUS  EXCEEDS  KLMT  UR  IF  THE  JET 

nouNpAHy  angle  chances  sigh. 

2.  THE  EXTERNAL  HOUNDAHI  RADIUS  IS  LESS  THAN  RLMT. 


dimension  BPTS(5«I0) 

CO  TO  (10.30),  nFLOn 

10  1F(HPTS(2.M-RLMT)  2U,50,S0 

20  IF(ePTS(4.N.l)«BPTS(4,N))  SO. SO. 40 

30  IF(nPfS(2.N)-RLMT)  so, SO. 40 

40  NSTMTnI 
CO  TO  bO 

SO  NSTMTB2 
bO  RETURN 
END 

subroutine  SLIP(ENSI .THETA 1 .GAMMA], EMS2,THETA2,CAnMA2, 

I  THETAS, NSTOP) 

THIS  subroutine  calculates  the  SLIPLINF.  ancle  for  THE 
OBLIQUE  SHOCK  PECOMPHESSIUN  SYSTEM  nHICH  OCCUNb  AT  THE 
IMPlNCMENr  POINT  OF  T«n  SUPERSONIC  STREAMS  IF  IT  EXISTS. 
ONE  or  THE  SHUCKS  MAI  BE  A  STHONC  SHOCK. 

SUBOUTiNt  RE0UNIRLU>««-STSHK 

•A44VARIABLES**** 


EMSI  a  MACH  STAH  OF  STHEAM  1. 

THETAl  a  FLOa  AuCLE  OF  STREAM  |. 

CAMMAl  a  ratio  OF  SPECIFIC  MEATS  F"N  STRFAM  1. 
DF.LTl  a  SHOCK  TURNING  ANCLE  FOR  STPVAM  i, 
OOELTI  a  IIIHNINC  AncLF  UEMIVATIVE  OF  STRFAM  I. 
LMS2  a  MACH  STAR  OF  SIREAM  2. 


no  ooo  non  oooon  ooooo  oooo  oooo  oooooooooo 


THKTA2  »  fLOii  AN'iLK  UA  SThtA**  2, 

UAMHAi  ■  RATlU  UF  SFFrtAlC  MaaTS  FUH  STWFA*«  2. 

UFLT2  >  SHUCK  TUHHINC  ANCI.L  FOR  STRbAH  2. 

DUeLT2  a  IuRNJnG  AnCLF.  ntHlVATlVt  OF  STHFAM  2. 

THFTAS  a  SLIPLlat  ANGLE 
NSTUH  at.  fur  A  ;>ULUTION 
a  ],  A  OR  HU  SULuriUN 

.^Utr  THAT  THETAl  IS  ASSUMEU  LAHGFH  THAN  THLrA2. 

i«NHSF(CHS,GA*tMA)aSUHT((2.0/(GAAMAAl  .0)  )«(t.HS««2)/ 

1 1  l.n-HGAHHA-1.0)/H.AHHA*l.0))*ttHS»»2))) 

CALCULATIUN  of  the  HAXTmUH  rUHNlNG  ANCLE  FUR  A  SIlfEN 
APPROACH  MACH  NUMHEH  AMO  GAMMA;  NACA  R- 1  MS , EO . 1 bH . 

SIMt<A2lEHN,GAHHA;a(».2S/(GAHHA*(tMM«*2)) )•( (GAMMAfl .0)* 

1  (LMM««2)*4.0*S'lMr((GAHMAAl.0)A|  (CAMMA«  !.())*  (EMN««4)A|t.O* 

2(CAMHA-1.0J*(CMN««2)Tlb.O))) 

S1NHA2  CALCULATES  THE  SINE  UF  THE  SHUCK  »AFE  ANCLE  SQUARED 
FOR  MAXIMUM  STHEAM  DEFLECTION  HfHlND  THE  SHUCK  (EON  |bb) 

OELTAH(EMN,GAMHA,SIN2MA|aArAN  C ( 2.0*SnRT ( ( 1 . 0-S I H 2* A ) /S1N2MA ) 
1«((EMN«*2  J*SlN2HA«|.on/(2.0«CEHH*A2)t((;AMiATl  .U>2.0*SlN2aA))) 

UELTAM  CALCULATES  THE  MAXIMUM  TUNNING  ANGLE  GIVEN  THE 
APPROACH  MACH  NUMHER  GAMMA  AND  THE  SINE  SUUAHEO  OP  THE  aAVE 
ANCLE, SlM2MA<FnH  THE  MAXIMUM  UEFLECTION  (EON  1I9A), 

PR(ISHR(tMN,SIN2MA,GAMMA)a(2.U*GANNA*(EHN**2>«SlN2MA-GAMMA 
1  *1.0)/(GAMMAAt.0) 

PHUSHK  CALCULATES  THE  STATIC  PRESSURE  RISE  PUR  AH  UHLIOUE 
SHOCK  GIVEN  THE  APPROACH  MACH  NUMMF.R,  THE  SINE  SjUARED  UF 
THE  aAVE  ANGLE  ANU  GAMMA  (COM  12b) 

NITaO 

NITMAX  a  1$ 

EHNlaEMNaSFlCMSI ,GAMMAI ) 

OELTIMaOELT AMICMMI ,GAMMA1,SINHA2(EMN| ,GAPMA1 I ) 

PRMAXIa  PHUSHK(CMN|,SIN«A2(EMN|,CAMMA| 1 ,GAMMAI  1 
EMN2  a  EMNMSr(EMS2,GAMMA2) 

OCLT2M  a  OCLTAN(eMH2,GAHMA2,SINNA2(EMN2,GAMMA2)  ) 

PRHAX2  a  PRUSHK(tMN2,SlRNA2(EHN2,GAMMA2),GAMMA2) 

IF< (THETA1-THCTA2)  .GT,  IDELT1MACELT2M ) )  GO  TU  bOO 

UETCRMIRE  HHICH  stream  IS  THE  Hi AK  STREAM 

IFIPRMAXl  •  PRMAX2)tOO,IOO,IIO 

STREAM  I  IS  THE  aEAK  STREAM 

too  PRSSA  a  PRMAXt 

PRl  a  PRSSA 

DELI  a  DLLTIH 

PRNS  a  (%.U*GAMMAt*EMNt**2  •  GAMHAt  *  1 . n  ) / ( GAMMA  I  A  I . D ) 

CO  TU  120 

STREAR  2  IS  THE  aEAKVR  STREAM 

110  PRSSA  a  PRmaX? 

PRl  a  PRSSA 

DELI  a  UELT2M 


PUNS  »  t  2.0*C*M''*2«t>'N^»«2-(.*»>*.A2»|  .<I)/(CA'1»<»2»  I  ,o» 


7  c 

CALCUtATL  TKt  MAXIMUM  ftirAL  STPfcAM  TUHtllUt;  tHGUf.  HY  A 
NUmCKICAL  INTLMACTtUN  PHOCPIXIkC 

t2U  ICUUkT  ■  I 

CALL  SrSHA (PKSSA, LMNl .OAMMAI .OAl.Tl , HOrLTl ) 

CALL  STSHA(PNSSA,t:Hl4  2,GAMMA2,()LLT2,nDrLT2) 

UCLTA  a  LLlTI  *  (>»:LT2 
OTHtlA  *  THRfAl  -  IMElAi 
DULTTA  a  DntLTl  *  r>UtLr2 
riODCL  a  DULTTA 

PHSSH  a  PHSSA  ♦  (PKKS  o’  PHSSA)/2.U 
IF  (UTHETA  -  DELTA)  22V, lb". 130 

MHCN  shuck  tunning  ANGLF.S  FOP  HAX.  DFFLFCTIUN  UF  THK  mFAK 
STREAM  AKF.t  (A)  CHEATER  THAN  THE  STRAMLINE  ANCLES  A  OEAR 
SHUCK  SOLUTION  EXISTS)  (R)  LESS  THAN  THE  STKEAHLINE  ANCLES 
A  STRONG  SOLUTION  SOUGHT. 

I  JO  ICUUNT  alcOUNT  ♦  1 

CALL  STSHK(PHSSB,EHN1,GAMNA1 .DELTl .DDELTi ) 

CALL  STSHk(PRSSH,EMN2,GAMNA2,DELT2,PDELT2) 

F200CL  >  DDELTI  *  0UELT2 
IF  (ABS(F20UEL)  -l.UE«5)  ISO. ISO. 140 
140  PRSSC  a  PRSSB-F20UEL*(PNSSB-PRSSA)/(F2nDEL'*FlUUEL) 

PRSSA  a  PPSSB 
PRSSH  a  PRSSC 
FlUUELa  F2()0EL 
IFdCOUNT  -IS)  130,130,170 
ISO  OELMAX  a  DF.LTl  ♦  DELT2 

THIS  completes  the  INTERACTION  LOUP  TO  DETERMINF.  THE  MAXIMUM 
TURNING  CAPARILITT  OF  THE  THO  STREASMS. 

IF(ObLMAX>DTHETA)bOO, 160,200 
160  thetas  a  THETAl  -  OELTl 
NSTOP  a  I 
RETURN 


170  HR1 TCI J,SOO)1COUNT.F2UOEL,PRSSB 

SOO  FORMAT(*l',SX,'*»*CON»ERCtNCC  ERROR  IN  SLIP.LUUPI  CIJ, 
1  SX.F10.«,SX,FI0.6,  •  )  ••*•} 

IFKDELTl  ♦  ObLT2  )  '•  DTHETA)  600,160,200 

BEGIN  ITERATION  LOOP  FUR  SLIP  LINE  ANGLE 

200  NIT  a  1 

IF(ABS(UUCLTlAnDELT2)  .LE.  l.OE'N)  GO  TU  160 
IF  (DELTA  •  OTHETA)  2IU,21'>.220 
210  FIDEL  a  DELTA  -  UTNFTA 
PR2  ■  PNl  •FIUEL/UDLTTA 
GO  TO  230 
220  NIT  a  1 

PRl  a  1,0  «  07HF.TA«(PRI>I.»)/0FI.1  A 
CALL  STSNkIPRI  ,  EMII  1  .CAmmAI  ,riFL1 1  .DDEITI  ) 

CALL  STSNKIPRI ,CMN2,CAMMA2,DEL12,nDFLT2> 
IFIABSIUdELTIaOUEL12)  .LE.  I.0F.-6)  GO  TU  160 
UDEL  a  DELTI  »  OtLI2  -  UTHFTA 
PH2  a  PR|  .  FIDAL/IUOEL'TI  .  t  lipKLT2) 

230  NIT  a  NIT  *  I 

CALL  STSHK  |PR2,EHM|  .CAmm  a|  ,  UK  IT  I  ,  DPKI.T  I  ) 


■. "V  '  "V  '  TT  V'*  '.T  '. 


C»l  L  SrSHK(»«H2,t<*«2,GAHi.A2,l>tl,t?,OOf.LT2) 
lF(ABS(l>ULl,Tl*PUeL1i)  .LE.  t.UE-6)  i.tl  TO  1»0 
F2neL  a  UFI.TI  ♦  Ot  LT2  •  UTHFTA 
lF(AHS(»20Fb)  -  t.ot- b)2SO,7Mi,24n 
24U  PH2  a  PH2or2l)FL/(UI>FLri  *  l>l<ei.T2| 

IFfNlT  -  N|T4AX)23(J,2iU,2t>M 
250  THFTAS  a  THtTAl  -  Otl.tl 
NSTOP  a  1 
KCTUPN 

c 

c 

2(>0  THFTAS  a  THCTAl  -  ULbTl 
NSrUP  a  I 

■  PITE(3,S10)  n1T,F2I)FL,PH2 

510  lOPaATl*  !•  ,5X,  •  ••♦COHKFPGEMCt  F.HHOH  IN  Rl,lH,  t',13,  2(5X.F|0.b 
1  )  ,  •  )•*••) 

HCTURN 

C 

bOO  NSTOP  a  3 

hRirE(3.70a) 

700  FURHAT(15X,4tH  **4SOLUTtUM  FOR  SLIPMHC  AtCte  DOESN'T  bXlST4*4 
3  //I 
RbTURN 
CND 

SUBROUTINE  8TSHK(PHSS,EMN,GAHNA.PrLT,DDEt.r) 

c 

C  THIS  SUBROUTINE  CALCULATES  THE  UBLIOUE  TURnINC  ANCLE  AND 

C  ITS  OEHlVAriVE  AS  A  FUNCTION  OF  THE  PKESSURE  RATIU 

c  accruss  the  shuck. 

c 

c  PRSS  a  PRESURE  RATIO  ACROSS  THE  SHOCK. 

C  F.MN  a  MACH  NUMHER  UPSTREAM  OF  THE  SHOCK. 

C  GAMMA  a  RATIO  OF  THE  SPECIEIC  HEATS. 

C  CELT  a  STREAMLINE  TUMNING  ANGLE. 

C  DDELT  a  OFRIVATlVE  OF  THE  TURNING  ANGLF 

C 

Aa(2.0«GAHHA4CMN«42«GAMMA4| .O'lGAMMA  * I . 0 ) aPRSS ) / ( ( GAMM A 
1  4l.n)*PKSS4  gamma  •  1.0} 

BaKFRSS  •  1.0)/(GAMMA4EMNaa2>PRSS  tI.0))*SUKT(A) 

DELT  a  ATAN(BI 
C 

c  COMPUTATION  or  the  DEHIVATIVC 

C 

CaA/(GANMA«CMNa«2-PR8S*t.O}  •  (GAMMA  *  I .U)4(PHSS>1 ,0)/(  (  ( 

1  GAMMA  tl.OiaRMSS  4  gamma  •  1.0)a42) 

DDELT  a  ( I .0/ ( I ,04B*a2))  •  ( (GANN AMEMNaaT )/ (GAMM A«Cmn«*2 
1  -PRSS4l.O))4(|.0/SORT(A))aC 
NETURN 
ENU 

function  PRSHKIEMSTAR,  OFLTA,  gamma) 

C 

CXAaaoBLlUUE  SHUCK  FUNCTION  (NEFEHENCt  NACA  R-1135) 

C 

C  THIS  FUNCTION  CALCULATES  THE  STATIC  PRESSURE  HATIJ  ACHUSS  AH 

C  OBLIUUC  SHOCK  (aEAK  SULUTIUN)  GIVEN  THF  APPHUACH  MACh  STAR  ANO 

C  THE  TURNING  ANCLF.  (In  NAUIANS). 

C 

C  4a*VARIABLFS44« 

C 

C  F.MSTAR  a  APPHUACM  hACH  STAN  (M4  a  V/C*). 

C  DELTA  a  TUMNIRG  ANGLE  ('IN  HAOIANb). 

C  gamma  a  RATIU  OF  SPECIFIC  HEATS. 

C  PRSHK  a  FINAL  TO  APPHuACH  STATIC  PPFS5UHF  NAriU. 


B-46 


nnnnnnnnn  nnrjn 


dimension  y(j) 

C**««*E()UAT1UN  COKFFIClKMl  FUNCTIONS. 

CUNSrS  (tMSUDfUFLTA.CAMMA)  a  -(FMSOn  *  2.0)/EMS0U  • 

I  UAMMAtlSIN  (0tLTA)*«2) 

CUNSTC  (LMSOO.DCLTA.dAHMA)  a  (2.0«EMS')r>  *  I  .  0  )  /  ( CMSUU*  •  2  )  ♦ 

I  (((GAMMA  *  I.U|**2)/4.U  *  (GAMMA  -  I .0)/LMSOD)»(S1N  (UELTA)«*2) 
CUNSTU  (i.MSuO,llEt.rA)  a  -(COS  ( Ot I.TA  )  •  » 2  ) / (  F  4S1JD •  * 2 ) 

F.MNSUD  (kMS,GAMMA)a(  ^.U/(GAMMk«t  .U)  )*(eMS**2i/(  1.0 
1  -((GAMM«.| .n)/(GAMMA«1.0) )•(LMS*•2)  ) 

C 

bM2atMNSU0  (EMSfAR, GAMMA) 

C**««*SOLUriUN  OF  CUHIC  bOIIATION  FUR  NAVE  ANCLF  SOUAKEO. 

A  a  (  1  .U/).U)*(  l.U«CUNStC  (C»'2, DELIA, GAMMA)  - 
I  (CONSTB  (tM2, DELTA, GAMMA))**2) 

B  a  (I .a/27.0)«(2.U*(CUNSTU  ( FM2 ,DELTA , G AMM A ) * • J )  - 

1  9.0«(CUNSTB  (f:M2, DELTA, GAMMA))*(CUNSTC  ( EM2 , DELTA , GAMMA ) )  * 

2  2T.O*CUhSTO  (CMI, DELIA)) 

C0SI>H1  a  (  .H/2.0)/SUR1(  •( Aaa )}/3 7 . O ) 

IF(AHS  (COSPHl)  •  1.0)  20,20,10 
10  PHSHK  a  O.u 
NETUHM 
C 

20  PHI  a  (ATAN  (SORTd.O  •  CUSPHJ  *»2  }/CUSPHl  )  ) 

IF  (AND  1,2,2 

1  CHI  a  PHI  ♦  J.14IS9) 

2  OU  )  lal,J 
AI  a  1 


C«a*«ay(n  is  THE  SINE  SUUARED  OF  THE  oAVE  ANGLE. 

J  T(l)  a  2.0«S0NT(-A/).O)*COS  (PHI/).0  ♦  ( A 1 • I .0 ) * 2 . 094 iOS)  - 
1  CUNSTB  (CM2, DELTA, GANMA)/3.0 

••**«THE  ROOTS  OF  THE  CUBIC  EON  alLL  NON  BE  ARRANGED  (N  ASCENDING 


ORDER,  that  is,  1(1)  LESS  THAN  <(2)  LESS  THAN  TO), 

no  «  lal,2 

N  a  I  ♦  1 
DU  i  JaN,} 

IF(Y(1)*T(J))  5,S,4 

4  SAVE  a  T(J) 
r(J)  a  Yd) 

Yd)  a  save 

5  CONTINUE 

6  CONTINUE 

•«*«*THC  ROUT  CORKESPUNUING  TO  THE  NEAA  SOLUTION  IS  Y(2)  AND 

THE  ROUT  Corresponding  to  the  strong  solution  is  yo). 

Yd)  IS  the  square  of  the  sine  or  THE  SHUCK  ANCLE  (SIGMA). 

1  a  2 

PRSHK  a  (3.0«CANMA«EM2«Y(1)  •  (GAMMA  •  l.0))/(GANMA  »  1.0) 

RETURN 

END 

SUBROUTINE  TCGRAL(PHID,CS0D,TRRO,FIIsl,EI ID,C1)J,CI JU.FTAJ, 

I  PHIv),Ct2J) 

«aaa«THlS  SUBROUTINE  CALCULATES  The  TURHULENT  JCI  NIXING  INTEGRALS. 

•  44VAR1  AHLEs'vva 

PHIO  a  DISCRIMINATING  STREAMLINE  VF:t,DCITV  RATIO. 

CSOO.  a  FREE-STREAM  CNOCCO  NUMMFP  SOUARKD. 

TRBU  a  BASF  TO  FREE-STREAM  STAGNATION  TFMPEHAruME  RATIO. 

CIIJ  a  MlAING  INDGRAL  I  FOR  J  .STREAMLINE. 


C  EIID  «  fixing  INrt.Ct<*L  I  »0R  D  .STHLAMM  nE  . 

C  EIJJ  «  mixing  iNlLtiNAI.  3  FtIH  J  STNEAMLlNt. 

C  EIJU  ■  mixing  INTEGMAL  3  F')!*  1>  STRFAMLir.t. 

C  RHll  s  UUMMI  VANlAMLt  lUR  THE  EUHCTIltN  STAfKMtMrS. 

C  TKOI  «  OUM«y  VAHlABLt  FOR  THb  FUrirTIUH  STAfaHMurS. 

C 

C 

UIMENSIUN  TNU(  3!>U)  .Ell  C3&0),F.12I  3M))  »F.(3(3!>0) 

CUMMUN  /EKFVP/  FHl(3bO) 

TJMlFIPHlI.CSQO.rHOll  «  PHI 1 / ( 1  HOI -tS»0» ( PH  1 1 • • 2 ) I 
TJMiFIPHII.CSOn.TRUI)  «  <PMH»*2)/crRlll-CSOO*(PHIl»»2)) 

TJM3F  (PHI  1  ,CSUU, thud  »  (PHII*lRt)l)/(TRUt-CSOO»lPHII**2)) 
C**«**ThE  error  FUNCTIUH  VEtUClTV  profile:,  PHI(I),  is  INlTlALlZbU  IN 
C  EBLUCK  OAfAE  AND  STORED  IN  LAHLLEO  COMnon  EeHryP*.  PHl(l)  IS 

C  GIVEN  FUR  1>I,3SU  VALUES  OF  ETA  IN  THE  RANGE  OF  CTA»3.S  TO 
C  ETAa3.5  IN  INCREMENTS  OF  UETAsO.02. 

C*****INCHCN£HT  size  and  initial  values  at  (ETA  RH)  ARE  SPECIFIED  HERE. 
EPSaEllO 
DETA  a  G.U2 
TKOd)  a  IRBO 
EIKII  a  0.0 
EI2(I)  a  0.0 
Et3(n  a  0.0 

C*****CALCULAT10N  OF  THE  MIXING  TAHLE  BT  THE  TRAPEZOIDAL  RULE. 

DO  2  lal,J44 

THU(1«1I  a  (TRBU  *  ( 1 .O'TKBO t aPhl ( I ♦ I ) 1 

E1I(I«1)  a  Flldi  ♦  U.SVdJNiriPHKI^n.CSOD.TROIUl))  * 

1  TJMIF(PHl(i),CSUO,TRU(l)))*DETA 

E12(l«ll  •  £12(1)  4  0.b«(TJM2F(PHI(l4l ),CSOD.TKO(l«l))  4 
1  TJM2F(PH1(I),CS0D.TH0(1)))*DETA 

E13(14I)  a  E13(I)  4  O.S* (TJm3P(PH1 ( 14 1) ,CSUD.TR0( 1 4 1) )  4 
I  T0M3F(PHl(l),CStiD,THU(l)))*0ETA 

J  a  |4t 

iriPHKJ)  .LT.  (0.25))  GO  TO  2 

ir(ABS(l.O-((ClUO)>F.I2(J))/(Ell(l)-E12(I)>)).LE.I.0C-04>  GO  TO  3 

2  CONTINUE 

C44444DFTCRR1NE  THE  J«  ANU  U*S1HEAML1NE  VALUES  OF  THE  INTEGRALS. 

3  EllJ  a  Ell(J)  -  CI2(J) 

C«««**TABLE  SEARCH  AND  IN  1 CNPULATION  FOR  CI30. 

DO  4  lat.J 

1F(EI1(1)  .GT.  EllJ)  CO  TO  B 

4  CONTINUE 

5  EllJ  a  e13(I<I)  4  ( (El 3( 1 )>EI3( 1-1 ) ) / (El I ( 1 )>C1 I ( 1-1 ) ) )* 

I  (ClIJ-ElKl-D) 

E12J  a  E12(I-t)  4  ((ll2(I)-ei2(l-l))/(Ell(l)-CIl(l-i)))* 

1  (EllJ-EIKl-ll) 

PHIjaPHI(l>t)  4((PH1(I)>PH1(1-I))/(C1I(1)-C1I(1-I)))« 

1  (EIlJ-EIKl-in 

ETAJaFLUA1(l«2)  «  UETA4 ( UETA/ (El U 1 ) -F 1 1 ( I - 1 ) )  ) * 

1  (EllJ-EIKl-l  ))-3.S 

C*«**aTABLE  SEARCH  AND  INTERPOLATION  FOR  EIID,  EI3D. 

00  S  lat.J 

IFIPHKll  .CT.  PHIU)  GO  TO  7 
S  CONTINUE 

7  EIID  a  Elld-n  4  ((Ell(l)-tll()-l))/(PHId)-PHld-|)))* 

1  (PHIO-RHId-D) 

E130  a  ElJd-l)  4  ((ll3(l)-EIld-l))/(PHld)-PHl(T-l)))» 

I  (PNI(t-PHl(  l-l  )  ) 

IF  (ABS(EPS).LT.I.E-b)  RETURN 

CIlDaEllU-FPS 

IF  (EIID)  H,8,N  . 

•  EllDaO. 

PHinaI.F.-t> 


kLTUtiN 

9  DU  lU  isl ,J 

IF  (blUD.GT.tllU)  GUTd  tl 

10  CUNTINUK 

11  PH10«PM1(1-1  J  ♦  ((«-|4l(l)-PHl(l-|))/<i'M(t)>Ell(l>ni)* 

1  (eild<Ell(I-t)) 

PCTUPN 

end 

SUSNUUriME  TfcGPSLCPHillfCliUK.THHU,!  llDfEl  3I>) 


>THIS  SUHNUliriNb  CALCULATES  THE  TUkkUl.FNT  JET  MIXIMC  IMTEGHALS 
KHFN  PHI-0  IS  OEFINEU  bT  THE  UbEHA  RPCUMPRESSIUn  CHITEKIA 


***VARlAbLES**« 

PHlb>DlSCRIH|NArrNG  STKEAPLINE  VEI nciTV  HATIU. 
CSOU«FhEE-STREAM  CRUCCO  NUHHER  SUliANEO. 

THPUabASL  TO  FREL-STRLAH  STAGNATION  TENPF NATURE  RATIO. 
F.IIJSNIXING  INTEGRAL  1  FUR  J  STREANLIUE. 

EllONNIXlNG  INTEGRAL  1  FOR  »  SIRtAbLINE. 

EIJJaMIXlNG  INTEGRAL  3  FUR  J  STREANI.INF. 

ElJUaMIXlNC  integral  J  FOR  0  STREAMLINE. 

PHIIaOUMMI  VARIAHLE  FUR  THE  FHbCTlUN  STATEMENTS. 
TROlaOUMNT  VARIABLE  FUR  THE  FUNCTION  STATEMENTS. 


OIMENSION  1RO(iSO),EII (3SU),E13( 350) 

COmmuN  /LRFVP/  PHI(3S0) 

TJNir (PHIl,CSO0,TMOl)apMll/(THUI-CSi)n»(PHI iaa2)  ) 
TJM3FlPHlI,CSOO,IR01)alPHiiaTHOI3/ITR01-CSOUa(pMIlA*2») 

M«««THC  ERRUR  function  VELUCITT  PROFILE.  PHI(I).  IS  INITIALIZEO  IN 
«HLOCK  OATA*  ANU  STOREO  IN  LARELEO  COMNUN  aCRFVPa.  PHI (I)  IS 
GIVEN  FUR  Ial,3S0  VALUES  OF  ETA  IN  THE  RANGE  UF  ErAa.J.S  TU 
CTAS3.S  IN  increments  UF  OETAao.02 

l*aaa INCREMENT  SIZE  And  INITIAL  VALUES  AT  (ETA  HR)  ARE  SPECIF lEO  HERE. 
OFTAaO.02 
TRO(I)bTRBO 
EIKDaO.U 
El3(l)a0.0 

M«a«CALCULATION  OF  THE  NIXING  TABLE  BY  THF.  TRAPF.ZUlUAL  RULE. 

DU  2.1a|,349 

TROllMUa  (IRBO  ♦  (  1  .  <i-THR0  )  aPHI  1 1  ♦  I  )  ) 

LlKlMUa  tIKI)  ♦  0.5a(TJN|r<M»l(Ial),CSi)i>,THo(ff)))  ♦ 

I  TJM|F(PHl(I),CSOU.TR(Ml)))aOF:TA 

EIKlAl)  a  C13(I)  t  0.5aiTJM3F(P»«I(lFl  ),CSOO,TRU(  )♦!  )  )  ♦ 

1  TJM3F( PHI ( I ) .CSUD.TRUI I ) ) ) aPETA 

IF  (PHKIAI  ).CE.PHIb)  GU  TO  3 

2  CONTINUE 

>***«UbrLRMlNE  THE  (••STREAMLINE  VALUF.S  UF  THE  INTEGRALS 

3  CIIU>EII(I)A((F.lt(ltI)-EII(l))/(PHl(Itl)-PHt(I)))« 

1  (PHI(;-PHI(I)) 

EIIUaEI3(I)«((F:i3(It|)-bI3(l))/(PMI(l*l)-PM|(l)))* 

1  (PHIU-PHId)) 

RF.TURN 

F.ND 

FUNCTION  SIGMAF  (CSUO.TRPU, gamma. »SIGMA) 

•  AtatTHIS  FUNCTION  CALCULATES  THF.  VALUF  UF  THfc  SIMILARln  PARA-^tTER, 
SIGRA.  FROM  THE  C'JRRFLATIUNS  UF  KURST  AND  TRIPP  :i>«  CHA'IMAPHAGADA 

NSICNAa  O.Sir.MA  AFTER  AOP&T  AN|  TRIPP 
NSIGNAb  I. sigma  after  CHANNAPRAI-AI'A 
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If  (NStG''«.Nf .0)  GUI II  I 

KOBSr  *N0  TRIHP  •»•«*»#•••**»••»••••»•*••**•»»♦•♦»*••»•••*»*•« 

SIGPAra  r.;,0«^.7bii«2.0*CSOI>/(  (  1 .0-CSO0)»(r,A»iM*-l  .0)  ) 

KtTURN 

c 

channahnagau*  •«*»♦»»•**•»•»•«**»•••»»•»»•*•»•*»•»•*»•*♦•••*•* 

1  ir  (C500-U.7)  2,J.i 

2  H«0.25 

cmo  4 

J  H«0.S*(CSQn-0.7)*0.75 

4  SlOaA^  a  l2.U/(R*(1.0*(t.0-CSOI<)/TRHI))) 

RfTURN 

KRO 

BLOCK  DATA 

*«*»*TNf  CHROR  fUNCrtnil  VCLUCITV  PRUriLC»  PHI(I),  IS  HUTlALI^bU  IR 
•LLUCK  UAIA*  and  STUHEO  IN  LAHfLbl)  CUMMnN  KbHFVP*.  PHl(I)  IS 
CIVfN  rOR  lal.JSO  VALUES  Of  fTA  IN  THf  NANCE.  Of  f.TAa-J.S  TO 
EIAaJ.S  1.V  InCRCNENTS  UF  UE:TAa0.n2. 

CUMHUN  /ERFVP/  At(4Si.A2(4S).A)(4S),A4(4S),AS(tSI,Ab(4S),i47(4S), 


1  A80)> 

UAIA  At 


/O. 000000 

O.OOOOOO 

• 

O.OOOOOO 

0 

0.01)0000 

0 

0,000000 

* 

0.000000 

0 

o.oonool 

0 

0.000001 

0 

O.IKlOOOl 

0 

0.000001 

0.000001 

• 

o.gotiooi 

0 

0.000002 

0 

0.1)110002 

0 

0.000O02 

0.00000 J 

0 

o.oouool 

0 

0.000004 

0 

0.000004 

0 

0.000005 

0 

o.oonoos 

0 

O.UUOOON 

0 

0.000007 

0 

0,000009 

0 

0.000009 

0. 000011 

0 

0.000012 

0 

0.000014 

0 

0.000016 

0 

0.000016 

0 

0.000020 

0 

0.000021 

0 

0.000026 

0 

0.000029 

0 

0.000011 

0 

0.000037 

0 

0.000U42 

0 

0.000047 

0 

O.OOOOSl 

0 

0,000059 

0 

0.000067 

0 

0.000075 

0 

0. 000094 

0 

0.000094 

0 

0. 000109 

/ 

DATA  A2 

/O.OOOllS 

0 

0.000111 

0 

0.000147 

0 

0.000164 

0 

0.000192 

0 

0,000201 

0 

0,000226 

0 

0.000291 

0 

0.000279 

0 

0.000110 

0 

0.000344 

0 

0.000101 

0 

0.000422 

0 

0.000467 

0 

0,000517 

0 

0.000971 

0 

0.0006 11 

0 

0,000696 

0 

0.000767 

0 

0.000849 

0 

0.000911 

0 

0.001024 

0 

0.001126 

0 

0.001217 

0 

0.001198 

0 

0,001409 

0 

0.001612 

0 

0.001799 

0 

0.001996 

0 

0.002140 

0.002330 

0 

0.002991 

0 

0.002796 

0 

0.001019 

0 

0.001110 

0 

0.001604 

0 

0.001921 

0 

0.004261 

0 

0.004611 

0 

0.009027 

0 

0.009494 

0 

0,009912 

0 

0.006404 

0 

0.006912 

0 

0.007496 

/ 

DATA  AJ 

/O. 000104 

0 

0.000791 

0 

0.009446 

0 

0.010169 

0.010960 

0 

0.011029 

0 

0.012729 

0 

0.011699 

0 

0.014706 

0 

0.015792 

0.016946 

0 

0.010172 

0 

0.019472 

0 

0.029951 

0 

0.022111 

0.0230ST 

0 

0.029491 

0 

0.027219 

0 

0.079041 

0 

0. 010967 

0 

0.01299b 

0 

0.019111 

0 

0.017192 

0 

0.019747 

0 

0.042211 

0 

0.044041 

0 

0.0479P2 

0 

0.090491 

0 

0.091460 

0 

0.096607 

0 

0.099099 

0 

0.061119 

0 

0.066910 

0 

0.070677 

t 

0.074561 

0.070692 

0 

0.092697 

0 

0.097291 

0 

0.091969 

0 

0.096620 

0 

0.101990 

DATA  A4 

0 

0.106661 

0 

0.111999 

0 

0.117434 

0 

0.121101 

/ 

/O. 120996 

0 

0. 1  19002 

f 

0.141219 

0 

0.147669 

0 

0.154292 

* 

0.161100 

0 

0,169119 

0 

0.179322 

0 

0.IR2/16 

0 

0.190105 

0.190004 

0 

0.206091 

0 

0.214705 

0 

0.272544 

0 

0,2  11069 

0.239769 

0 

0.249641 

0 

0.297666 

0 

0.266904 

0 

0.27626  1 

0.209022 

0, 

0.299914 

0 

0.109154 

0 

0.115116 

0 

0.125457 

0.319709 

1 

• 

n. 146062 

0 

0.196972 

0 

0.167171 

0 

0.  17797*1 

0.100673 

0 

0,199997 

0 

0.410519 

0 

0.421552 

0 

0.412647 

0.441799 

0 

0.494999 

0 

0.466217 

0 

0.477472 

0 

0.466746 

0.900029 

0 

0.911111 

0 

V. 922599 

0 

0.911640 

0 

0.545069 

/ 

DATA  AS 

• 

/O. 996261 

0 

0.967409 

0 

0.976504 

0 

<1. 96651b 

0 

0  ,  *,0O4«>> 
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O.bl 1  )M2 

0.672179 

9 

». 632061 

9 

0.64  3400 

0.653911 

l>.6(<4  3  44 

0.6  14593 

0.604717 

9 

0.694695 

0.  104534 

0.714226 

0.723763 

9 

0.7331*1 

9 

0. 742 356 

0 

0. 751 40J 

0. 760276 

9 

0.760977 

9 

0.777497 

9 

0.705434 

0 

0.793986 

O.Hl;l«S4 

9 

0.6097  31 

9 

0.617317 

9 

0.074712 

9 

0.631915 

0.636923 

0.045739 

9 

0.057361 

9 

0.050789 

9 

0.865026 

0.97107U 

0.076925 

9 

9 

0. *00068 

9 

0.693361 

0.690471 

0.903400 

9 

0.900151 

9 

0.912126 

9 

0.917130 

/ 

DAT*  *6 

/n. 921 364 

0.925437 

9 

0.92«337 

9 

0.9 3  3063 

9 

0.936614 

9 

0.940113 

9 

0.943404 

9 

0.946550 

9 

0.949557 

0 

0.952427 

9 

n.9b5l6b 

9 

a.s'tm* 

9 

0 

O.U62674 

9 

0,964871 

9 

0.967009 

0.969037 

9 

0.970961 

9 

0.972785 

9 

0.97*51 1 

9 

0.976146 

9 

0.977691 

9 

0,979151 

0 

0.900529 

0 

0.9NI 829 

# 

0.9630S4 

9 

0.984208 

9 

0.965293 

0 

0.906314 

0 

0.987214 

9 

0.906174 

0.969018 

9 

0.909610 

0 

0.99055! 

9 

0.991245 

9 

0.991094 

0.992500 

9 

0,99J065 

0 

0.993593 

9 

0.994065 

9 

0.994543 

9 

0.994969 

9 

0.995365 

0 

0.995733 

9 

0.99607b 

/ 

DATA  *7 

/O. 996392 

0.996606 

9 

0.996958 

0 

0.991210 

9 

0.997442 

0.997657 

9 

0.997H56 

9 

0,99*039 

0 

0.998206 

9 

0.99*363 

0,990506 

9 

0.998636 

9 

0.99*758 

0 

0.999869 

9 

0,996971 

0.999064 

0.999149 

9 

0.999227 

0 

0.999299 

9 

0.999364 

9 

0.999424 

0.999478 

9 

0.999527 

0 

0.999577 

0 

0.999613 

0.999651 

0.999685 

9 

0.999715 

0 

0.999743 

0 

0,999766 

9 

0.999791 

9 

0.999812 

9 

0.999831 

0 

0.999848 

9 

0.999663 

§ 

0.999677 

0.999889 

9 

0.999900 

0 

0.999910 

9 

0.999919 

9 

0,999927 

9 

0.9999Jb 

9 

0.999941 

0 

0.999947 

9 

0. 999952 

f 

DATA  A8 

/O. 999957 

0.999961 

9 

0.999965 

0 

0.999968 

9 

0.999971 

0.999974 

0.999976 

9 

0,999976 

9 

0.999960 

9 

0.999982 
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